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1. Introduction 



The goal of the lectures on lattice QCD (LQCD) is to provide an overview 
of both the technical issues and the progress made so far in obtaining 
phenomenologically useful numbers. The lectures consist of three parts. 
My charter is to provide an introduction to LQCD and outline the scope 
of LQCD calculations. In the second set of lectures, Guido MartincUi will 
discuss the progress we have made so far in obtaining results, and their 
impact on Standard Model phenomenology. Finally, Martin Liischer will 
discuss the topical subjects of chiral symmetry, improved formulation of 
lattice QCD, and the impact these improvements will have on the quality 
of results expected from the next generation of simulations. 

QCD is the regnant theory of strong interactions. It is formulated in 
terms of quarks and gluons which we believe are the basic degrees of free- 
dom that make up hadronic matter. It has been very successful in pre- 
dicting phenomena involving large momentum transfer. In this regime the 
coupling constant is small and perturbation theory becomes a reliable tool. 
On the other hand, at the scale of the hadronic world, ^ 1 GeV, the 
coupling constant is of order unity and perturbativc methods fail. In this 
domain lattice QCD provides a non-perturbative tool for calculating the 
hadronic spectrum and the matrix elements of any operator within these 
hadronic states from first principles. LQCD can also be used to address 
issues like the mechanism for confinement and chiral symmetry breaking, 
the role of topology, and the equilibrium properties of QCD at finite tem- 
perature. Unfortunately, these latter topics were not covered at this school, 
so I will give appropriate references to compensate for their omission. 

Lattice QCD is QCD formulated on a discrete Euclidean space time grid. 
Since no new parameters or field variables are introduced in this discretiza- 
tion, LQCD retains the fundamental character of QCD. Lattice QCD can 
serve two purposes. First, the discrete space-time lattice acts as a non- 
perturbative regularization scheme. At finite values of the lattice spacing 
a, which provides an ultraviolet cutoff at n/a, there are no infinities. Fur- 
thermore, renormalized physical quantities have a finite well behaved limit 
as a — > 0. Thus, in principle, one could do all the standard perturbative 
calculations using lattice regularization, however, these calculations are far 
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more complicated and have no advantage over those done in a continuum 
scheme. Second, the pre-eminent use of transcribing QCD on to a space- 
time lattice is that LQCD can be simulated on the computer using methods 
analogous to those used for Statistical Mechanics systems. (A brief review 
of the connection between Euclidean field theory and Statistical Mechanics 
is given in Section 5.) These simulations allow us to calculate correlation 
functions of hadronic operators and matrix elements of any operator be- 
tween hadronic states in terms of the fundamental quark and gluon degrees 
of freedom. 

The only tunable input parameters in these simulations arc the strong 
coupling constant and the bare masses of the quarks. Our belief is that 
these parameters are prescribed by some yet more fundamental underlying 
theory, however, within the context of the standard model they have to 
be fixed in terms of an equal number of experimental quantities. This is 
what is done in LQCD. Thereafter all predictions of LQCD have to match 
experimental data if QCD is the correct theory of strong interactions. 

A very useful feature of LQCD is that one can dial the input parameters. 
Therefore, in addition to testing QCD we can make detailed predictions of 
the dependence of quantities on as and the quark masses. These predictions 
can then be used to constrain effective theories like chiral perturbation 
theory, heavy quark effective theory, and various phenomenological models. 

My first lecture will be devoted to an overview of the scope of LQCD and 
to showing that simulations of LQCD are a step by step implementation of 
field theory. The second lecture will be devoted to explaining the details of 
how to transcribe the quark and gluon degrees of freedom on to the lattice, 
and to construct an action that, in the limit of zero lattice spacing, gives 
continuum QCD. I will also spend some time on issues of gauge invariance, 
chiral symmetry, fcrmion doubling problem, designing improved actions, 
the measure of integration, and gauge fixing. 

Numerical simulations of LQCD are based on a Monte Carlo integration 
of the Euclidean path integral, consequently, the measurements have statis- 
tical errors in addition to the systematic errors due to lattice discretization. 
In order to judge the quality of lattice calculations it is important to under- 
stand the origin of these errors, what is being done to quantify them, and 
finally what will it take to achieve results with a given precision. These is- 
sues will be covered in the third lecture along with an elementary discussion 
of Monte Carlo methods. 

The fourth lecture will be devoted to the most basic applications of 
LQCD - the calculation of the hadron spectrum and the extraction of 
quark masses and a^. Progress in LQCD has required a combination of 
improvements in formulation, numerical techniques, and in computer tech- 
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nology. My overall message is that current LQCD results, obtained from 
simulations at a ~ 0.05 — 0.1 fermi and with quark masses roughly equal 
to TOg, have already made an impact on Standard Model phenomenology. 
I hope that the three sets of lectures succeed in communicating the excite- 
ment of the practitioners. 

In a short course like this I can only give a very brief description of 
LQCD. In preparing these lectures I have foimd the following books [1-4], 
and reviews [5-8] very useful. In addition, the proceedings of the yearly 
conferences on lattice field theory [9] are excellent sources of information 
on current activity and on the status of results. I hope that the lectures at 
this school and these references will allow any interested reader to master 
the subject. 

2. Standard Model of Particle Interactions 

The Standard Model (SM) of particle interactions is a synthesis of three 
of the four forces of nature. Thc^so forces are described by gauge theories, 
each of which is characterized by a coupling constant as listed below. 

STRONG INTERACTIONS ~ 1 

ELECTROMAGNETIC INTERACTIONS Uem « 1/137 

WEAK INTERACTIONS Gp « 10"^ GeV-^. 

The basic constituents of matter are the six quarks, u, d, s, c, b, t, 
each of which comes in 3 colors, and the six leptons e, Ve, H, v^, t, Vr- 
The quarks and leptons are classified into 3 generations of families. The 
interactions between these particles is mediated by vector bosons: the 8 
gluons mediate strong interactions, the and Z mediate weak interac- 
tions, and the electromagnetic interactions are carried by the photon 7. 
The weak bosons acquire a mass through the Higgs mechanism, and in 
the minimal SM formulation there should exist one neutral Higgs particle. 
This has not yet been seen experimentally. The mathematical structure 
that describes the interactions of these constituents is a local gauge field 
theory with the gauge group SU{?>) X SU{2) x U{1). The parameters that 
characterize the SM are listed below; their origins are in an even more 
fundamental but as yet undiscovered theory. 
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Pareimeters in the 

Parameters Number 



Standard Model 

Comments 



Masses of quarks 



6 u, d, s light 

c, b heavy 
i = 175 ± 6 GeV 



Masses of leptons 



Mass of 

Mass of Z 

Mass of gluons, 7 

Mass of Higgs 

Couphng ag 

Couphng ttem 

Couphng Gp 

Weak Mixing Angles 
CP Violating phase 

Strong CP parameter 



e, fj,, T 



80.3 GeV 
91.2 GeV 

(Gauge symmetry) 
Not yet seen 

« 1 for energies < 1 GeV 
1/137 (=1/128.9 at Mz) 

10-5 

^ij ^2, ^3 
5 

e =? 



The question marks indicate that these values are tiny and not yet mea- 
sured. No mixing angles have been put down for the lepton sector as I have 
assumed that the neutrino masses arc all zero. A status report on exper- 
imental searches for neutrino masses and their mixing has been presented 
by Prof. Sciulli at this school. The structure of weak interactions and the 
origin of the weak mixing angles has been reviewed by Daniel Treille, while 
David Kosower covered perturbative QCD. 

In the SM, the charged current interactions of the W-bosons with the 
quarks are given by 

nw = \{ J^W^ + h.c.) (2.1) 
where the current is 

■/m = ^ («. 2^)741-75)1^ I s I (2.2) 
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and the Fermi coupling constant is related to the SU{2) coupling as 
Gf/V2 = gl/8M^. V is the 3x3 Cabibbo-Kobayashi-Maskawa {CRM) 
matrix that has a simple representation in terms of the flavor transforma- 
tion matrix 



This matrix has off-diagonal entries because the eigenstates of weak inter- 
actions are not the same as those of strong interactions, i.e. the d, s, b 
quarks (or equivalently the u, c, t quarks) mix under weak interactions. 

For 3 generations the unitarity condition = imposes constraints. 
As a result there are only four independent parameters that can be ex- 
pressed in terms of 4 angles, 9i,92, 6*3 and the CP violating phase 5. The 
representation originally proposed by Kobayashi and Masakawa is 



where — cos 9i and Si — sin 6i for i = 1, 2, 3. A phcnomenologically more 
useful representation of this matrix is the Wolfenstein parameterization 
(correct to 0(A'')) 



where A = sin 6'c = 0.22 is known to 1% accuracy, A w 0.83 is known 
only to 10% accuracy, and p and rj are poorly known. The elements Vtd 
and Vub of the CKM matrix given in Eq. 2.5 are complex if 77 ^ 0, and 
as a result there exists a natural mechanism for CP-violation in the SM. 
The phenomenology of the CKM matrix and the determinations of the four 
parameters have been reviewed by Profs. Buras and Richman at this school 



A list of processes that are amongst the most sensitive probes of the 
CKM parameters is shown in Fig. 1. Lattice QCD calculations are be- 
ginning to provide amongst the most reliable estimates of the hadronic 
matrix elements relevant to these processes. Details of these calculations 
and their impact on the determination of p and r] parameters will be cov- 
ered by Prof. Martinelli at this school. My goal is to provide you with the 
necessary background. 




(2.3) 




(2.4) 




[10,11]. 
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CKM 



D Klv K ^ttIu B ^ -kIv 



D K*lv 



B plu 



I-AV2 ^ A 
-A 1 - AV2 

AX^{1 - p-iii) -AX'^ 

/ /\ 



I 



A\-^{p - ir]) \ 




Mixing 



Mixing ( ^^'^ \ B^D*l. 
{BbJ V Decays; 



Fig. 1. The CKM matrix in the Wolfenstein parameterization. I show examples of phys- 
ical processes that need to be measured experimentally and compared with theoretical 
predictions via the "master" equation to estimate various elements of the CKM matrix. 



3. Scope of Lattice QCD Calculations 

The primary goal of LQCD is to test whetiier QCD is tlie correct tiieory 
of strong interactions. Most particle physicists believe that the answer is 
yes because of its successes in describing processes with large momentum 
transfer, its mathematical elegance, and uniqueness. However, quantitative 
confirmation is lacking. The problem with the conventional perturbative 
approach (expansion in a small parameter) when analyzing hadronic pro- 
cess at scales < 1 GeV is that the strong coupling constant ~ 1. Thus, 
perturbation theory in as is not reliable. As a result we cannot calculate 
the masses of mesons and baryons from QCD even if we are given as and 
the masses of quarks. This is illustrated in Fig. 2 where I show a com- 
parison between the binding energy of a hydrogen atom and that of the 
proton. It turns out that almost all the mass of the proton is attributed to 
the strong non-linear interactions of the gluons. 

Simulations of lattice QCD have six unknown input parameters. These 
are the coupling constant as and the masses of the up, down, strange, charm 
and bottom quarks (the top quark is too short lived, 0.4 x lO"^"* seconds, 
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QED 




Hydrogen Atom 



Me =0.5MeV 
Mp =938MeV 
Ending = 13.6 eV 
(EM force) 



M u ~ 3 MeV 

M d ~ 6 MeV 
M p =938 MeV 
(Strong force) 

Fig. 2. An illustration of the difference in the binding energy between the electron and 
the proton in a hydrogen atom (interacting via electromagnetic forces) and the quarks 
in a proton (strong force). 

to form bound states and is thus not suitable for lattice studies). Once 
these have been fixed in terms of six precisely measured masses of hadrons, 
thereafter, the masses and properties of the hundreds of other particles 
made up of these quarks and gluons have to agree with experiments. For 
example, consider the chiral perturbation theory (xPT) relation for the 
pion mass 

{M+f = B^{mu + nid) + C^(to„ + nidf + .... (3.1) 

Lattice calculations will provide a check of this relation and fix the con- 
stants Bjr, Ctt, . . ., i.e. constrain the chiral lagrangian. In addition, having 
measured M^r for a number of values of the quark mass, one can invert 
Eq. 3.1 to determine the current quark masses. I shall call such tests of 
QCD and the determination of the six unknown parameters of the SM a 
direct application of lattice QCD. 




Proton 
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Fig. 3. The Feynman diagram for the semi-leptonic decay D~ — > K e~V^. The QCD 
corrections are illustrated by the various gluons being exchanged between the initial and 
final hadrons. The leptonic vertex can be calculated reliably using perturbation theory, 
whereas the hadronic vertex requires non-perturbative methods. 

The second area in which perturbative estimates are not reUable is in 
the calculation of the matrix elements occurring in the weak decays of 
hadrons. The reason is that the non-perturbative QCD corrections to the 
basic weak process can be large due to the exchange of soft gluons between 
the initial and final states. This is illustrated in Fig. 3 for the case of 
the semi-leptonic decay D Klv. In this case the matrix element of 
the weak interactions Hamiltonian between the initial D meson and final 
kaon receives large corrections from QCD which cannot be estimated using 
PQCD. 

The result of the experiment is a number, i.e. the decay rate. To check 
whether the Standard Model prediction is in agreement with this number 
theorists must derive an analytical expression for the same process. This 
expression, in general, consists of a product of three parts: a combination 
of the parameters in the S'M, matrix elements (ME) of the interaction 
Hamiltonian between the initial and final states, and kinematic factors. 
(This is the well-known Fermi's golden rule.) Schematically one can write 
this "master" equation as follows. 

Expt.^ — (SM parameters) (matrix elements) (kinematic factors). (3.2) 

Thus, for each such experimental number one gets a constraint on one 
particular combination of the SM parameters provided the AIE are known. 
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By using several precision experiments and calculating the necessary M E 
one can over-constrain the model and therefore test it. I shall call lattice 
calculations of ME an "indirect" probe of Standard Model parameters 
as one usually needs additional phenomenological/experimental input to 
extract the unknown parameters via the master equation. 

Physicists have been busy calculating and measuring these parameters 
for the last 25 years, so one already has reasonable estimates. It is, there- 
fore, useful to keep in mind the level of precision that lattice calculations 
have to aim for to make an impact. Non-lattice estimates, taken from the 
1996 Review of Particle Properties and from Prof. Buras's lectures [10], 
are 

Parameters Value Comments 



mu(MS,2GeV) 
mrf(MS,2 GeV) 
m,(M5,2 GeV) 
mc{MS, nic) 
mb{MS,mb) 



3-4.5 MeV 
5-8 McV 
100-140 MeV 

1.0- 1.6 GeV 

4.1- 4.6 GeV 



XPT, sum-rules 

XPT, sum-rules 
XPT, sum-rules 
J ftp sprectra 
T spectra 



M-Higgs 

as{Mz) 

\Vud\ 
\Vus\ 
\Vub\ 

\Vcb\ 



0.118 ±0.003 

0.9736 ±0.0010 

0.2196 ±0.0023 

(3.3 ±0.2 ±0.4 ±0.7)10-3 
0.040 ± 0.003 



Not Found 



World average 
n ^ peve and 



K+ 
B - 



— > Tr'^e+i'e and 

ttIv; B ply 
Xclv and 
D*W 



\e\ 

We'\ 



2.26 X 10-3 
- 30 X 10-^ 



K - 



-> mixing 
TTTT decays 



Note that only the magnitude of the CKM elements can be determined from 

experiments. Clearly, the most enticing opportunity for lattice QCD is to 
help determine the quark masses, ag, and the essentially unknown Wolfen- 
stein parameters p and r]. For other elements of Vckm one has reasonable 
estimates, however here too further progress requires significant improve- 
ment in the estimates of the hadronic matrix elements. Phenomenologists 
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are relying on LQCD to provide improved estimates for these. 

Let me end this section by mentioning that in addition to vaUdating 

QCD as the correct theory of strong interactions by reproducing the mea- 
sured spectrum of mesons, baryons, and gluebaUs, and the calculation of 
weak matrix elements, simulations of LQCD will also allow us to 

— Investigate the topological structure of the QCD vacuum, and the mech- 
anisms responsible for confinement and spontaneously broken chiral sym- 
metry. For a recent review see [12]. 

— Calculate hadronic properties like wavefunctions, decay constants, form 
factors for electromagnetic, weak, and rare decays, and structure functions 
for deep inelastic scattering. Estimate the pion nucleon sigma term and 
the proton spin carried by the quarks by calculating the matrix elements 
< N\qq\N > and < N\q^f^^5q\N >. Some of these topics will be covered by 
Prof. Martinelli, for others a good reference is the proceedings of LATTICE 
conferences [9]. 

— Analyze the behavior of QCD at high temperatures and address questions 
like the nature of the phase transition, the transition temperature Tc, the 
equation of state, and the collective excitations near and above Tc- The 
topic of QCD at finite temperature has been covered by J-P Blaizot at 
this school, so I will only list for you my favorite pedagogical reviews of 
the subject [13]. The lectures by Profs. Martinelli, Liischer, and me will 
be confined to understanding and probing the properties of QCD at zero 
temperature. 



4. Overview of the Lattice Approach 

LQCD calculations arc a non-pcrturbativc implementation of field theory 
using the Feynman path integral approach. The calculations proceed ex- 
actly as if the field theory was being solved analytically had we the ability 
to do the calculations. The starting point is the partition function in Eu- 
clidean space-time 




(4.1) 



where S is the QCD action 




(4.2) 
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and M is the Dirac operator. The fermions are represented by Grassniann 
variables ip and ip. These can be integrated out exactly with the result 

Z = jvA^ detM e!"'^ i-k^.-P"") , (4.3) 

The fermionic contribution is now contained in the highly non-local term 
detM, and the partition function is an integral over only background 
gauge configurations. One can write the action, after integration over the 
fermions, as 5 = Sgauge + Sguarks = J ct^x (jF^^Ff^") - J2^ log(DetMi) 
where the sum is over the quark flavors, distinguished by the value of the 
bare quark mass. 

It is expedient to define the "quenched" approximation (QQCD) from 
the very start due to the central role it has played in simulations. It con- 
sists of setting DetM = constant which corresponds to removing vacuum 
polarization effects from the QCD vacuum. Details and some consequences 
of this approximation are discussed in Section 16.8. 

Results for physical observables are obtained by calculating expectation 
values 

P) = VA^ O e-^ . (4.4) 

where O is any given combination of operators expressed in terms of time- 
ordered products of gauge and quark fields. (Expectation values in the 
path integral approach correspond to time-ordered correlation functions.) 
The quarks fields in O are, in practice, re-expressed in terms of quark 
propagators using Wick's theorem for contracting fields. In this way all 
dependence on quarks as dynamical fields is removed. The basic building 
block for fermionic quantities is the Feynman propagator, 

SF{y,3,h;x,%,a) = {M-^f^, (4.5) 

where M~^ is the inverse of the Dirac operator calculated on a given back- 
ground field. A given element of this matrix (M~^)^''^'^ is the amplitude for 
the propagation of a quark from site x with spin-color i, a to site-spin-color 

So far all of the above is standard field theory. The problem we face 

in QCD is how to actually calculate these expectation values and how to 
extract physical observables from these. I will illustrate the second part 
first by using as an example the mass and decay constant of the pion. 

Consider the 2-point correlation function, {0\T[J2^O f{xJ,)O,,{o,0)]\0) 
with t > 0, where the operators O are chosen to be the fourth component of 
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the axial current O f = Oi ^ — ^7475'0 as these have a large coupling to 
the pion. This 2-point correlation function gives the amplitude for creating 
a state with the quantum numbers of the pion out of the vacuum at space- 
time point by the "source" operator Of, the evolution of this state to 
the point {x,t) via the QCD Hamiltonian (or more precisely the transfer 
matrix to be discussed in Section 5); and finally the annihilation by the 
"sink" operator Of at {x,t). The rules of quantum mechanics tell us that 
Oi will create a state that is a linear combination of all possible eigenstates 
of the Hamiltonian that have the same quantum numbers as the pion, i.e. 
the pion, radial excitations of the pion, three pions in J = state, .... The 
second rule is that on propagating for Euclidean time t, a given eigenstate 
with energy E picks up a weight e~^*. Thus, the 2-point function can be 
written in terms of a sum over all possible intermediate states 

(O|E0/(a^,*)^.(O)|O) =E^^^^^^^ e-"*. (4.6) 

X n ^ 

To study the properties of the pion at rest we need to isolate this state from 
the sum over n. To do this, the first simplification is to use the Fourier 
projection as it restricts the sum over states to just zero-momentum 
states, so En Mn. (Note that it is sufficient to make the Fourier pro- 
jection over either Oi or Of.) The second step to isolate the pion, i.e. 
project in the energy, consists of a combination of two strategies. One, 
make a clever choice of the operators O to limit the sum over states to a 
single state (the ideal choice is to set O equal to the quantum mechanical 
wave- functional of the pion), and two, examine the large t behavior of the 
2-point function where only the contribution of the lowest energy state that 
couples to Oi is significant due to the exponential damping. Then 

X ^ 

The right hand side is now a function of the two quantities we want since 
(01^4 (p = 0)|7r) = M^f„. In this way, the mass and the decay constant 
are extracted from the rate of exponential fall-off in time and from the 
amplitude. 

Let me now illustrate how the left hand side of Eq. 4.7 is expressed in 
terms of the two basic quantities we control in the path integral - the gauge 
fields and the quark propagator. Using Wick contractions, the correlation 
function can be written in terms of a product of two quark propagators 
Sp, 



-{0\J2i'{x,th475i^{x,t) ^-(0, 0)7475^(0,0)10) 
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(A) Local Interpolating operators 



(( Gauge 
Yj Links 

Sf (0,0) 
(B) Extended Interpolating operators 

Fig. 4. A scliematic of tlie pion 2-point correlation function for (A) local and (B) non- 
local interpolating operators. 

= (0|^5F(0;(r,t)74755f (a^,i;0)7475|0). (4.8) 

X 

This correlation function is shown in Fig. 4A. To illustrate this Wick con- 
traction procedure further consider using gauge invariant non-local oper- 
ators, for example take O = ■0(a;, t)7475(7^e^= '^^^^^''^^^)'tp{y,t) where V 
stands for path-ordered. After Wick contraction the correlation function 
reads 

(0| Sf{0; X, i)7475(^ei'r '<^^-)Sf{z, t; y, 0)7475 (^e/o" '3^-)|0). (4.9) 

and involves both the gauge fields and quark propagators. This correla- 
tion function would have the same long t behavior as shown in Eq. 4.6, 
however, the amplitude will be different and consequently its relation to 
/tt will no longer be simple. The idea behind strategies for improving 
the projection of O on to the pion is to construct a suitable combina- 
tion of such operators that approximates the pion wave-function. For ex- 
ample, if 4>{x) is the pion wave-functional, then the "smeared" operator 
O = '0(0)7475 (/ d^x(j){x)^lj{x)) is the ideal choice. 

Finally, it is important to note that to get the 2-point function corre- 
sponding to the propagation of the physical pion we have to average the 
correlation function over gauge configurations as defined in Eq. 4.4, i.e. do 
the functional integral. 
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An overview of how such calculations of expectation values can be done 
consists of the following steps. 

• Defining a finite dimensional system: The Yang-Mills action for gauge 

fields and the Dirac operator for fcrmions have to bo transcribed onto a dis- 
crete space-time lattice in such a way as to preserve all the key properties 
of QCD - gauge invariance, chiral symmetry, topology, and a one-to-one 
relation between continuum and lattice fields. This step is the most diffi- 
cult, and even today we do not have a really satisfactory lattice formulation 
that is chirally symmetric in the = limit and preserves a one-to-one 
relation between continuum and lattice fields, i.e. no doublcrs. In fact, the 
Nielson-Ninomiya theorem states that for a translationally invariant, local, 
hermitian formulation of the lattice theory one cannot simultaneously have 
chiral symmetry and no doublcrs [14]. One important consequence of this 
theorem is that, in spite of tremendous effort, there is no practical formu- 
lation of chiral fermions on the lattice. For a review of the problems and 
attempts to solve them see [15-17] and [18]. 

• Discretization Errors: This problem is encountered when approximating 
derivatives in the action by finite differences, 

dx4>{x) Aa;(/)(x) = + a) - 0(a; - a)) 

= -^(e"^- -e-"^-)<^(a;). (4.10) 

As is obvious this introduces discretization errors proportional to the square 

of the lattice spacing a. These errors can be reduced by either using higher 
order difference schemes with coefficients adjusted to take into account 
effects of renormalization, or equivalently, by adding appropriate combina- 
tions of irrelevant operators to the action that cancel the errors order by or- 
der in a. The various approaches to improving the gauge and Dirac actions 
are discussed in Sections 12 and 13. In my lectures I will concentrate on 
the three most frequently used discretizations of the Dirac action - Wilson 
[19], Sheikholeslami-Wohlert (clover) [20], and staggered [21], which have 
errors of 0(a), 0{asa) — 0{a'^) depending on the value of the coefficient 
of the clover term, and 0{a^) respectively. Reducing the discretization er- 
rors by improving the action and operators increases the reliability of the 
results, and for dynamical simulations may even prove necessary due to 
limitations of computer power. At the same time it is important to note 
that there exists a well defined procedure for obtaining continuum results 
with even the simplest formulation, Wilson's original gauge and fermion 
action [19]. In this case, however, it is necessary to quantify and remove 
the discretization errors by a reliable extrapolation to a = 0. 
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• Generating background gauge fields: The Euclidean action S = 
J d'^x{ jF^^F^^ — Tr log M) for QCD at zero chemical potential is real 
and bounded from below. Thus e"*^ in the path integral is analogous to 
the Boltzmann factor in the partition function for statistical mechanics 
systems, i.e. it can be regarded as a probability weight for generating con- 
figurations. Since S is an extensive quantity the range of S is enormous, 
however, the distribution is very highly peaked about configurations that 
minimize the free energy. Configurations close to the peak dominate the 
functional integral. Such important configurations can be generated by 
setting up a Markov chain in exact analogy to, say, simulations of the Lsing 
model. They are called "importance sampled" as they are generated with 
probability of occurrence given by the weight e"*^. For a discussion of the 
methods used to update the configurations see [1] or the lectures by Creutz 
and Sokal in [2]. 

• Calculation of the quark propagator on a given background gauge field: For 
a given background gauge configuration, the Feynman quark propagator Sp 
is a matrix labeled by three indices - site, spin, and color. A given element 
of this matrix gives the amplitude for the propagation of a quark with 
some spin, color, and space-time point to another space-time point, spin, 
and color. Operationally, Sp is simply the inverse of the Dirac operator 
M . Once space-time is made discrete and finite, the Dirac matrix is also 
finite and its inverse can be calculated numerically. Inverting the Dirac 
operator is the key computational step in LQCD and consumes over 90% 
of the CPU cycles. To construct the correlation functions requires Sp only 
from a given source point to all other sites of the lattice. This corresponds 
to 12 columns (or rows) of M~^, one each for the 12 spin and color degrees 
of freedom. 

• Correlation functions: The correlation hmctions are expressed as path 
ordered products of quark propagators and gauge fields using Wick contrac- 
tion of the operator fields. There are two kinds of correlation functions that 
are calculated gauge invariant correlation functions as discussed above 
for the case of the pion, or those in a fixed gauge. 

• Expectation values: The expectation value is simply an average of the 

correlation function evaluated over the set of "importance sampled" con- 
figurations. Even on a finite lattice the set of background gauge config- 
urations is infinite as the variables are continuous valued. Thus, while it 
is possible to calculate the correlation functions on specified background 
gauge configurations, doing the functional integral exactly is not feasible. 
It is, therefore, done numerically using monte carlo methods. 



20 



Rajan Gupta 



My recommendation for the simplest way to understand the computa- 
tional aspects of LQCD calculations is to first gain familiarity with the 
numerical treatment of any simple statistical mechanics system, for exam- 
ple the Ising model. The differences arc: (i) the degrees of freedom in 
LQCD are much more complicated - SU(3) link matrices rather than Ising 
spins, and quark propagators on these background configurations evaluated 
by inverting the Dirac operator; (ii) The action involves the highly nonlocal 
term Ln Det M which makes the update of the gauge configurations very 
expensive; and (iii) the correlation functions are not simple products of spin 
variables like the specific heat or magnetic susceptibility, but complicated 
functions of the link variables and quark propagators. 

The subtleties arising due to the fact that LQCD is a renormalizable 
quantum field theory and not a classical statistical mechanics system come 
into play when analyzing the behavior of the correlation functions as a 
function of the lattice spacing a, and in the quantum corrections that renor- 
malize the input parameters (quark and ghion fields and their masses) and 
the composite operators used in the correlation functions. At first glance 
it might seem that we have introduced an additional parameter, the lattice 
spacing a, in LQCD, however, as will be shown in Section 14, the coupling 
g and the cutoff a are not independent quantities but are related by the 
renormalization group as 

Aqcd = lim ie-i/2/5o«=(») {Po9'ia))-^^^"'° , (4.11) 
a— >o a 

where Aqcd is the non-perturbative scale of QCD, and /3o and /3i are the 
first two, scheme independent, coefficients of the /3-function. In statistical 
mechanics systems, the lattice spacing a is a physical quantity - the inter- 
molecular separation. In a quantum field theory, it is simply the ultraviolet 
regulator that must eventually be taken to zero keeping physical quantities, 
like the renormalized couplings {Aqcd and quark masses), spectrum, etc, 
fixed. 

The reason that lattice results are not exact is because in numerical 
simulations we have to make a number of approximations as discussed in 
Section 16. The size of the associated uncertainties is dictated by the 
computer power at hand. They are being improved steadily with computer 
technology, better numerical algorithms, and better theoretical understand- 
ing. I shall return to the issue of errors after discussing the formulation of 
LQCD. 
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5. Connection between Minkowski and Euclidean Field Theory 
and with Statistical Mechanics 



This section is a brief recapitulation of technical issues regarding the Eu- 
clidean path integral and extracting physics from Euclidean correlation 
functions. For details see [3,4] or the texts on statistical field theory, for 
example [22,23]. 

A £)-dimensional Minkowski field theory {D — 1 spatial dimensions and 
one time dimension) is connected to a D-dimensional Euclidean field theory 
through analytical continuation. Under Wick rotation 

xo=t ^ —ixi = —ir , 

Po = E ^ip4. (5.1) 
The Euclidean convention is 



2 \ ^ 2 2 _ J.2 _ 



Pl = i2p^=p'-E' = -Pl,. (5.2) 

i=l 

The connection between statistical mechanics and a Euclidean field the- 
ory in £)-dimensions is most transparent when the field theory is quantized 
using the Feynman path integral approach. This connection forms the ba- 
sis of numerical simulations using methods common to those used to study 
statistical mechanics systems. A synopsis of the equivalences between the 
two is given in Table 1. 

Let me start with a quick summary of the Euclidean the- 
ory. Consider the generic 2-point correlation function rM{t,k) = 
(0|/ (fxe-'''''T[Tix,t)J'\0,0)]\0) in Minkowski space. The Heisenberg 
operator T{x,t) can be written as 

jr(x, t) = e^nt-ipxjr^-mt+ipx (53^ 
so that 



TM{t,k) = {0\ J d3xe-*'= "(^e'^*-'f^jre-'«*+'f^jrt^|o) 



(5.4) 



where, in dropping the time-ordered product, it is assumed that t > 0. 
Using the invariance of the vacuum under time and space translations and 
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Table 1 

The equivalences between a Euclidean field theory and Classical Statistical Mechanics. 



Euclidean Field Theory 


Classical Statistical Mechanics 


Action 

unit of action h 

Fcynman weight for amplitudes 
Vacuum to vacuum amplitude 
Vacuum energy 

Vacuum expectation value ^o|o|o^ 

Time ordered products 

Green's functions (o T[C'i . . . On] o) 

Mass M 

Mass-gap 

Mass-less excitations 
Regularization: cutoff A 
Renormalization: A — + oo 
Changes in the vacuum 


Hamiltonian 

units of energy /3 = 1/feT 
Boltzmann factor e~^^ 

Partition function . e^l^^ 

^ — ' con J . 

Free Energy 

Canonical ensemble average 

Ordinary products 

Correlation functions (^Oi . . . On) 

correlation length ^ = 1/M 

exponential decrease of correlation functions 

spin waves 

lattice spacing a 

continuum limit a —* 

phase transitions 



assuming = T for brevity, this becomes 

= (0|(5(p - k) (|jre-'«*jr^ |0) . (5.5) 

Since the time dependence is exphcit, we can Wick rotate to Euclidean 
space using t ^ —it and keeping H unchanged 

TE(t) = {0\6{p - k) (^Te-'^^T^ |0) . (5.6) 

To get the path integral representation we divide the time interval into N 
steps T = Na and insert a complete set of states |<?!')(0| at each step. 

TE{t)=d{p-k) (0|j-|</.i)(<^i|e-'^'^|</)2>... 

(j>l...<j>N 

...(0;v-i|e-«"|(/)iv)(<^iv|:^|O) (5.7) 

The discussion of how, by this introduction of the complete set of states, 

one trades the operator form of the quantum mechanical amplitude by 
the sum over paths is standard and thus not repeated. On this "lattice" 
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the transfer matrix T = exp(— TYa) plays the role of the time translation 
operator, from which the continuum Hamiltonian is obtained as 

n=\im -- logT. (5.8) 
o^o a 

In order for 7i to be a self-adjoint Hamiltonian it is necessary that T is 
a symmetric, bounded, and positive operator acting on a Hilbert space 
of states with a positive norm. The necessary condition for this, is the 
Osterwalder-Schrader reflection positivity [24]. That this condition holds 
for LQCD has been shown by an explicit construction by M. Liischer [25]. 
Therefore, the correlation functions in the desired physical Minkowski the- 
ory can be obtained by analytical continuation of their Euclidean counter- 
parts. 

Assuming that the conditions for a sensible Euclidean theory have been 

established, the vacuum to vacuum amplitude (also called the generating 
functional or the partition function), in presence of a source J, is 

Z[J] = j #e-'5«+-^"^ (5.9) 

where 5*^; is the Euclidean action. Correlation functions can be obtained 
from this in the standard way by differentiating logZ[J] with respect to J, 
for example (o|(A(a;)(A(0)|o) = logZ[J]|j=o. 

Lastly, I summarize the relations between Euclidean and Minkowski 
quantities. The basic tool is the repeated insertion of a complete set of 
momentum states in Eq. 5.7. I use the relativistic normalization of states 

{k\p)=2ELHk^j,JkypJk^p. "--^ 2E{2'Kf5\p-k). (5.10) 

Then, for a spectrum of stable isolated states \n) with masses M„ that 
couple to an interpolating field T, 

n 

Now it is easy to obtain the connection between these M„ and the physical 
spectrum. A fourier transform in Euclidean time followed by a rotation to 
Minkowski space gives 

_e-M„|r| 1 1 



/ 



2M„ " 2M„(M„ - ipi) 2M„(M„ + ipi) ' 
1 

mI+pV 
^1 
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The Mn are precisely the location of the poles in the propagator of |n) . This 
is why one can simply short step the necessary analytical continuation as 
the masses measured from the exponential fall-off are the pole masses of 
the physical states. For r ^ oo only the lowest state contributes due to the 
exponential suppression, and this provides a way of isolating the lightest 
mass state. 

This simple connection holds also for matrix elements. Consider, for 
T2 >> Ti >> 0, the three-point function 



d'xd've-"''e-''<' (0|T[;F,(x,T5)C>(s,,T,):f,(0,0)]|0) = 

a b 

The two factors on the left and right of (ajOlfe) are each, up to one decay 
amplitude, the same as the 2-point function discussed in Eq. 5.11. They, 
therefore, provide, for large T2 — n and n, the isolation of the lowest states 
in the sum over states. These can be gotten from a calculation of the 2-point 
function. Consequently, ^ajOl^) measured in Euclidean simulations is the 
desired matrix element between these states and no analytic continuation 
to Minkowski space is necessary. 

These arguments break down once the sum over states is not just over 
stable single particle states. Consider, for example, an that is the p 
meson interpolating field. In that case there are also the 27r, An intermediate 
states with a cut beginning aX E = ^/2M^. The physical characteristics 
of p meson (mass and width) can be represented by a complex pole in the 
propagator that, in Minkowski space, lies on the second sheet under the cut. 
This physical pole can no longer be reproduced by measuring the "energy" 
given by the exponential fall-off of Euclidean correlation function at large 
T [26]. It requires measuring the correlation function for all intermediate 
states and momenta (all t), and then doing the fourier transform in r 
followed by the analytical continuation p4 —iE. 

For the same reason there does not exist a simple way for calculating 
n-particle scattering amplitudes. A general 2-particle scattering amplitude 
(4-point correlation function) is complex precisely because the intermedi- 
ate states involve a sum over resonances and a continuum with all values 
of momenta. In principle this can be obtained from Euclidean correlation 
functions but only by measuring them with essentially infinite precision for 
all possible r and then doing the analytical continuation. In practice this 
is a hopeless possibility. The one exception is the amplitude at threshold 
which is real. M. Liischer has shown how the 2-particle scattering am- 
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plitudcs at threshold can be calculated from the energy shift in a finite 
volume, the difference between the 2-particle energy and the sum of the 
two particles in isolation [27-30]. 



6. Formulation of Lattice Gauge Theories 



In 1974 Wilson [31] formulated Euclidean gauge theories on the lattice as 

a means to study confinement and carry out non-pcrturbativc analysis of 
QCD. The numerical implementation of the path integral approach requires 
the following five steps. 

— Discretization of space-time. 

— The transcription of the gauge and fermion degrees of freedom. 

— Construction of the action. 

— Definition of the measure of integration in the path integral. 

— The transcription of the operators used to probe the physics. 

Of these, the construction of the action and the operators is the most 
intricate. Rather than give you the most recent and improved version of 
the lattice action, which will be the subject of Martin Liischer's lectures, I 
have chosen to follow the original path of discovery. 



6.1. Discrete space-time grid 

There are a number of possible ways to discretize space-time in 4 Euclidean 
dimensions. These include hypercubic, body-centered cubic [32], and ran- 
dom [33] lattices. Of these the simplest is the isotropic hypercubic grid 
with spacing a = as = ax and size Ns x Ns x Ns x Nt- Since this grid 
is used almost exclusively, it will be the only one discussed in these sets 
of lectures. In the case of the random lattice even the free field analysis 
becomes difficult [33]. In fact very little exploration of variant grid ge- 
ometries has actually been done, as there is little expectation of significant 
physics gains to offset the additional computational complexity. Recently, 
proposals involving anisotropic lattices - a hypercubic grid with different 
spacing in space and time directions - have shown promise for reducing the 
discretization errors. I shall come back to this development in Section 18.3 
when discussing the glueball spectrum. 
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6.2. Lattice transcription of field variables i}{x) and An{x): 

The lattice transcription of the matter fields il>{x) is straightforward. The 
quark field is represented by anticommuting Grassmann variables defined 
at each site of the lattice. They belong to the fundamental representation 
of SU(3). The rules for integration over a pair of fermion variables tp and 
tjj are those standard for Grassmann variables 



Since the fermion action is linear in both ip and ip, these rules can be used 
to integrate over them, and the path integral reduces to one over only the 
gauge degrees of freedom. Thus, it turns out that in practice one docs not 
have to worry about transcribing grassmann variables on the lattice and 
implementing the Pauli exclusion principle. 

As an aside I would like to mention that the rules of integration for 
Grassmann variables are sufficient to do the fermion path integral. Con- 
sider {tp{xi) . . . 'tp{xn)tp{xi) . . . 'ijj{x„)). To get a non-zero expectation value 
one needs to saturate each site with a ip^j. The sites not covered by 
the fields in the operator have to come from the term in the ex- 

pansion exp5/ = l + 5/+<Sj/2 + ..., where V is the total number of 
sites. Since <S/ = ■tp{x){]p + m)tp{x) is a sum over sites, the number of 
terms in Sj~" poses an exponentially hard combinatorial problem. Fur- 
thermore, each term contributing to the connected diagram will contain a 
path-ordered product of gauge links (coming from the discretized Ip term 
in Sf. For example, see Eq. 7.5) connecting the sites in the operator, while 
the fermionic integration is ±m^~" with the sign given by the number of 
times the ip and tp are anti-commuted. You should convince yourself that 
the fermionic functional integral is the sum over all possible paths touching 
each site a maximum of one time. Having enumerated all such paths and 
their weights, the gauge functional integral now requires integrating over 
each link in each path. This is clearly a very formidable approach, and in 
practice rarely considered. 

The construction of gauge fields is less intuitive. In the continuum, the 
gauge fields A^{x) carry 4- vector Lorentz indices, and mediate interactions 
between fermions. To transcribe them Wilson noted that in the continuum 






(6.1) 
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a fcrmion moving from site a; to y in presence of a gauge field picks 
up a phase factor given by the path ordered product 

^(y) = pJ^^sA,{.)d., ^^^^ ^g_2) 

The eight types of gluons that mediate interactions between quarks are 
written in terms of the matrix A^{x) = A'^^{x) ■ Aa/2, where the group 
generators Aa are normahzed to TrXaXb = "^Sab- Eq. 6.2 suggested that 
gauge fields be associated with links that connect sites on the lattice. So, 
with each link Wilson associated a discrete version of the path ordered 
product 

U{x,x + fi) = Ui,{x) = e^s^'-^^+t), (6.3) 

where, for concreteness, the average field is defined at the midpoint of 
the link, and [/ is a 3 x 3 unitary matrix with unit determinant. Also, the 
path ordering in eqn. (2.3) specifies that 

U{x,x-fi) = U_^{x) = e-'^'sAA^-^) = U\x-ij.,x) . (6.4) 

6.3. Discrete Symmetries of the lattice theory 

The symmetry group of the continuum theory Poincarc invariance - is 
reduced to a discrete group. On a hypercubic lattice rotations by only 90° 
are allowed so the continuous rotation group is replaced by the discrete 
hypercubic group [34]. Translations have to be by at least one lattice unit, 
so the allowed momenta are discrete 

2-Kn 

k = — — n = 0,1, . . . L 
La 

or equivalently 

fc = ± ^ n = 0,l,...L/2. 

La 

On the lattice momentum is conserved m,odulo 2tt. 

In addition to the local gauge symmetry and Poincare invariance, the 
lattice action is invariant imder parity (V), charge conjugation (C) and 
time reversal (T). The behavior of the field variables under these discrete 
transformations is given in Table 2. 

6.4. Local Gauge Symmetry 

The effect of a local gauge transformation V{x) on the variables tlj{x) and 
U is defined to be 

tp{x) V{x)t{j{x) 
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Table 2 

The behavior of the gauge and fermion degrees of freedom under 
the discrete transformations of parity, charge-conjugation, and time- 
reversal. The charge-conjugation matrix C satisfies the relation 
C7pC~^ = — = —7^) and shall be represented by C = 7472- 
Note that in the Euclidean formulation T is a linear operator. 





V 


C 


r 


(74(33, r) 
Ui{x, t) 


U4{-x, r) 
U-i{-x, t) 


UI{x,t) 

U*{x,t) 


U-4{x, -r) 
Ui{x, -t) 


il>{x, t) 
tpix,r) 


74i/'(— a;, t) 
V'(-a;,T)74 


C'lp'^{x,T) 

-i/'^(a;,r)C-l 


7475'i/'(a:i -t) 
■)p{x, -r)7574 



U^{x) ^ V{x)U^{x)V\x + il) (6.5) 

where V{x) is in the same representation as the U , i.e., it is an 5C/(3) ma- 
trix. With these definitions there are two types of gauge invariant objects 
that one can construct on the lattice. 

— A string consisting of a path-ordered product of links capped by a fermion 
and an antifermion as shown in Fig. 5a. A simple example is 

Tr ^(.x) U^{x) U,Ax + fj.) . . . Up{y - p) i>{y) (6.6) 
where the trace is over the color indices. I will use the word string for a 
generalized version; a single spin and color trace as defined in Eq. 6.6, or 
a product of elementary strings that are path-ordered and do not have a 
spatial discontinuity. On lattices with periodic boundary conditions one 
does not need the -0, "0 caps if the string stretches across the lattice and is 
closed by the periodicity. Such strings are called Wilson/Polyakov lines. 

— Closed Wilson loops as shown in Fig. 5b. The simplest example is the 
plaquette, a 1 x 1 loop, 

Wl^^ = Re Tr {U^{x) UAx + ft) Ul(x + u) Ul(x)) . (6.7) 
Unless otherwise specified, whenever I use the term Wilson loops I will as- 
sume that the real part of the trace. Re Tr, has been taken. For SU {N > 3) 
the trace of any Wilson loop in the fundamental representation is com- 
plex, with the two possible path-orderings giving complex conjugate values. 
Thus, taking the trace insures gauge invariance, and taking the real part 
is equivalent to averaging the loop and its charge conjugate. 

Exercise: Convince yourself that there do not exist any other independent 
gauge invariant quantities. 
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Fig. 5. The two gauge invariant quantities, a) An ordered string of XJ' s capped by a 
fermion and an anti-fermion and b) closed Wilson loops. 



Exercise: Convince yourself that to preserve gauge invariance one has to 
formulate the lattice theory in terms of link variables J/^ rather than the 

directly. 

A gauge invariant action, therefore, has to be built out of foops and 
strings. However, note that these two basic gauge invariant objects can be 
of arbitrary size and shape. Furthermore, they can be taken to he in any 
representation of color SU (3). The only necessary constraint is that it give 
the familiar continuum theory in the limit a ^ 0. The reasons for wanting 
to preserve gauge invariance at all a is that otherwise one would have 
many more parameters to tune (the zero gluon mass, and the equality of 
the quark-gluon, 3-gluon, and 4-gluon couplings are consequences of gauge 
symmetry) and there would arise many more operators at any given order in 
a (these are irrelevant operators, but at finite a give rise to discretization 
errors as discussed in Section 8). In fact, one would like to preserve as 
many symmetries of the theory as possible at all values of a and not just 
in the continuum limit. It turns out to be relatively simple to maintain 
local gauge invariance, however, maintaining chiral symmetry, which has 
proven so important in the phenomenology of light quarks, and a one-to-one 
relation between Dirac and lattice fermion degrees of freedom has eluded 
us. This is discussed in Section 9. Also, discretization de facto breaks 
continuous rotational, Lorentz, and translational symmetry. The key idea 
of improved discretization schemes, to be discussed in Sections 8, 12, and 
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13, will be to implement the breaking of these symmetries by operators 
that come with sufficiently many powers of a so that their effect is very 
local, and thus small on physical (long-distance) quantities. 



7. Simplest Formulation of the Lattice Action 

7.1. Gauge Action 

The gauge action can be expressed in terms of closed loops. I will outline 
the construction for the simpler case of the abelian U{1) model in which the 
link variables arc commuting complex numbers instead of SU(3) matrices. 
Consider the simplest Wilson loop, the 1x1 plaquette W^^^: 

Expanding about x + gives 



tag[A^{x + ^)+A^(x+fi,+ ^)-A^{x+0+^)-A^{x+^)] -^^ 



e^p[ia''g{d^A,-d.A^) + '±l{dlA, - df^A^) + . . .] 

= 1 + iaVM- - ^F^.F'^- + 0(a6) + ... (7.2) 
The real and imaginary parts of the loop give 

(j4 2 

ReTr(l - W^^^) = —^F^^Fi^" + terms higher order in a 

lm{Wli^)=a''gF^, . (7.3) 

So far the indices /i and v are uncontracted. There arc 6 distinct positively 
oriented plaquettes, {ji < v}, associated with each site. Summing over 
/z, I/, and taking care of the double counting by an extra factor of |, gives 

l5]^ReTr(l-W^,ir) = ^^^F,^F^-^ \ j d^xF.^F^^ . 

X IX<V X ^,1^ 

Thus, to lowest order in a, the expansion of the plaquette gives the 
continuum action. The steps in the derivation for non-abelian groups are 
identical and the final expression (up to numeric factors) is the same. The 
result for the gauge action for SU(3), defined in terms of plaquettes and 
called the Wilson action, is 
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For historic reasons the lattice calculations arc mostly presented in terms 
of the coupling /3 = 6/(?^. Since /? is used for many other quantities like 
1 /kT or the /3-function, I shall try to be more explicit in the use of notation 
{g versus /?) if there is a possibility for confusion. 

Exercise: Repeat the Taylor expansion for the physical case of SU(3). Show 
that the non-abelian term in F^^ = d^A^ — d^A^ + A,y] arises due to 
the non-commuting property of A matrices. The details of the derivation 
for SU(2) are given in section 3 of Ref. [8]. 

There are four important points to note based on the above construction 
of the lattice action. 

1) The leading correction is O(a^). The term ^F^i,{d^^Au — df,Af^) is 
present in the expansion of all planar Wilson loops. Thus at the clas- 
sical level it can be gotten rid of by choosing an action that is a linear 
combination of say 1x1 and 1x2 Wilson loops with the appropriate 
relative strength given by the Taylor expansion (see Section 12). 

2) Quantum effects will give rise to corrections, i.e. o? X{g'^)a^ where in 
perturbation theory X{g'^) = 1 + Cig"^ + . . ., and will bring in additional 
non-planar loops. Improvement of the action will consequently require 
including these additional loops, and adjusting the relative strengths 
which become functions of g^. This is also discussed in Section 12. 

3) The reason for defining the action in terms of small loops is compu- 
tational speed and reducing the size of the discretization errors. For 
example the leading correction to 1 x 1 loops is proportional to 
whereas for 1 x 2 loops it increases to 5a^/12. Also, the cost of simula- 
tion increases by a factor of 2 — 3. 

4) The electric and magnetic fields E and B are proportional to F^^. 
Eq. 7.3 shows that these are given in terms of the imaginary part of 
Wilson loops. 

7.2. Fermion Action 

To discretize the Dirac action, Wilson replaced the derivative with the 
symmetrized difference and included appropriate gauge links to maintain 
gauge invariance 

^^V= ^i^{x) iM^M^ + f^) - Ul{x-jl)'4,{x-ij)\ . (7.5) 
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It is easy to sec that one recovers the Dirac action in the hmit a ^ by 
Taylor expanding the and + fi) in powers of the lattice spacing a. 
Keeping only the leading term in a, Eq. 7.5 becomes 

^V^(x)7^[ (1 + iagAi,{x + |) + . . . ) {.jjix) + a^'{x) + ...)- 

(1 - iagA^ix - f )) + ••• ) (i'ix) - a^j:'{x) + ...)] 

= ^{x)^^{d^+'^dl + ...)i;{x) 

+igii{x)^M^^ + '^{\dlA^ + [dMd^^ + A^dl) + ■ ■ -W^) , 

which, to O(a^), is the kinetic part of the standard continuum Dirac action 
in Euclidean space-time. Thus one arrives at the simplest (called "naive" ) 
lattice action for fermions 

= mg y^^ijj{x)ip{x) 

X 

+ ^ ^i'{xhi^[Ufj,{x)il>{x + fi) - Ulix - fi)ilj{x - fi)] 

X 

= Yj>{x)M^ym{y) (7.6) 

X 

where the interaction matrix is 

M^.[C/] = m,5i, + - (7.7) 

The Euclidean 7 matrices are hermitian, 7^ = 7!^, and satisfy {7^,71/} = 
25 The representation I shall use is 

. - {.I %) . = (i \) -. = (;;). (-) 

which is related to Bjorken and Drell conventions as follows: 7i = i'ysDi 
74 = 7b _D) 75 = 1%D- representation 71,73 are pure imaginary, 

while 72,74,75 are real. 

The Taylor expansion showed that the discretization errors start at 
O(a^). For another simple illustration consider the inverse of the free- field 
propagator m + i/ a^^^^sin{p^a). Set p = and rotate to Minkowski 
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space {p4 iE, i.e. sinp4a isinhEa). Then, using the forward propaga- 
tor (upper two components of 74), gives 

mP°'^a = sinhEa (7.9) 

for the relation between the pole mass and the energy. This shows that, 
even in the free field case, the continuum relation E{p = 0) = m is violated 

by corrections that are O(a^). 

Symmetries: The invariance group of the fermion action under rotations 
in space and time is the hypercubic group. Full Euclidean invariance will 
be recovered only in the continuum limit. The action is invariant under 
translations by a and under V, C, and T as can be checked using Table 2. 
The naive action tjj^M^yipy has the following global symmetry: 

i){x) lp{x)e-'^ (7.10) 

where ^ is a continuous parameter. This symmetry is related to baryon 
number conservation and leads to a conserved vector current. For rUg = 
the action is also invariant under 

^{x) e**''=i/'(x) 

^(x) Jj{x)e''''^'^ (7.11) 

Having both a vector and an axial symmetry in a hard-cutoff regularization 
would imply a violation of the Adler-Bell- Jackiw theorem. It turns out that 
while the naive fermion action preserves chiral symmetry it also has the 
notorious fermion doubling problem as explained in Section 9. The chiral 
charges of these extra fermions are such as to exactly cancel the anomaly. 
In fact the analysis of S'^ lead to a no-go theorem by Nielsen-Ninomiya [14] 
that states that it is not possible to define a local, translationally invariant, 
hermitian lattice action that preserves chiral symmetry and does not have 
doublers. 

We will discuss two fixes to the doubling/anomaly/chiral symmetry 
problem. First, include an additional term in the action that breaks chi- 
ral symmetry and removes the doublers (Wilson's fix). Second, exploit 
the fact that the naive fermion action has a much larger symmetry group, 
C/y (4) (g) J7a(4), to reduce the doubling problem from 2'' = 16 — > 16/4 while 
maintaining a remnant chiral symmetry (staggered fermions) . In both cases 
one regains the correct anomaly in the continuum limit. Before diving into 
these fixes, I would like to first discuss issues of the measure of integra- 
tion and gauge fixing, and give an overview of Symanzik's improvement 
program. 
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7.3. The Haar Measure 

The fourth ingredient needed to complete the formulation of LQCD as a 

quantum field theory via the path integral is to define the measure of inte- 
gration over the gauge degrees of freedom. Note that, unlike the continuum 
fields Aj^, lattice fields are SU(3) matrices with elements that are bounded 
in the range [0, 1]. Therefore, Wilson proposed an invariant group measure, 
the Haar measure, for this integration. This measure is defined such that 
for any elements V and W of the group 



where f(U) is an arbitrary function over the group. This construction is 
simple, and has the additional advantage that in non-perturbative studies 
it avoids the problem of having to include a gauge fixing term in the path 
integral. This is because the field variables arc compact. Hence there are 
no divergences and one can normalize the measure by defining 



For gauge invariant quantities each gauge copy contributes the same 
amount. Consequently the lack of a gauge fixing term just gives an overall 
extra numerical factor. This cancels in the calculation of correlation func- 
tions due to the normalization of the path integral. Note that gauge fixing 
and the Fadcev-Papov construction does have to be introduced when doing 
lattice perturbation theory [36]. 

The measure and the action are invariant under C, V, T transformations, 
therefore when calculating correlation functions via the path integral we 
must average over a configuration U and the ones obtained from it by the 
application of C, V, T. In most cases this averaging occurs very natu- 
rally, for example. Euclidean translation invariancc guarantees V and T 
invariance. Similarly, taking the real part of the trace of a Wilson loop is 
equivalent to summing over U and its the charge conjugate configuration. 
A cook-book recipe for evaluating the behavior of correlation functions af- 
ter average over configurations related by the lattice discrete symmetries 
is given in Section 17. 

One important consequence of these symmetries is that only the real 
part of Euclidean correlation functions has a non-zero expectation value. 
This places restriction on the quantities one can calculate using LQCD. For 
example, simulations of LQCD have not been very successful at calculations 
of scattering amplitudes, which are in general complex [48] . The connection 




(7.12) 




(7.13) 
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between Minkowski and Euclidean correlation functions was discussed in 
Section 5. To get scattering amplitudes what is required is to first calculate 
Euclidean amplitudes for all values of the Euclidean time r, followed by a 
fouricr transform in r, and finally a Wick rotation to Minkowski spacx;. For 
such an procedure to work requires data with phenomenal accuracy. Unfor- 
tunately, simulations give Euclidean correlation functions for a discrete set 
of T and with errors that get exponentially magnified in the rotation. The 
best that has been done is to use Liischer's method that relates scattering 
phase shifts to shifts in two particle energy states in finite volume [27,28]. 
So far this method has been tested on models in 2-4 dimensions, and used 
to calculate scattering lengths in QCD [29,30]. 



7.4- Gauge Fixing and Gribov Copies 

Gauge invariance and the property of group integration, / dU U = 0, are 
sufficient to show that only gauge invariant correlation functions have non- 
zero expectation values. This is the celebrated Elitzur's theorem [37]. It 
states that a local continuous symmetry cannot be spontaneously broken. 

In certain applications (study of quark and gluon propagators [38], 3- 
gluon coupling [39], hadronic wavcfunctions [40], and determination of 
renormalization constants using quark states [41]) it is useful to work in 
a fixed gauge. The common gauge choices are axial (r/^^/i = where 
77 in Minkowski space is a fixed time-like vector. The common choice is 
A4 = 0), Coulomb {diAi = 0), and Landau {d^An = 0). Their lattice 
analogues, assuming periodic boundary conditions, are the following. 

Maximal Axial Gauge: On a lattice with periodic boimdary condi- 
tions one cannot enforce ^4.4 = by setting all links in the time direction 
to unity as the Wilson line in each direction is an invariant. The gauge 
freedom at all but one site can be removed as follows. Set all time- like 
links ?74 = 1 on all but one timeslice. On that time-slice all U4 links are 
constrained to be the original Wilson line in 4-direction, however, one can 
set links in z direction 1/3 = 1 except on one z-plane. On that z-plane set 
U3 links to the Wilson line in that direction and all links in y direction 
1/2 = 1 except on one y-line. On this y-line set all U2 links to the value 
of y- lines and links in x direction Ui = 1 except on one point. While such 
a construction is unique, the choice of the time-slice, z-plane, y-line, and 
x-point and the labeling of x, y, z, t axis are not. Also, this gauge fixing 
does not respect Euclidean translation invariance and it is not smooth. 

Coulomb Gauge: The lattice condition is given by the gauge transfer- 
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mation V{x) that maximizes the function 

F[V] = ^^'Re TrV(x) (Uiix)V^x + i) + U}{x- i)V^{x - iU (7.14) 

separately on each time-shce and i runs over the spatial indices 1 — 3. 
Landau Gauge: This is defined by the same function as in Eq. 7.14 

except that the sum is evaluated on the whole lattice and i = 1 — 4. 

The maximization condition, Eq. 7.14, corresponds, in the continuum, 
to finding the stationary points of F[V] = \\A^\\ = J d'^xTT{A'^)'^{x). 
What wc would like is to find the global minimum corresponding to the 
"smoothest" fields A^, or the copies in case it is not unique. For any gauge 
configuration there is an infinite set of gauge equivalent configurations ob- 
tained by applying the set of transformations {^(.t)}. This set is called the 
gauge orbit. The first question is whether the gauge fixing procedure gives 
a unique solution, i.e. whether the gauge-fixing algorithm converges to the 
same final point irrespective of the starting point on the gauge orbit. For 
the maximal axial gauge, the construction is unique (the configuration {U} 
and {VUV^} for any {V} give the same final link matrices), so the answer 
is clearly yes. For the Coulomb and Landau gauge choices the answer is 
no - there are a number of solutions to the gauge-fixing condition. These 
are the famous Gribov copies [42]. For a recent theoretical review of this 
subject see [43]. 

The algorithms for fixing to either the Coulomb gauge or the Landau 
gauge are mostly local. Local algorithms generate, for each site, a gauge 
matrix Vi{x) that maximizes the sum, defined in Eq. 7.14, of the link ma- 
trices emerging from that site. This gauge transformation is then applied 
to the lattice. These two steps are iterated until the global maximum is 
found. The ordered product Yli 'Vi{x) gives V{x). For all practical purposes 
we shall call this end-point of the gauge fixing algorithm a Gribov copy. It 
is, in most cases and especially for the local gauge-fixing algorithms, the 
nearest extremum, and not the global minimum of A set of Gri- 

bov copies can be generated by first performing a different random gauge 
transformation on the configuration and then applying the same gauge- 
fixing algorithm, or by using different gauge-fixing algorithms. For recent 
studies of the efficacy of different algorithms see [44,45] . 

The Monte Carlo update procedure is such that it is highly improbable 
that any finite statistical sample contains gauge equivalent configurations. 
So one does not worry about gauge copies. In the gauge-fixed path integral 
one would a priori like to average over only one Gribov copy per config- 
uration and choose this to be the "smoothest". Finding the V{x) which 
gives the smoothest configuration has proven to be a non-trivial problem. 
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The open questions are: (i) how to define the "smoothest" solution of the 
gauge-fixing algorithm (i.e. defining the fundamental modular domain), (ii) 
how to ascertain that a solution is unique (lies within or on the boundary 
of the fundamental modular domain), and (iii) how to find it. 

The ultimate question is - what is the dependence of the quantities one 
wants to calculate on the choice of Gribov copy? A number of studies have 
been carried out [46], however the results and their interpretation are not 
conclusive. Firstly, it is not known whether the smoothest solution lying 
within the fundamental modular domain can be found economically for all 
configurations and over the range of gauge coupling being investigated. In 
the absence of the smoothest solution it is not clear whether (i) the residual 
gauge freedom simply contributes to noise in the correlation function, and 
(ii) one should average over a set of Gribov copies for each configuration to 
reduce this noise. In my opinion these remain challenging problems that 
some of you may wish to pursue. 

The problem of Gribov copies and the Faddcev-Papov construction (a 
gauge fixing term and the Faddeev-Papov determinant) resurfaces when 
one wants to do lattice perturbation theory or when taking the Hamilto- 
nian limit of the path integral via the transfer matrix formalism. This 
construction is analogous to that in the continuum [36]. Finally, let me 
mention that in a fixed gauge, the lattice theory has only a BRST invari- 
ancc [27]. While BRST invariance is sufficient to prove that the theory 
is renormalizable, the set of operators of any given dimension allowed by 
BRST invariance is larger than that with full gauge invariance. In the 
analyses of improvement of the action and operators it is this enlarged set 
that has to be considered [47]. 

8. Continuum limit and Symanzik's Improvement program 

The Taylor expansion of the simple Wilson action for the gauge and fermion 
degrees of freedom showed the presence of higher dimensional operators. 
Since these are suppressed by powers of the lattice spacing, they are ir- 
relevant operators in the language of the renormalization group, i.e. they 
vanish at the fixed point at a — > 0. Another way of expressing the same 
thing is as follows. The lattice theory with a hard cut-off at ir/a can be re- 
garded as an effective theory [49]. Integration of momenta from oo n/a 
generates effective interactions. Thus, in general, we can write the action 
as 




(8.1) 



n 



I 
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where Sc{Sl) is the contimmm (lattice) action, and the first sum is over all 
operators of a given dimension, and the second over all dimensions greater 
than four. Similarly, the operators used to probe the physics can be written 
as 

Oc = OLi-Jr/a) + E E (^""DfOn^M (8.2) 

n i 

where in this case quantum corrections can induce mixing with lower di- 
mension operators also. The coefficients C" and -D" are both functions 
of the coupling a^. Parenthetically, one important reason for preserving 
gauge, Lorentz, and chiral symmetries at finite a is that they greatly restrict 
the set of possible operators of a given dimension. 

Now consider the expectation value of a given operator 

{Oc) = J dUOce-^^ 

n i 

= ( dU{OL + Y^Y^a''D'^0'l]{l-Y^^a'^CY'AT + ■..]e-^'^ 

n i mi 

= J dUOLB-^'^ 

+ j (i[7(0i{EE«"^r^r + •••}) e"^" 

+ / ^^t^(EE«"^"'^"{i-EE«'"^™'4r + ---})e-^"-(8-3) 

n i mi 

Note that the operators A^, which appear in the action, are summed over 

all space-time points. This leads to contact terms in Eq. 8.3 which I shall, 
for brevity, ignore as they necessitate a proper definition of the operators 
but do not change the conclusions. 

The basis for believing that LQCD provides reliable results from sim- 
ulations at finite a is that the contributions of the lattice artifacts (last 
two terms in the last expression) are small, can be estimated, and thus 
removed. At this point I will only make a few general statements about 
the calculation of expectation values which will be embellished on in later 
lectures. 

— The contribution of the last two terms in Eq. 8.3 vanishes in the limit a — > 
unless there is mixing with operators of lower dimension. In such cases 
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the mixing cocfBcicnt has to be determined very accurately (necessitating 
non-perturbative methods) otherwise the a ^ hmit is divergent. 

— Calculations have to be done at values of a at which contributions of 
irrelevant terms organized in successive powers in a decrease at least geo- 
metrically. 

— To achieve continuum results simulations are carried out at a sufficient 
number of values of a to allow reliable extrapolation to a = 0. The "fine- 
ness" of the lattice spacing a at which to do the calculations and the number 
of points needed depends on the size of the corrections and how well we 
can determine the fimctional form to use in the extrapolation. 

— To improve the result to any given order in a, one has to improve both 
the action and the operators to the same order. Such an approach, based 
on removing lattice artifacts organized as a power series expansion in a, is 
called the Symanzik improvement program. 

Eq. 8.3 shows clearly that there is fair degree of flexibility in constructing 
LQCD. We are free to add any irrelevant operator with a sensible strength 
and still recover QCD in the continuum limit. The only difference between 
various constructions of the action is the approach to the continuum limit. 
What we would like to do is improve the action and operators so as to 
get continuum results on as coarse lattices as possible. On the other hand 
improvement requires adding irrelevant terms to both the action and the 
operators. This increases the complexity of the calculations and the sim- 
ulation time. Thus, the strategy for minimizing discretization errors is a 
compromise - optimize between simplicity of the action and operators (de- 
fined in terms of the cost of simulation, code implementation, and analyses) 
and the reduction of the discretization errors evaluated at some given fixed 
value of the lattice spacing, say a w 0.1 fermi. The remaining errors, pre- 
sumably small, can then be removed by extrapolating the results obtained 
from simulations at a few values of a to a = 0. 

The number of floating points operations needed in the simulations of 
LQCD scale roughly as in the quenched approximation and L^*^ with dy- 
namical fermions. Current lattice simulations of QCD are done for lattice 
spacing a varying in the range 2 GeV < < 5 GeV. For these pa- 
rameter values the corrections due to the 0{a), 0{m,a),0{pa),0{a^), . . . 
terms in the Wilson action are found to be large in many observables. The 
importance of improving the lattice formulation (action and operators) is 
therefore self-evident, especially for full QCD. Thus, the search for such 
actions is a hot topic right now. Before discussing these let me first sum- 
marize what one hopes to gain by, and the criteria by which to judge, such 
an improvement program. 
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1. Improved scaling, i.e. discretization errors are reduced. Consequently 
one can work on coarser lattices (smaller L = Ns) for a given accuracy. 
Even a factor of 2 reduction in L translates into 2^° w 1000 in computer 
time! 

2. Better restoration of rotational and internal symmetries like chiral, 
staggered flavor, etc.. 

3. The trajectory along which the continuum limit is taken should not 

pass close to extraneous phase transitions that are artifacts of the lat- 
tice discretization (see Section 15 for a discussion of phase transitions 
in the lattice theory). The problem with such transitions is that in 
their vicinity the desired continuum scaling behavior can be modified 
for certain observables, i.e. the artifact terms in Eq. 8.3 may not have a 
geometric convergence in a and may thus be large. 

4. The study of topology on the lattice is particularly sensitive to short 
distance fluctuations that are artifacts of discretizing the gauge action. 
For the Wilson gauge action one finds that the topological charge, cal- 
culated using Lilscher's geometric method [50], jumps from ^ 1 at a 
core radius of the instanton of pcU « 0.7, whereas the action of such 
instantons, Sp^ , is significantly less than the classical value S**^' = Stt^//?. 
The entropy factor for such small instantons with S < 5^' overwhelms 
the Boltzmann suppression, leading to a divergent topological suscepti- 
bility in the continuum limit. These short distance artifacts are called 
dislocations and it is not a posteriori obvious how to separate the phys- 
ical from the unphysical when p « pc. It is, therefore, very important 
to start with a lattice formulation in which such artifacts are excluded, 
and there exist only physical instantons, i.e. those with mean core size 
that has the correct scaling behavior with g. One finds that dislocations 
are suppressed by improved actions, i.e. improving the gauge action in- 
creases Sp^ [51] and better definition of the topological charge [52]. 

5. Lattice perturbation theory becomes increasingly complicated as more 
and more irrelevant operators are added to the action. Even though, 
with the advent of non-perturbative methods for determining renor- 
malization constants, the reliance on lattice perturbation theory has 
decreased, perturbative analyses serve as a valuable guide and check, 
so we would like to retain this possibility. What one expects for an 
improved action is that the lattice and continuum value of a,,, and the 
coefficients of the series expansion in as for a given observable, are much 
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closer. In that case a much more reliable matching between the two the- 
ories is possible. At present the 1-loop matching between the lattice and 
continuum theories is one of the larger sources of uncertainties. 

6. The Adler-Bell-Jackiw (ABJ) Anomaly: In the continuum the flavor 
singlet axial current has an anomaly. This shows up as part of the 
UV regularization scheme. In the Pauli-Villars scheme one introduces 
a heavy mass as the cutoff. This breaks 75 invariance, while in dimen- 
sional regularization, 75 is not well-deflned away from 4 dimensions. 
Naive fermions, as shown below, preserve the 75 invariance but at the 
expense of flavor doubling. A consequence of the doublers is a pair- 
wise cancellation and no net anomaly. A viable lattice theory should 
reproduce the ABJ anomaly. 

Before discussing such "improved" actions, it is instructive to revisit the 
construction of the naive Dirac action to highlight its problems. 

9. Fermion Doubling Problem in the "naive" action 

The problem with the naive discretization of the Dirac action, specified 
by the operator in Eq. 7.7, is that in the continuum limit it gives 
rise to 2'^ = 16 flavors rather than one. This doubling problem is readily 
demonstrated by the inverse of the free field propagator (obtained by taking 
the fourier transform of the action with all Ufj,{x) = 1). 



which has 16 zeros within the Brillouin cell in the limit — !■ 0. Defining 
the momentum range of the Brillouin cell to be {— 7r/2, 37r/2}, the zeros 
lie at = and tt. As discussed later, this proliferation holds under very 
general conditions specified by the Nielsen- Ninomiya theorem [14], and 
is intimately related to simultaneously preserving chiral symmetry. The 
inclusion of the gauge fields does not solve the doubling problem. 

Before discussing partial fixes, it is useful to investigate the properties 
of these extra zero modes under chiral transformations. Let us define a 
set of 16 4- vectors = {(0, 0, 0, 0), (tt, 0, 0, 0), . . . , (tt, tt, tt, tt)} with A = 
{1...16}, and consider the expansion of the massless propagator about 
these points. Then 




(9.1) 
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Tabic 3 

The chiral sign and degeneracy factor for the doublcrs. 



A 


Degeneracy 




(0,0,0,0) 


1 


+1 


(1,0,0,0) 


4 


-1 


(1, 1,0,0) 


6 


+1 


(1,1,1,0) 


4 


-1 


(1,1,1,1) 


1 


+1 



= ^^7^ sin k^a 

= -^7^ sinfc^a , (9.2) 

where iS^ = {+1,-1} depending on whether the fj, component of 11^ is 
{0, n}. In the last expression a new representation of gamma matrices has 
been defined by the similarity transformation 'y^ = YA7^y| where 

Ya = n(7M75)< (9.3) 

and are {0, 1} depending on whether the momentum expansion in the 
direction is about or tt. Now, 75 ~ YaJsY^ — SfS2S^S^j5 = X^'y^, 
with the signs given in table 3. Thus, the sixteen species break up 
into two sets of 8 with chiral charge ±1, and render the theory anomaly- 
free. Consequently, this "naive" discretization is phenomenologically not 
acceptable. The presence of doublers can, in fact, be traced to a larger 
flavor symmetry of the naive fermion action as I now discuss. 

In the continuum the mass-less Dirac action is invariant under the 
U{l)v ® U{1)a flavor transformations 

i^R ^ Vni^R, ^ V^^y^ . (9.4) 

Since the left and right handed fields arc only coupled through the mass 
term, the vector symmetry Vl = Vr = e*^ holds for all m, whereas the 
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axial symmetry Vl = = e*^ (or cquivalcntly Vl = Vr = e^^^-') holds 
only for m = 0. The naive fermion action iSjv for a single flavor has a much 
larger symmetry group, U{4)v (S) U{4:)a, under 

V^(x) ^ V^(x)r, (vl^{l + e(a:)75) + V^^{l - 6(a;)75))rl (9.5) 

if one considers translations by 2a. is given in Eq. 11.1, the phase e{x) 
in Eq. 11.4, and V^^^^^ = t{x)^^. Here Vh,R are general 4x4 unitary 
matrices acting on the spinor index of the fermions. Even though there is, 
in Eq. 9.5, a dependence on x through t{x), this is a global symmetry as 
specifying it at one point fixes it at all points. 

To see how doubling arises, consider the vector (Vl = Vr) discrete 
subgroup with the thirty two elements {±^4} defined in Eq. 9.3. Eq. 9.2 
shows that the Ya shift the momenta hy -K/a along the direction specified 
by A 

YASF{k)Yl = Spik + ^A) . (9.6) 

Since this is a similarity transform, the sixteen regions of the Brillouin 
zone related by this transormation are physically equivalent. This analysis, 
which relied only on the spinor structure of the action, makes it easy to 
see why the doublers exist even for the interacting theory. In Section 11 
this symmetry will be exploited to construct staggered fermions. In this 
formulation the sixteen doublers are reduced to four, and at the same time 
a U{l)v (8) U{l)e symmetry (corresponding to when Vl,r are just phases) 
is retained. The U^l),, plays the role of a chiral symmetry at finite a and 
in the continuum limit becomes the U{1)a- 

9.1. Generalities of Chiral Symmetry 

In the massless limit of the continuum theory, the flavor symmetry at the 
classical level is Uv{^) x Ua{^) x SUv{nf) x SUA{nf), where n/ is the 
number of light fermions. The unbroken Uv{i) gives rise to baryon number 
conservation. The Ua{^) is broken by instanton contributions and the flavor 
singlet axial current is anomalous. The SUA{nf) is spontaneously broken 
and the quark condensates acquire a non-zero expectation value in the 
QCD vacuum. Associated with this spontaneously broken symmetry are 
(riy — 1) goldstone bosons, i.e. the pions. Chiral symmetry has played a 
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very useful role in the continuum theory. For example, one can (i) classify 
operators under distinct representations of SUL{nf) x SUairif), and chiral 
symmetry prohibits mixing between operators in different representations, 
(ii) derive relations between matrix elements of a given operator between 
states with different number of pions in the final state (for example K 
TTTT -i-^ K ^ tttttt), (hi) exploit the operator Ward identity, and (iv) derive 
relations between correlation functions, i.e. the Ward-Takahashi identities. 
It is therefore desirable to preserve chiral symmetry at finite a and not 
just in the continuum limit. To determine if this is possible one needs to 
understand the relationship between doublers and chiral symmetry in the 
lattice theory. 

The chiral symmetry defined in Eq. 7.11 is realized by the lattice theory 
provided 75 M + M75 = 0. It is easy to verify that the naive lattice action 
satisfies the hermiticity property 

75M75 = Mt (9.7) 

Thus, in the massless limit, the Euclidean lattice action has to be anti- 
hermitiau for chiral symmetry, Eq. 7.11, to hold. This is true of 
given in Eq. 7.7, and was accomplished by simply taking the symmetric 
difference for the derivative. The important question is - what are the 
general conditions under which a lattice theory with this realization of 
chiral symmetry exhibits doublers. The answer is given by the Nielsen- 
Ninomiya theorem [14] . Consider a generalized action such that S~^{p, m = 
0) = iFip). Then, if 

• the function F{p) is periodic in momentum space with period 27r/a. 
(This is a consequence of translation invariance.) 

• the lattice momenta are continuous in the range {0, 27r}. This is true in 
the i — > 00 limit at all values of the lattice spacing. 

• F{p) is continuous in momentum space. This is guaranteed if the inter- 
actions are local. 

• F{p) Pn'Jij. for small p^ and as a —> to match the continuum theory. 

• the action has chiral symmetry, i.e. it satisfies the hermiticity property 
75M75 = and the massless Dirac action is anti-hermitian. 

the theory will have doublers. It is easy to understand this from the 

schematic one dimensional plot of F{p) shown in Fig. 6 which highlights 
the fact that F{p) has to have an even number of zeros under the above 
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-3t/2 




37t/2 



Fig. 6. A schematic 1-d representation of the massless inverse propagator F(j)). 

conditions. The ± slopes imply that one state is left moving while the other 
is right moving. In l-dimcnsion the direction is the same as handedness, 
consequently left-handed and right-handed fields come in inseparable pairs, 
and together form a Dirac fermion. 

The key consequence of the Nielscn-Ninomiya theorem is that one cannot 
define a lattice theory with chiral fermions. For example if one writes 
down for a left-handed Dirac fermion, then its doubler will be right- 
handed but have the same transformation under L $iR. Consequently, the 
two "states" coexist to form a Dirac fermion. Assuming that the Nielsen- 
Ninomiya theorem is air-tight, then how does one proceed to (i) define a 
lattice theory without doublers and which recovers chiral symmetry in the 
continuum limit, and (ii) define chiral gauge theories on the lattice? 

The answer to the first question will be given by explicit reconstruction of 
Wilson's and staggered fixes. An alternate very exciting new development, 
which I will not do justice to, is the resurrection [53] of the Ginsparg- Wilson 
condition for a chirally symmetric formulation [54] 



Liischer has recently presented a clear analysis of how this condition al- 
lows the construction of a lattice theory without doublers and with chiral 

symmetry, avoids the NN theorem, and gives the ABJ anomaly [55]. Un- 
fortunately, the bottleneck with this alternative is the lack of a closed form 



aD-f^D . 



(9.8) 
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expression for the lattice Dirac action that is sufficiently local to make 
simulations feasible. 

The second question - construction of chiral gauge theories - is still 

unresolved. For those interested in this very important unsolved problem, 
a good starting point are the reviews [15-17] which contain some of the 
recent efforts. 

10. Wilson Fermions (WF) 

In Section 8, it was shown that the lattice action is not unique. One has the 
freedom to add an arbitrary number of irrelevant operators to the action, as 
these do not change the continuum limit. Wilson's solution to the doubling 
problem was to add a dimension five operator ar^Oip, whereby the extra 
fifteen species at = w get a mass proportional to r/a [56]. The Wilson 
action is 

= TUg ^ii{x)llj{x) 

X 

+ ^ ^'4}{x)li^[Ui^{x)'il){x + ji) - Ul{x - fi)^{x - fi)] 

X,fJ, 

- ^ ^?(a;)[?7;.(x)V'(a; + A) -2^p{x)+ Ul{x - jl)i){x - jx)] 

X,fJ. 

{rUga +4r) 

= — ^ 2^i;{x)^p{x) 

X 

+ ^ '^i'iiln - r)Ui^{x)il>{x + fi)- (7a( + r) U^ix - fi)il>{x - fij] 

X 

= j^^tMZi^;; (10.1) 

where the interaction matrix is usually written as 

with the rescaling 

K= l/(2m,a + 8r) 
= ^ruga + 4r V = i^l^J^k . (10.3) 
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In this form the Wilson action has the foUowing interpretation as a Sta- 
tistical Mechanics model. The local term acts like an inertia keeping the 
fermion at the same site while the non-local term makes the fermion hop 
to the nearest neighbor site with strength k. During this hop the fermion 
acquires a twist by (7^ — r) in spin space and by 1/^ in color space. The 
case r = 1 is special as the spin twist 1 ± 7^ is a projection operator of 
rank 2, 

T>(^2») = 2 . (10.4) 

Later we shall sec that for r = 1 the dispersion relation has a single branch, 
i.e. there are no doublers in the free field limit. 

From Eq. 10.3, one notes that for the free theory the quark mass is given 
in terms of the lattice parameters k and r as 

rUga = ^ -Ar= ^ - (10.5) 

with a zero at k = Kc = l/8r. For the interacting theory (C/^(x) ^ 1) 
we will continue to define m^a = 1/2k — 1/2kc with the proviso that Kc 
depends on a. The renormalization of Kc independent of n implies that the 
quark mass has both multiplicative and additive renormalizations. This is 
a consequence of the explicit breaking of chiral symmetry by the irrelevant 
term proportional to r. Wilson's fix for doublers comes with a high price 
tag - a hard breaking of chiral symmetry at 0{a). 

The free propagator in momentum space for Wilson fermions is 

Sp{p) = M^\p) = --^ ; : . . (10.6) 

1 - 2k 2^^[r cosp^a-«7p smp^aj 

Due to the dimension 5 operator (terms proportional to r) the 15 extra 
states at p^a = TT get masses of order 2rn/a, where n is the number of 
components oi p^a = tt. 

The operator M satisfies the following relations: 

75Mly75 = Mw, 

lhSp{x,y)-iB = SF{y,x) , 

Ml{K,r) = Mwi-K,-r) . (10.7) 

The first two state that the "hermiticity" property of is preserved. The 
second equation relates a quark propagator from x ^ y to the antiquark 
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propagator from y ^ x. This important identity, called liermitieity or 75 
invariance, leads to a significant simplification in numerical calculations as 
discussed in Section 17.1. The adjoint in Sp is with respect to the spin 
and color indices at each site. The last relation shows that is not 

anti-hermitian due to the Wilson r term. 

An analysis similar to Eq. 7.9 for the naive fermion action, shows that 
the pole mass derived from Eq. 10.6 is different from the bare mass and 
given by 

mP°'«a = r(cosh£;a - 1) + sintiEa . (10.8) 

This shows, as expected, that the discretization corrections in spectral 
quantities occur at 0{a). 



10.1. Properties of Wilson Fermions 



I give a brief summary of the important properties of Wilson fermions. 

• The doublers are given a heavy mass, 2rn/a, and decouple in the con- 
tinuum limit. 

• Chiral symmetry is broken explicitly. The derivation of axial Ward iden- 
tities, using the invariance under the transformation given in Eq. 7.11, 
have the generic form 

For WF, the variation of the action under an axial transformation is 
d5S 



de 



d^^Af, - 2mP + raX, (10.10) 



where X is an additional term coming from the variation of the Wilson 
r term [57]. Thus, in general, all relations based on axial WI will have 
corrections at order ra and involve mixing with operators that would 
normally be absent due to the chiral symmetry. 

• Calculation of the matrix element (0|T[A^V;y Vp] |G) shows that the ABJ 
anomaly is correctly reproduced in the continuum limit [58]. Karsten 
and Smit show, by an explicit 1-loop calculation, that while each of the 

extra 15 states contribute terms proportional to r, the total contribution 
from all sixteen states is independent of r and equals the ABJ anomaly. 
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The quark mass gets both additive and multipUcative renormaUzation. 

The zero of the quark mass is set by Kc in Eq. 10.5. There are two 
ways to calculate Kc at any given a. (i) Assume the chiral relation 
oc niq, calculate the pion mass as a function of 1/2k, and extrap- 
olate it to zero. The value of k at which the pion becomes massless 
is, by definition, Kc- (ii) Calculate the quark mass through the ratio 
{df,A^{x)P{0))/{P{x)P{0)) (based on the axial Ward identity) as a 
function of 1/2k and extrapolate to zero. The two estimates can dif- 
fer by corrections of 0{a). Both these calculations requires a statistical 
average over gauge configurations and an extrapolation in 1/2k. Conse- 
quently, on an individual configuration the zero mode of the matrix M 
will not occur at Kc but will have a distribution about Kc- 

In the full theory detM suppresses configurations with zero modes, how- 
ever, this protection is lost in QQCD. Nevertheless, the problem of zero 
modes at k < Kc on certain configurations, called "exceptional" config- 
urations, needs to be addressed in both theories. 

Near these zero-modes the quark propagator becomes singular even 
though one is in the physical region k < Kc- Consequently, in a large 
enough statistical sample one will hit "exceptional" configurations, and 
hadronic correlations functions on these will show fluctuations by or- 
ders of magnitude (depending on small the quark mass is) from the rest. 
Since this pathology is a feature of Wilson-like actions and not of the 
gauge configurations themselves, one cannot simply remove these con- 
figurations from the statistical sample. The problem of such zero-modes 
increases with decreasing quark mass and (3. Two practical approaches 
used so far have been to (i) simulate at niq > m* where m* is sufficiently 
large such that the spread in the position of the zero modes docs not give 
rise to large fluctuations in the correlation function, and (ii) for fixed 
niq increase the lattice volume as it decreases the probability of hitting 
exceptional configurations. These fixes do not help if one is interested 
in simulating QCD near the chiral limit, i.e., for physical light quarks. 
For this new ideas are needed. A more detailed analysis of this problem 
and a possible solution is given in [59,60]. 

The spin and flavor degrees of freedom are in one-to-one correspondence 
with continuum Dirac fermions. Thus, the construction of interpolating 
field operators is straightforward, for example, tjj'js'ip and t/j-yitp are in- 
terpolating operators for pseudoscalar and vector mesons just as in the 
continuum (see Section 17). 

The Wilson term changes the discretization errors to 0{a). 
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Feynman Rules and details of gauge fixing for the plaquette gauge 

action and Wilson fcrmions are given in [36] . 

Symmetries of Sp: The transformation of the propagator under the 
discrete symmetries is the same as for naive fermions. For a given back- 
ground configuration U these are 



V 
T 
C 

n 



SF{x,y, [U]) ^ 7475 SF{x^\y'\ [U^)) 7574 
Sf{x, y, [U]) ^ 7472 Sj{x, y, [C/^]) 7274 
SF{x,y,[U])^'y5 sUy,x, [f/]) 75 



crn : Sf{x, y, [u]) ^ C74 sUy"^, x^, [C/^]) 74C-1 (1011) 

These relations are very useful in determining the properties of correla- 
tion functions (real versus imaginary, even or odd in p, etc.) as discussed 
in Section 17. 

The dispersion relation, Eq. 10.8, has the following two solutions 



gEa ^ (ma + r)± y^jl + 2mar + m?a?) 

For r = and m small (naive fermion action), the two solutions are 

E = m + 0(771^0^) 

E=-m+ — +Oim'^a^), (10.13) 
a 

where the second solution corresponds to the pole in the Euclidean prop- 
agator at p4 = TT and is the temporal doubler. For r = 1 this doubler is 
removed, as expected, and the single root Ea = log(l +ma) is associated 
with the physical particle. For a general r ^ 1 the second root is given 
by 

E = -m + ^[iT^ + log ^^^^^ + ... . (10.14) 

Such branches which are not associated with the physical state or with 

doublcrs are called ghost branch. I use a very restrictive definition of 
a doubler a solution which has \E\ = 0form=p = 0. Again 
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from Eq. 10.14 it is clear that this ghost branch is pushed to infinity 
for r = 1. In general, if the action couples (n + 1) time slices, then 
the dispersion relation has n branches. The roots can be real, complex- 
conjugate pairs, and/or have an imaginary part in/a as illustrated above. 
Such additional branches can be a nuisance for improved actions, for 
example the 1)234 action discussed in Section 13.2. The goal is to obtain 
a doublcr and ghost free theory This can be done by a combination of 
tuning the parameters (for example choosing r = 1 in the Wilson case) 
and/or going to anisotropic lattices at <C flg as advocated by Lepage 
and collaborators [107]. 

• Conserved Vector Current: The Wilson action is invariant under 
the global transformation defined in Eq. 7.10. To derive the associated 
conserved current (for degenerate masses) we use the standard trick of 
calculating the variation of the action under a space-time dependent 
9{x) and associating the coefficient of 6{x) in SS with 9^ J^. The change 
in the action under such an infinitesimal transformation is 

SS = V(a;)(7^ - r)U^{x)^l){x + jl) exp{i9{x) - i9{x + fi)) 
- k'^iIj{x + iJ,){jn + r) Uj^{x)ip{x) exp{iO{x + fi) - i6{x)) 

which, to first order in 6, is 

-Y. Wi.x){-i^-r)U^{x)^P{x + jl) + 

i>{x + /x)(7^ + r) Ul{x)^{x)] [i—] . 

The conserved current, obtained after integration by parts, is 

= V^(a;)(7^ - r)U^{x)^{x + A) + i^{x + + r) Ul{x)^{x) . 

(10.15) 

is hermitian and reduces to the symmetrized version of the 1-link 
vector current for r = 0. In many applications like decay constants one 
uses the local (flavor) currents defined by 

2y^(a;) = V^(a;)7^y [/^(a;)V(a; + fi) + h.c. 

2Af,{x) = ip{x)j5jfj, — Ufj,{x)ilj{x + fi)+ h.c. , 



(10.16) 
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which are not conserved, and consequently have associated non-trivial 
renormalization factors Zy and Za which have to included when calcu- 
lating matrix elements. 

• The chiral condensate: As explained in [57], the chiral condensate 
(il^ip) is not simply related to the trace of the Wilson quark propaga- 
tor (Si? (0,0)). The breaking of chiral symmetry by the r term intro- 
duces contact terms that need to be subtracted non-perturbatively from 
Sf{0,0). This has not proven practical. The methods of choice are to 
use either the continuum Ward Identity 



{^p^}"^^ = {0\Sf{0,0)\0) = lim mg d'^x{0\P{x)P{0)\0), (10.17) 



or the lattice version of the Gell-Mann, Oakes, Renner relation [61] 



A comparison of the efficacy of the two methods is discussed in [62] . 

• Operator Mixing: There is, in general, a mixing of operators with 
different chirality in the Callen Symanzik equations. This is true even 
at Kc and is a consequence of the explicit breaking of chiral symmetry. 
This mixing of operators poses serious problems in the the calciilation 
of matrix elements of the weak interactions Hamiltonian, i.e., the lattice 
analyses of Ward identities and obtaining amplitudes with the correct 
chiral behavior is more involved. For example the chiral behavior of 
matrix elements of the tree-level 4-fermion operators (relevant to the 
extraction of Bk, Br, and Bs) is not the same as in the continuum 
theory. The reason being that the Ward Identities associated with the 
spontaneously broken ^[/^(n/) receive corrections at 0{a) in the lattice 
theory. This issue will be discussed in more detail by Prof. Martinelli. 

In Section 13 I will discuss how the 0{a) errors in WF can, with very 
little extra effort, be reduced to O(a^) or higher. For the time being let 
mc use this drawback of Wilson fermions to motivate the formulation of 
staggered fermions. This formulation has 0{a^) errors, retains enough 
chiral symmetry to give the correct chiral behavior for amplitudes involving 
Goldstone pions, but at the expense of a 4- fold increase in the number of 
flavors. 



GMOR 




(10.18) 



INTRODUCTION TO LATTICE QCD 



53 



11. Staggered fermions (SF) 

The 16-fold doubling problem of the naive fermion action given in Eq. 7.6 

can be reduced to 4 by the trick of spin-diagonalization [21]. Staggered 
fermions x s-re defined by the transformation 

V(x) = r,xW V^(x) = x(x)rt = 7r72"^7f 7r- (H-i) 

In terms of x the action can be written as [21] 

Ss = rUg T,xX{x)x{x) + lT,x,nXxVx,n{U^,xXx+n - Ulx_p,Xx-ix) 

^E.,yX{x)Mgx{y) (11-2) 
with the matrix given by 

M'^[U]x,y = mg6^y + ^ XI ^=='1^ [^^^-A,!/-/* - -(ll-S) 

ft 

The 7 matrices are replaced by the phases rix,ij,- It is convenient to define 
the following phase factors 

Vx,^ = i-)^''<>'"'' 

g _ ^„'J2;i+X2+X3+X4 

Sr{x) = ^ (1 + Ci^C - VnV^ + Cu-C^Vi^riu) ■ (11.4) 

From Eqs. 11.2,11.3 it is easy to sec that the different spin components of 
X are decoupled as the phase factor rjx,^ depends only on the site index and 
direction and do not have a spinor index. One can therefore drop the spin 
index on x leaving only color degrees of freedom at each site. This reduces 
the original 2''-fold degeneracy of naive fermions by a factor of four. The 
mass term in Eq. 11.3 is hermition, while the Ip term is anti-hermitian, 
which as I discussed before is required to realize chiral symmetry. 

The SF action has translation invariance under shifts by 2a due to the 
phase factors r]x,n- Thus, in the continuum limit, a 2^ hypercube is mapped 
to a single point, and the 16 degrees of freedom reduce to 4 copies of Dirac 
fermions. That is, for every physical flavor, the staggered discretization 
produces a 4-fold degeneracy denoted as the staggered flavor. At finite a 
the gauge interactions break this flavor symmetry and the 16 degrees of 
freedom in the hypercube are a mixture of spin and flavor. This is one 
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of the major drawbacks of staggered fcrmions. Construction of operators 
and their interpretation in terms of spin and flavor is non-trivial as I show 
below for the case of bilinears. I begin the discussion of staggered fermions 
by listing its symmetries. I will then outline the construction of bilinear 
operators and end with a discussion of Ward identities. 
Symmetries of SF: The analogue of 75 invariance for SF is 



e^Ml{x,y)ey = Ms{x,y), 
<^xSp{x,y)€y = SF{y,x) . 
and the symmetries of this action are [63] 



Translations 

Rotations 

Inversion 

Charge Conjugation 

Vector(Baryon#) 

Axial 



I 

C : 
Uil)v : 
U{1)a : 



X(x) 
Xix) ■ 



(11.5) 



Cf^ix)xix + ^J■) 

e{x)x{x), 
-e(a;)x(a;) 

e^'xix), 



X{x) 



where the phases ri,(,e and Sr are defined in Eq. 11.4. It is the last 

invariance, the U{1)a, or also called the C/(l)t invariance, that becomes a 
flavor non-singlet axial symmetry in the continuum limit, and is the main 
advantage of staggered fermions. 

Operators with Staggered fermions: There are two equivalent ways 
of understanding flavor identification of staggered fermions: (i) momentum 
space [64], and (ii) position space [65,66]. Here I review the position space 
approach. 

For the purposes of constructing hadronic operators it is convenient to 
represent the 16 degrees of freedom in a hypercube by a hypercube field Q, 
which in the continuum limit represents the four spin components of four 
degenerate flavors [67] . This is done by dividing the lattice into hypercubes 
identified by a 4- vector y, i.e. y is defined on a lattice of spacing 2a, and 
the set of 16 vectors {A} to points within the hypercube. These are called 
hypercube vectors, and addition of two such vectors is defined modulo 2, 
= {Af^ + Bfj)mod2- Then, in the free field limit (setting gauge links 
U'x.ii = 1 with the understanding that for the interacting theory non-local 
operators are to be connected by a path ordered product of links that 
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makes them gauge invariant, or these operators are to be evaluated in a 
fixed gauge) 

W = -^J2^AX{y + A) (11.6) 

^ {A} 

where Nf = 4 and 

= it'li'li'lt'- (11.7) 
A general bilinear operators can be written as 

OsF = Tt{QjsQiI) 

= ^x{y + A)Tr(r^7srB^F)x(y + b) 

A,B 

= ^ x(y + ^)(7s o ^F)ABx{y + B) 

A,B 

= QilS<^^F)Q. (11.8) 

where the matrices 75 determines the spin and S^p the flavor of the bilinear. 
Each of these can be one of the standard sixteen elements of the Clifford 
algebra. In the last expression, the field Q is expressed as a 16 component 
vector. It is easy to see from Eq. 11.8 that local operators are given by 
7S = ^_F. The identification of meson operators is as follows [65,68,69]. 

= 0(7V75 ^ ^f)Q 

P ^Q{l5®iF)Q 

S =Q{l®iF)Q. (11.9) 

In the continuum the mesons lie in the flavor SU(4) representations 1 and 
15, while on the lattice these representations break into many smaller ones. 
For example, the continuum 15-plet of pions breaks into seven lattice rep- 
resentations (four 3-d and three 1-d). 

15 (75 (g) 7i) ® (75 (g) 7i74) ® (75 O 7i75) ® (75 ® 7i7475) © 

(75 O 74) © (75 ® 75) ® (75 O 7475) , 

1 ^ (75 1) • (11.10) 

In the first line, « = 1,2,3, so the representations are 3-dimensional. The 

other representations are 1-dimensional. Similarly, the continuum rho 15- 
plet (times three spin components) breaks into eleven representations (four 
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6-d and seven 3-d). 

15 {ji O 7j) e (7i O IjH) ® in O 7j75) ® (7i O 7j7475) ® 
(7i ® 7i) ® (7i O 7i74) ® (7i ® 7i75) ® (7i ® 7i7475) ® 
(7i O 74) ® in O 75) ® (7i O 7475) , 
1 ^ (7i«'l) • (11-11) 

Of the sixteen pions, the Goldstone pion has flavor ^5. The remaining 15 
pions are given by the other 15 members of the Chfford algebra and are 
heavier due to the staggered flavor symmetry breaking, and their masses 
do not vanish at rriq = 0. In fact the mass for each representation can be 
different. The degree to which the staggered flavor symmetry is broken can 
be quantified, for example, by comparing the mass of the Goldstone pion 
versus those in the other representations. 

The above examples make it clear that the mixing of spin and flavor 
complicates the analysis and is a major disadvantage of staggered fermions. 
I will not go into details of calculations using staggered fermions. A useful 
list of references is as follows. The construction of baryon operators is 
given in [70] . For an illustrative calculations of the spectrum see [69] . The 
transcription of 4-fermion operators arising in effective weak Hamiltonian 
and their matrix elements within hadronic states are discussed in [63,71]. 
Chiral Symmetry and Ward Identities: The staggered fermion action 
can be written in terms of the 16 component fields Q, defined in Eq. 11.6, 
as [4] 

'5s = TOg^0(y)(i® i)g(y) 

y 

+ EW[(7M®lR + (75®CM?5)nM]Q(j/). (11.12) 

This action, for rUq = 0, is invariant under the following abelian transfor- 
mation 

Q{y) ^ e*''(^=««=)0(y) 

Q{y) - 0(y)e^'(^=®«'^) . (11.13) 

Thus one preserves a non-trivial piece of the original ?7(4) x U (4) symmetry 
(called U{l)f or U{1)a) whose generator is 75 (S> £,5 and the associated 
Goldstone boson operator is (5(75 <?> ^5)<5- This axial U{1)a symmetry, 
along with the U{l)v, allows, for appropriate transcription of operators, 
Ward Identities for lattice amplitudes similar to those in the continuum. 
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The best studied example of this is the lattice calculation of the K^T{^ 
mixing parameter Bk [72]. 

The basic tools for constructing the Ward Identities are exemplified by 
the following two relations involving qnark propagators and the zero 4- 
momentum insertion of the pseudo-Goldstone operator x(a^)x(a;)(— 1)"^ and 
the scalar density x{^)xi^) [63] 

{-ir'G2{xf,Xi) +Gi{xf,x,){-ir' = 

(mi +m2)^Gi(x/,a;)(-l)"G2(:c,a;i) (11.14) 

X 

G2{xf,Xi) - Gi{xf,xi) = 

{mi - m2)^Gi{xf,x)G2{x,Xi) (11.15) 

X 

where Gi{xf,Xi) is the propagator from Xi to Xf for a quark of mass of 
mass mj. The derivation of these identities, using the hopping parameter 
expansion, is given in Ref. [63]. Note that these relations are valid for all 
a and at finite m. An example of the usefulness of these relations is in 
extracting the chiral condensate as shown below. Eq. 11.14 gives 

(XX) = G{xi, Xi) = m 5^(-l)--*G(a;i, ar)(-l)-G(a;, Xi) 

X 

= mY,\G{xi,x)\'', (11.16) 

X 

where the last term on the right is the zero 4-momentum sum of the pion 
correlator. Similarly Eq. 11.15 gives 



dm 



J2 Gixi, x)Gix, Xi) = J2 \Gi^i^ a;)p(-l)---* . (11.17) 



Combining these two equations one can extract the value of the condensate 
in the chiral limit 

(Xx)(m = 0)= (1-m— )(xx)- E \Gi{xi,x)\^ (11.18) 

even 

by summing the Goldstone pion's 2-point function over sites with even 

(x - Xj). 

11.1. Summary 

Even though most of the simulations done until recently have used cither 
the staggered or Wilson formulation of fermions, neither is satisfactory. 
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They have their relative advantages and disadvantages which have been 
evaluated on a case by case basis. The rule of thumb is that in the ex- 
traction of observables where chiral symmetry plays an essential role and 
the external states are Goldstone bosons, staggered fermions do better. 
Otherwise Wilson fermions are preferred due to their correspondence with 
Dirac fermions in terms of spin and flavor. The focus of current research is 
to improve the Wilson formulation and thus have the best of both worlds. 
This is being done in a variety of ways as discussed in Section 13. 



12. Improved gauge actions 

There have been a number of methods proposed to improve the gauge 
action beyond 0{a^). These include the Symanzik approach, mean-field 
improvement, and renormalization group inspired perfect actions. Let me 
begin this discussion by reviewing the Symanzik program. 

12.1. Liischer-Weisz Action 

The leading order term in the expansion of all Wilson loops is O^^^ = 
'^fj^^ Ffj.i^Ffiv and corrections begin at 0{a^) as there are no dimension 
5 operators. Thus, any lattice action written as a linear combination of 
Wilson loops will have the correct continuum limit with corrections at 
O(a^). There are three dimension 6 operators, which in continuum notation 
are: 



of = ^ Trf^^i^^p^.i^^pj . (12.1) 

Also there are only 3 six-link loops that one can draw on the lattice. These, 
as shown in Fig. 7, arc the planar df\ twisted £3^' and the L shaped 
£3^'' respectively. Thus, classical improvement of the lattice action, i.e. 
removing the O(a^) term, can be achieved by taking a linear combination 
of the plaquette and these three six-link loops. Each of these loops has the 
expansion 

C = ^(4) 0(4) + ^(6) 0(6) ^ ^(6) 0(6) ^ ^(6) 0(6) + . . _ (12.2) 
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^ 1 



(6) 



(6) 



Fig. 7. The three types of six link loops. 



Table 4 

The coefficients of the continuum opera- 
tors of dimension 4 and 6 in the classical 
expansion of Wilson loops with 4 and 6 
links. 
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and Liischer and Weisz have shown an elegent way of calculating the ex- 
pansion coefficients [73]. Their results are summarized in Table 4. 
The lattice gauge action can be written as 

= iU'\9')C^''+j:'^^'(s')^^'] (12-3) 

in terms of the plaquette and these three 6-hnk loops. Using the results in 
Table 4 one gets the normalization condition 

cW(5^) + 8c^^\g^) + 8cf + 16cf (5^) = 1 (12.4) 

such that in the continuum limit the action reduces to -F)i,yF^y/4. Also, 
from Table 4 one can see that tree level improvement can be obtained by 
the choice 

c(4) + 20cl') = 0; = 0; cf = . (12.5) 
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This removes Ol by cancelling the contribution from the two planar loops. 
Of the three dimension six operators, oj^^' alone breaks Lorentz (rotational) 
invariance. Thus, rotational invariance of the heavy qq potential can be 
used to non-perturbatively tune c^^ and to compare the various improved 
actions as discussed later. 

To include quantum effects. Curd, Menotti, and Paffuti did a 1-loop 
calculation. They imposed the criteria that all Wilson loops be improved 
at the leading logarithm level [74,75]. This corresponds to improving the 
string tension. (To kill the O(a^) terms at all length scales instead of only in 
the long distance heavy quark potential requires including the finite pieces 
in the 1-loop calcTilations. The fully improved action, as we shall sec in the 
final Lilscher-Weisz form, will lead to 0{g'^) corrections to c^^^ and c'f^ and 
the need for 0^2^.) The result of their calculation required the improved 
action to have the same form as given in Eq. 12.5. This agreement with the 
classical result is fortuitous. Symanzik's on-shell improvement program, in 
general, prescribes that one has to derive three relations to fix the cf'^ such 
that on-shell quantities do not have O(a^) corrections. Liischer and Weisz, 
in fact, found that improving all large loops produced only two constraints 
at tree-level [76] 



Thus, they were left with a 1-parameter family of improved actions. The 
reason for this freedom is the presence of a redundant operator O^^ gener- 
ated by the field redefinition + a'^g^d{g)D,^F,^^ [73,76]. Since one 
has the liberty to choose the coefhcient of an interaction that is part of a 
redundant operator without affecting on-shell improvement, the convenient 
choice C3^^ = simplifies the action and again leads us back to Eq. 12.5. To 
summarize, Eqs. 12.6 are the most general set of constraints that can be 
derived based on improving large Wilson loops. The criteria employed by 
Curci et.al, improving the large distance part of qq potential, is a subset. 
Since this involves a spectral quantity one can use the redundant operator 
to get the simpler form Eq. 12.5. 

The relation cf^ ~ holds both at tree and 1-loop level as it comes 
from exploiting the redundant operator. Combining this with the first 
constraint in Eq. 12.6 and the normalization condition Eq. 12.4 leads to a 
single condition, which at tree-level is 



c(4) + 204'^)-44'^)+4cf =0 



(12.6) 




,.(fi) L 

- 12' 



(12.7) 
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and is all that can be gotten from Wilson loops. To get a second condition 
Liischer and Weisz had to set up other spectral quantities to be improved. 
Failing to find a suitable condition in the real theory, they looked at a 
situation in which the gluons are massive and appear as asymptotic states. 
Using the mass and scattering amplitude of these states they derived the 
1-loop improved action [76], which goes under the name Liischer-Weisz 
action and has leading correction at 0{g'^a^), 

cf\9') = ~-0.02521g' 

(5') = -0.0044l52 
cf(5^)=0. (12.8) 

In this form for the Cj, the normalization condition, Eq. 12.4, is already 

incorporated. 

The bottom line on the Liischer-Weisz action is that since the coefficient 
is small, the action, to a very good approximation, can be improved by 
keeping just and C^^^ as in the classical case. This 1-loop result can 
be mean-field improved without much effort as discussed next. 



12.2. Tadpole/Mean- field Improvement 

The representation of gauge fields by the exponentiated form 

2 2 

C/^W = e^»sA,(.) ^ 1 ^ iagA^ix) - ^A^x) + ... (12.9) 

gives rise to local quark-gluon vertices with one, two, . . . gluons. Of these 

vertices, all but the lowest arc lattice artifacts. Contracting the two glu- 
ons in a^(/^A^(a:)/2 gives rise to tadpole diagrams. Parisi, and Lepage 
and Mackenzie [77,78] recognized that these artifacts, though naively sup- 
pressed by powers of a, arc in fact only suppressed by powers of due to 
the ultraviolet divergences generated by the tadpole loops. They reasoned 
that these tadpole contributions were responsible for (i) the poor match 
between short distance quantities (like the link in Landau gauge, small 
Wilson loops, Kc, ■ ■ .) and their perturbative estimates, and (ii) large coef- 
ficients in the perturbative lattice expansions. They, therefore, proposed a 
mean- field improvement (also called tadpole improvement (TI)) to get rid 
of these artifacts. 
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The mean-field improvement starts by assuming that the lattice fields 
can be split into a UV and IR parts, and that the UV part should be 
integrated out [77,78], i.e. 

^ragA^{x) ^ gms(4«(x)+Ajf^(a;)) ^ yQe*"sA;,«(a.) ^ uoU^{x) . (12.10) 

The rescaling of links by an overall constant factor uq < 1 keeps the theory 
gauge invariant. Under this rescaling (each U uqU) the Wilson gauge 
and fermion action becomes 



-.w 



E 



/3uq Wii + V'V' ~ KUotpDtp 



(12.11) 



Thus, a simplified recipe for mean-field improvement is to use the bet- 
ter expansion parameters = g^/u^, k = kuq, U = U/uq, and ip = 
'ipcont I ^2kuq. Note that in terms of the tilde variables the theory looks 
the same as the original theory except that all links are scaled by uq. 
Lepage and Mackenzie [78], and subsequent work, have shown that mean- 
field improvement of 1-loop estimates gives, in many cases, results that are 
very close to non-perturbative estimates, and the expansion coefficients are 
much smaller. Tadpole improvement is, therefore, a useful thing to do, es- 
pecially since it requires negligible extra effort once the 1-loop calculation 
has been done. The only limitation is that it does not work equally well in 
all cases and there is no real understanding of why it docs not work well 
in some cases. To illustration tadpole improvement I will apply it to the 
Liischer-Weisz action. 

There are two common choices of the tadpole factor uq : (i) the fourth 
root of the plaquette, and (ii) the expectation value of the link in Landau 
gauge. So far most calculation have used the first choice, whereas Lepage 
and collaborators now advocate the link in Landau gauge as a better choice 
[79]. To implement tadpole improvement requires both the 1-loop expres- 
sions for these variables and their non-perturbative values. For the Wilson 
action 

uo = {U)lg = 1 - 0.0772^2 (12.12) 
and for the Liischer-Weisz action, Eq. 12.5, it is [80] 

u, = (Wn)^/^ = 1 - (i)(^)(^^^^)5^ = 1 - 0.06104V • (12-13) 
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Denoting the 1-loop tadpole factor as uq = 1 — Xg^, the tadpole- 
improved Liischer-Weisz (TILW) action is [81] 



The non-perturbative value of uq has to be determined self-consistently as 
the relative co-efficients in the action themselves depend on uq. Note that 
which comes in only as a quantum effect, has, after TI, a coefficient 
is reduced to a negligible value. 

An important point to note about the TILW action is the strength of 
Cf\ Since uq = 0.85 — 0.90 for lattice parameters currently in use, the 
strength of C^^^ is larger than the classical value 1/20, i.e tadpole im- 
provement suggests an over-correction. Later on we shall see that this 
over-correction is also characteristic of improved actions based on renor- 
malization group methods. 

A simplified description of TI applied to the above example is that the 
coefficient of C^^^ was multiplied by v^{pert) / v^{non — pert), the rela- 
tive factor of two extra links between the plaquette and Cf ^ . The factor 
UQipert), which differs from UQ{non — pert) at 0{g'^), is absorbed in the 
1-loop expansion. This, as is obvious, makes the coefficient of small. 
The correction is, therefore, dominated by uo{non — pert) which can be 
measured with very high precision. There, however, is the residual uncer- 
tainty due to the freedom in the choice of uq. To summarize, the imderlying 
philosophy of TI is that the large tadpole contributions are present in the 
perturbative expansion of all quantities, and can be absorbed by factoring 
UQ^pert) raised to the appropriate power depending on the number of links. 
This makes the expansion coefficients small. The modified perturbative ex- 
pansion is then in terms of an improved = g^/u^, and the correction is 
dominated by overall factors of UQ{non — pert), which should be chosen to 
be the expectation value of the link in Landau gauge. 

One of the first tests to convincingly demonstrate improvement using 
the TILW action was restoration of rotational invariance in the qq poten- 
tial. The results from [82] are reproduced in Fig. 8. Since then tests of 
improvement using the scaling of the finite temperature transition, Tc, and 
topological susceptibility have been carried out. These results, as will be 
discussed later, all make a strong case for improvement. 





(12.14) 
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Fig. 8. Static-quark potential computed on 6^ lattices with a ^ 0.4 fm using the Wilson 
action and the TILW action. The dotted line is the standard infrared parameterization 
for the continuum potential, V(r) = Kr — 7r/12r + c, adjusted to fit the on-axis values 
of the potential. Figure reproduced from [82]. 
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12.3. The Renormalization Group and the Renormalized Trajectory 

This is a brief summary of Wilson's formulation of the renormalization 

group and the concept of the renormalized trajectory. For details see [83- 
90]. I will start with a discussion in momentum space. Consider the effec- 
tive theory formulated at some scale Q in terms of the generalized action 

S{Q)=^ K^^\Q) So. (12.15) 

a 

where 5^ are all possible operators (interactions) and Ka^ are the cor- 
responding coupling constants. Now consider an exact renormalization 
group transformation (RGT) Tf, that integrates all momenta between Q 
and Q/b= Q/2. Here 6 = 2 is the scale factor of the transformation. The 
renormalized theory at scale Q/2 is 

SiQ/2) = Y, K^a\Qm (12.16) 

a 

where K^^ are analytical functions of 

/^«=T6[i^(0)] . (12.17) 

The sequence of theories defined by repeated use of % defines a flow tra- 
jectory in the space of couplings. The physical theory is preserved under 
each Tf,, and thus along the whole trajectory. Correlation lengths, mea- 
sured in lattice units, change trivially ^/a — > ^/ab. Fixed points of this 
transformation are defined by 

i^W=Tb[i^W], (12.18) 

i.e. points at which the theory reproduces itself under 71. Since all corre- 
lation lengths (momenta) scale by b under %, therefore at the fixed points 

£| . , = or oo. (12.19) 

^ I fixed point ^ 

Fixed points with ^ = arc trivial fixed points. Examples of these are 
the T = and T = oo limit of most statistical mechanics models. We 
are interested in non-trival fixed points with ^ = oo. At these points the 
correlation length ^/a = l/Ma, measured in imits of lattice spacing, is 
infinite. This allows us to take the continuum limit of lattice field theories 
— set a = holding M fixed. 

Now consider a point {Kp} close to the fixed point {K'^} and a RGT with 
an infinitesimal b. Under the assumption that there are no singularities in 
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the space of coupling constants, the new couphngs K'^{Kf^) generated by 
Tb can be written as a Taylor expansion about the fixed point 

dK' 

K{K0) = K*) + ^^AK0 + ... . (12.20) 

In that case one can define a "linear region" in the neighborhood of {K'^} 
in which 



K' 



AK(3 (12.21) 



where AKp = Kp — K"^. The linearized transformation matrix (also called 
the stability matrix) 



(12.22) 



controls the RG flows locally. In field theory discussions, the stability 
matrix Tap is called the /9-function. Note that, because of the truncation 
of Eq. 12.20 all statements based on the analysis of Tap have the implicit 
assumption that they are approximately true in the linear region close to 
the fixed point and for small &, and become exact only as {Kp] — » {K*p}. 

Wilson formulated the real space renormalization group as an alternate 
to the momentum space renormalization group discussed above. This ap- 
proach is characterized by a constrained integral over the position depen- 
dent original degrees of freedom s to produce an effective theory, 



JvsT{s',s)e^''^''^'K (12.23) 



Here T(s', s) is the constraint that specifies how the s in some local neigh- 
borhood (called the block cell) are to be averaged and replaced by an 
effective s'. A necessary condition on T{s',s), that preserves the partition 
function, is 



/ 



Vs' T{s',s) = 1. (12.24) 



One can see that if the block cell in f— dimensions is 6^, and the con- 
straint T{s' , s) is equivalent to integrating momenta in the range Q/2 — Q, 
then the two approaches are equivalent. In general, an exact integration 
of the upper half momenta is not feasible and T{s', s) is constructed as 
some sensible "average" of the original spins in the block cell. Other than 
the condition Eq. 12.24, the definition of a sensible "average" is that it 
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preserves the physical properties of the theory. Because of this dual de- 
scription, momentum versus real space formulations, the RGT is also called 
a blocking transformation a usage more common in statistical mechan- 
ics and numerical implementations. 

A brief summary of the properties of the matrix 7^^, evaluated using 
either the momentum or real space formulations, and its associated eigen- 
values Aj and eigenvectors $8 obtained by solving = Aj$i is as follows. 

1. Tajs is, in general, not a symmetric matrix. Its eigenvalues are, however, 
real. 

2. The eigenvalues depend on the scale factor b. Since the theory obtained 
after two successive transformations, 6 (8) 6, should be the same as the 
one obtained after a single scale change by 6^, the general form of the 
eigenvalues is A; = 6^* where the Aj are independent of % and are called 

the critical exponents. 

3. The eigenvectors are by definition the linear scaling fields. Consider 
a small deviation of size u from the fixed point along a given eigenvector 

Then under a RGT, the renormalized value of the coupling, u', is 
given by 

= TbU^ = 6^*u$ . (12.25) 

Note the crucial use of the assumption of linearity in deducing this simple 
scaling of u. The eigenvectors depend on the precise form of %. 

4. Eigenvalues Aj > 1 lead to flows away from the fixed point as deviations 
along these $i grow under Tf,. These eigenvalues and the corresponding 
eigenvectors are called relevant. 

5. Operators corresponding to eigenvalues A^ < 1 die out geometrically 
with the number of blocking steps. These operators do not contribute to 
the long distance properties of the system and are consequently called 
irrelevant. The associated exponents control corrections to scaling. 

6. Eigenvalues that are exactly unity are called marginal. Deviations 
along marginal directions do not change under RGT. To ascertain that 
an eigenvalue is truly marginal, one has to go beyond the linear approxi- 
mation. Typically marginal operators develop logarithmic corrections at 
higher order. A quick examination of the /3- function of asymptotically 
free field theories should convince you that g is a marginal coupling at 
tree-level and becomes relevant due to quantum corrections. 
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7. There is an additional class of eigenvectors, called redundant opera- 
tors, that are not physical. Redundant operators are manifestation of 
the invariance of the action under a redefinition of the fields, i.e., they 
define relations between the different operators such as the equations of 
motion. The corresponding eigenvalues can be < 1, = 1, or > 1. These 
eigenvalues depend on the choice of the RGT, i.e. different RGT with 
the same or different scale factor b will give different A^. Since different 
non-linear RGT correspond to different rescaling of the fields, redundant 
operators are specific to the transformation. In fact, one way to isolate 
them is to repeat the calculation of Tap with a different RGT: the expo- 
nents Aj that change with the RGT correspond to redundant operators. 
Even though they carry no physics information it is important to isolate 
these operators, especially in numerical calculations. For example, when 
improving the action one wants to use all possible redundant operators 
to simplify the action and thus minimize the computational cost. 

8. Under certain conditions Tap is block diagonal, i.e. the couplings break 
up into sub-sets that are closed under the RGT. For example, the spin 
flip symmetry s — > — s of the Ising Hamiltonian at zero magnetic field 
causes the interactions to break up into odd and even sectors. One 
can, therefore, calculate Tap for each sector separately as the eigenvalue 
spectrum is independent. 

In this generalized space of couplings {Ka}, there exist points, called 

critical points {K^}, with ^ = oo. Starting from a critical point, a RGT 
produces another since ^' = = oo. Thus the set of critical points define 
a hypersurface in the infinite dimensional space {Ka}. The RG flows on 
this surface can (a) meander randomly, (b) go to some limit cycle or a 
strange attractor, or (c) converge to a fixed point {K^}. We shall only 
be concerned with possibility (c). At the fixed points the renormalized 
coiiplings are exactly equal to the original couplings, K'a = Ka, and the 
theory reproduces itself at all length scales. This is in distinction to critical 
points where only the long distance behavior is reproduced. The location of 
the fixed point depends on the RGT and in general a given RGT will have 
more than one fixed point with ^ = oo. Each such fixed point has a basin 
of attraction i.e. the set of {Ka} that converge to it under the RGT. This 
hypersurface is thus orthogonal to all relevant eigenvectors of 7^. The basin 
of attraction defines the universality class since the long distance behavior 
of all theories is the same and governed by the fixed point theory. 

The flow along a relevant scaling field, whose end-point is the fixed 
point is called the renormalized trajectory (RT). Just as the fixed point 
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is an attractor of all critical points in its universality class, the RT is the 
attractor of all flows terminating in the critical points that lie in the basin 
of attraction of the fixed point. This is illustrated in Fig. 9. For simplicity 
I have assumed that there is only one relevant coupling, as is the case for 
QCD. Along the RT there are no scaling violations as, by construction, 
all the irrelevant scaling fields are zero. Thus, simulations done using an 
action along the exact RT, will reproduce the continuum physics without 
discretization errors. The RT is thus the holy grail of attempts to improve 
the lattice formulation, i.e. the perfect action. 

Simulations on finite lattices arc done in a region of coupling space where 
all ^ <C and not at the fixed point or at the critical points. Current typ- 
ical values lie in the range ^/a = = 5 — 15. Even for quenched 
simulations this range will not change in the near future. Points in this 
region lie on flows that terminate at critical points that may lie quite far 
from the fixed point. Corrections to scaling along these flows at ^ « 10 
may, therefore, be large. One can reduce these corrections by adjusting the 
action to lie on another flow that starts closer to the fixed point and at 
^ « 10 is closer to the RT. Thus, to improve the action a non-perturbative 
approach is to estimate the RT. Since flows are attracted by the RT, its 
location can be estimated by starting with even the simple plaquette ac- 
tion and following the flows generated under some blocking transformation. 
Evaluating the blocked action along the flow is therefore the basis of all 
renormalization group methods for improving the action. Methods used 
for determining the blocked action are reviewed in [86]. 

A standard assumption underlying the validity of the renormalization 
group is that the flxed point action is local (the strength of the couplings 
fall off exponentially with the size of the Wilson loops and the separation 
between points that are coupled by the Dirac operator), nevertheless, a lo- 
cal action still involves an "inflnite" number of couplings. Keeping all these 
when evaluating the action generated under blocking, and in simulations, 
is impractical. So one is forced to to tnmcate the action to a finite (small) 
number of interactions. Also, the RT at ^ ~ 10 may be sufficiently different 
in form from that near the fixed point. Consequently, the questions one has 
to face when designing improved actions are (i) what is the minimal set of 
terms to keep, (ii) are all the redundant interactions removed from this set, 
and (iii) how to compare such improved actions evaluated using different 
RGT. Also note that the projection of the RT on to a truncated subspace 
may not be the action with the smallest scaling violations that one can 
design in that space. I shall explore attempts to answer these questions 
next. 
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Fig. 9. A schematic of the critical hypersurface = 00), the fixed point, the renormal- 
ized trajectory (RT) flowing towards the infrared under blocking transformations, and 
flows attracted by them. Also shown is the simple Wilson action and flows starting from 
it. 



1 will summarize current activity in designing improved gauge actions 
by reviewing those that are being used to calculate phenomenologically 
relevant quantities. Unfortunately, there is no standard notation used by 
different collaborations. To facilitate comparison I will therefore try to 
write them in a common form as well as maintain contact with the original 
notation. Wilson loops will be denoted by W^^^ where the index i will be 
used for both the type of loop and its representation. Unless otherwise 
stated, it will be implicit that the real part of the trace of Wilson loops 
is to be taken. The W^^ are normalized such that the maximum value of 
the trace is unity. As before (3 is the overall coupling and the Ci give the 
relative strengths of different loops. 

12.4^. Iwasaki Action 

Iwasaki studied block transformations in perturbation theory. His criterion 
was to obtain a gauge action that after a few blockings comes close to the 
RT. His Ansdtz is [91] 




(12.26) 
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where W^^^ is the rectangular loop C\ shown in Fig. 7. (Iwasaki used the 
notation (1 - 8co)W^^^ + CqW^^'^ with Cq = -0.331 which is the same as 
Eq. 12.26 with a rescaled (3.) This action shows better rotational invariance 
and scaling of Tc/^/a. It is being pursued by the JLQCD collaboration [92]. 

12.5. b = \/3 blocking transformation 

The -s/S blocking transformation has certain very nice properties for 4- 
dimensional gauge theories in addition to being a smaller scaling factor 
than the conventional 6 = 2 [89]. The block cell consists of a site and 
its eight neighbors. As a result the transformation has much higher rota- 
tional symmetry for both fermions and gauge fields [89]. Gupta and Patel 
have studied the flow under blocking using Monte Carlo Rcnormalization 
Group methods. They investigated the effect of (i) truncation and (ii) vari- 
ation in the Ci as a function of /3 [87] . For phenomenological studies they 
approximate the RT by a linear trajectory in a four parameter space 

= p Re{ 4>^i^TVW^(i><i) 

+ 4^i^[^(TrT4^(i><i))2-iTrM/(i><i)] 

+ c^^El^lTrW^^^^^^r-^] }, (12-27) 

with 

cixV4>^i = -0.12; 4^Vcf ' = -0.12; 4>^V4>^^ = -0.04. (12.28) 

This particular choice of couplings (truncation of the action) was made 
for the following reasons. The coupling for the plaquette in the adjoint 
representation, Cg^^, was included to avoid the phase structure in the 
fundamental-adjoint plane (see Section 15); that in the 6-dimcnsional rep- 
resentation, CgX^, since it differs from Cg^^ only at 0{a'^) and higher, and 
the 1x2 rectangle in the fundamental representation, c]^^, as it is the 
simplest loop that kills the O(a^) terms a la Symanzik. 

Classical expansion, see Eq. 7.3, of the loops gives the normalization 
condition 

+ Scp""^ = 1 , (12.29) 

where I have lumped together the first three terms coming from W^^^^^ 
as the effective contribution of the plaquette. Using the ratios given in 
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Eq. 12.28, one gets a very high weight c^^^ = 1/0.11, which as shown in [51] 
is correlated with a suppression of dislocations versus physical instantons. 
Similarly, the coefficient of the leading discretization correction of^ , using 
Table 4, is 



1 

24 



1X1 I 5 1X1 I 9 ly, 



5 1x2 

6 



F 



(12.30) 



For the ratios given in Eq. 12.28 this is not zero, but over-corrected by 
the proposed strength of W^^^'^\ At 1/a w 2 GeV, this over-correction is 
larger than that for the tadpole improved Liischer-Weisz Ansdtz, similar to 
Siwasaki, and smaller than that in Sqcdtaro- 

Early spectrum studies and measurement of the glueball spectrum [93] 
to test improvement were inconclusive. In hindsight the failure to quantify 
improvement is due to poor statistics and not looking at quantities like 
rotational invariance at strong coupling, where the lattice artifacts show 
up clearly for the simple Wilson action. Recently this action has also been 
used to study lattice topology and the results show that short distance 
artifacts are suppressed [51]. 

12.6. 6 = 2 blocking transformation 



The QCDTARO collaboration has also investigated the flow of the action 
under blocking [94,95]. They use a 6 = 2 blocking transformation and eval- 
uate the couplings on blocked lattices using Schwinger-Dyson equations. 
Their overall conclusion is that the blocked action, to a very good approx- 
imation, can be represented by the two couplings W^^^^ and W^^^. Their 
estimate of the RT cannot be approximated by a linear relation between 
the coefficients. An approximate Ansdtz in the range of couplings simulated 
is 

Sqcdtaro = E E E (K^' ' ^-^^K^^') (12-31) 

Note that the coefficient of VF^^^ is larger than the tadpole improved 
Luscher-Weisz, Iwasaki, or the 6 = \/3 actions. 

12.1. Classically Perfect Action 

Hasenfratz and Niedermayer showed that one can derive the classically 
perfect action for asymptotically free-field theories using a saddle-point in- 
tegration of the RGT about .g = [96,97]. They use the 6 = 2 blocking 
transformation proposed by Swendsen [98]. The practical issue here is 
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Table 5 

The current coefficients of the truncated "perfect" action with 4 and 
6 Unk loops. 



Loop 


Cl 


C2 


C3 


C4 


1x1 

Twisted 6 link 


3.248 
-0.2810 


-1.580 
0.0051 


0.1257 
0.0049 


0.0576 
-0.0096 



whether the saddle-point integration gives an action that can be approxi- 
mated by a few local terms. The application of this method to SU(3) and 

the tuning of the truncated version is still under investigation [99]. The 
current version of the truncated classical fixed-point action (CFP) is [100] 

ScFP = ^^{cl(l-ReTrH^*)-h...-h4(l-ReTrl^')'*}. (12.32) 

where the sum over i involves two kinds of loops, the plaquette and the 6- 
link loop that twists around a 3-cube {C^^ in Fig. 7). Higher powers of the 
loops arise naturally in the saddle-point integration. Since the strengths 
decrease rapidly, the authors choose to truncate at the fourth power. The 
choice of the particular 6-link loop is motivated by the desire to keep the 
action local in time and to help improve rotational invariancc. The values 
of coefficients ci — C4 for the two types of loops are given in Table 5. 

At this point full control over stability with respect to truncation has not 
been achieved. The second issue is whether this Ansdtz, obtained at a = 0, 
continues to be a useful parameterization at 1/a ~ 2 GeV where simula- 
tions are done, i.e. is the renormalized trajectory approximately linear. To 
test this requires a non-perturbativc tuning about 1/a ^ 2 GeV using the 
same 6 = 2 blocking transformation used in the saddle-point integration. 
The calculations necessary to do this tuning have not yet been initiated. 
As a result this Ansdtz is considered as work in progress. Nevertheless, 
the tests for improvement in the scaling of Tc, string tension, topological 
susceptibility, and in the restoration of rotational symmetry are already 
very encouraging [100,52]. 

12.8. Summary 

The question - which is the most improved gauge action - is still far from 
settled. As should be obvious, the drawback of the renormalization group 

improved actions is that the expansion is in the size of the loops or the 
dimensionality of the representation and not in the dimensionality of the 
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irrelevant operators. Since each loop contributes to operators of all dimen- 
sions it is not clear how to parameterize the residual error (O(a^) + 0{a'^), 
or O(a^), or higher). Also, implementing the equations of motions (remov- 
ing the freedom associated with redundant directions) has not been easy. 
So there is not yet a consensus on the obvious minimal set of operators that 
should be kept. Different operators have been added based on (i) their im- 
portance in perturbative expansions, (ii) size of the coefficients in blocked 
actions (keep in mind that the freedom of the choice of the blocking trans- 
formation and consequently the associated RT effects these coefficients), 
and (iii) specific incorporation of non-pcrturbative effects; for example the 
inclusion of a negative adjoint coupling helps avoid the phase structure in 
the fundamental-adjoint plane (see Section 15), and adding loops mimick- 
ing the F'^'^Ff^,^ operator help give scale invariant instantons [100,101]. One 
hopes that a more precise method for tuning the gauge action will emerge 
in the future. Fortunately, in spite of this ambiguity, the situation is very 
encouraging. The improvement seen in the restoration of rotational invari- 
ance (measured using qq potential as shown in Fig. 8), glueball spectrum, 
and finite temperature studies is remarkable. To a good approximation, 
most of these actions are very well represented by adding just the W^^^^'' 
term with a coefficient larger than the value suggested by tadpole improve- 
ment, and keeping a few higher representations of the plaquette. Such a 
truncated space keeps the computational problem tractable. 

The bottom line is that even though there is not a commonly accepted 
improved gauge action, the improvements seen with any of the above ac- 
tions arc significant. At present the more pressing challenge, in terms of 
the size of the errors and the computational cost, is improving the Dirac 
discretization. This is discussed next. 



13. 0{a) improvement in the Dirac Action 

The first attempt to improve the Wilson fermion action by removing the 
0{a) corrections was by Hamber and Wu [102]. They added a two link 
term to the Wilson action 



■>HW = 
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^Y.Mx)Mg'^My)- (13.1) 

It is easy to check by Taylor expansion that, at the classical level, the 0{a) 
artifact in Sw is cancelled by the HW term. One can include quantum 
corrections by calculating the relative coefficient in perturbation theory. 
It turns out that this action has not been pursued - in the 1980's due to 
the lack of computer power, and in the 1990's due to the merits of the 
Sheikholeslami-Wohlert (called SW or the clover) action [20]. 

13.1. The Sheikholeslami-Wohlert (clover) Action 

SW proposed adding a second dimension five operator to Sw, the magnetic 
moment term, to remove the 0{a) artifacts without spoiling Wilson's fix 
for doublers. The resulting SW action is 

Ssw = Sw ^ tp{x)a^i,Fi^^ip{x) . (13.2) 

The fact that this action is 0(a) improved for on-shell quantities cannot 
be deduced from a Taylor expansion. To show improvement requires an 
enumeration of all dimension 5 operators, the use of equations of motion 
to remove some of them, and incorporating the remaining two {m^ijjijj and 
mF'^" Fp_i,) by renornializing g'^{\ + bgma) and m —> m(l + bmma). 

The construction of the SW action has been explained in detail by Prof. 
Liischer at this school, so instead of repeating it, I direct you to his lectures. 
I will simply summarize results of tuning Csw using both perturbative and 
non-perturbative methods. A list of perturbatively improved values of Csw 
that have been used in simulations is as follows. 

Csw = 1 tree-level improvement 

Csw = 1 + 0.2659(1)5^ 1-loop improved [103] 

Csw = 1/^0 tadpole improved 

Csw = 1/^*0 (1 + 0.015951) 1-loop tadpole improved 

In the 1-loop TI result wq was chosen to be the fourth root of the plaque- 

ttc. As an exercise convince yourself that using the expectation value of the 
link in Landau gauge gives a value closer (1.6) to the non-perturbatively 
improved result (1.77 at 5 = 1) given below. For each of the above values of 
Csw, the errors are reduced from 0(a) (Wilson action) to at least O(asa). 

The recent exciting development, also reviewed by Liischer, is that the 
axial Ward identity Csw can be used to non-perturbatively tune Csw 
and remove all 0{a) artifacts provided one simultaneously improves the 
coupling g^, the quark mass mq, and the currents [104]. Estimates of 
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Csw have been made at a number of values of g for the quenched theory 
with Wilson's gauge action [104,105]. These data can be fit by [105] 



for 6/^2 > 5.7. 

The advantages of the SW action is that it is local and leaves perturba- 
tion theory tractable [106]. Adding the clover term is only a ^ 15% over- 
head on Wilson fermion simulations. An example of improvement obtained 
using the SW versus Wilson action for fixed quark mass (M^/Mp = 0.7) 
is that the deviations from the a — values in Mj^j^fa and Myector/ 
at a = 0.17 fermi are reduced from 30 — 40% 3 — 5% [105]. Further 
tests and evaluation of the improvement are being done by a number of col- 
laborations. I expect a number of quenched results to be reported at the 
LATTICE 98 conference, and these should appear on the e-print archive 
by late 1998. 

13.2. D234c action 

As mentioned before, there are two sources of discretization errors: (i) 
approximating derivatives by finite differences, and (ii) quantum correc- 
tions, for example those due to modes with momentum > tt/o that are 
absent from the lattice theory. Alford, Klassen, and Lepage have proposed 
improved gauge and Dirac actions that reduce (i) by classical improve- 
ment and (ii) by tadpole improvement [107,79]. Their current proposal for 
isotropic lattices, called the D234c action as it involves second, third, and 
fourth order derivatives, is 



Here Aj," is the n*^ order lattice covariant derivative (symmetric and of 
shortest extent). The term proportional to C3 is precisely what is needed 
to kill the leading discretization error in a|/'* as shown in Eq. 7.6. The 
terms proportional to r are generated by a field redefinition and thus rep- 
resent a redundant operator. They, therefore leave unalterated the 0{a) 
improvement of the naive action for on-shell quantities. The three terms 
in it are respectively the Wilson term, that is again used to remove the 



1 - 0.6084^2 - 0.2015,g4 + 0.03075^^ 
1 - 0.8743^2 



(13.3) 
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doublcrs, clover term (for F^^ they use an O(a^) improved version) as in 
the SW action needed to remove the 0{a) error introduced by the Wilson 
fix, and an correction. To include quantum corrections to the classical 
values C3 = Cf = Ci ~ I they resort to tadpole improvement, i.e., they 
divide all links appearing in M0234C by uq, which is taken to be the mean 
link in Landau gauge. 

Quantum corrections generally also induce additional operators, how- 
ever, in perturbation theory these start at 0{as). After tadpole improve- 
ment the coefficients of for such operators are all expected to be small, 
so one hopes that their contributions can be neglected. The philoso- 
phy/motivation is analogous to dropping the small c'^^ term in the tadpole 
improved Luscher-Weisz pme gauge action given in Eq. 12.14. 

For two of the coefficients, Cp and C3, the authors have carried out 
non-perturbative tests of the accuracy of tadpole improvement [79]. Since 
only the operator A|j'^' breaks rotational symmetry at O(a^), therefore C3 
can be tuned non-perturbatively by studying rotational invariance of dis- 
persion relations. Similary, the clover term and the Wilson term break 
chiral symmetry at 0(a), thus their relative coefficient Cp can be tuned by 
studying the axial Ward identity. Results of their tests show that errors in 
the spectrum and dispersion relations for momenta as high a,s p — 1.5 /a 
are « 10% even at o = 0.4 fermi! The authors are extending their anal- 
yses of the D234 action to anisotropic lattices which, as will be discussed 
in Section 18.4, have already yielded rich dividends in the calculation of 
glueball masses. 

13.3. NRQCD action 

In the simulation of heavy quarks one has to worry about O^mna) errors 
in addition to 0{pa) and 0{AQCDa). For charm and bottom quarks these 
errors are 0(1) in the range of lattice spacings one can reasonably simulate. 
Lepage and collaborators, recognizing that the heavy quark in heavy-heavy 
or heavy-light systems is, to a good approximation, non-relativistic, devel- 
oped a non-relativistic formulation (NRQCD) valid for niHa > 1 [108]. In 
this approach, the action for the two-component Pauli spinors is obtained 
by a Foldy-Wouthuysen transformation, and relativistic and discretization 
corrections are added as a systematic expansion in 1/M and p/M. One 
version of the NRQCD action is [109] 

r=?tV't (13.5) 

-^.(-^).('-^)>'('-^):.('-^)./'- 
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where tp is the two-component Pauli spinor and Hq is the nonrelativistic 
kinetic energy operator, 

A(2) 

6H includes relativistic and finite-lattice-spacing corrections, 



8(Mo)3 24Mo 16n(M( 



0) 



Here E and B arc the electric and magnetic fields, A and A^^^ are the 
spatial gauge-covariant lattice derivative and laplacian, and A^**) is a lattice 
version of the continuum operator ^Df. The parameter n is introduced 
to remove instabilities in the heavy quark propagator due to the highest 
momentum modes [108]. Again, the tree-level coefficients Ci = 1 can be 
tadpole- improved by rcscaling all the link matrices by Uq. 

This formulation has been used very successfully to calculate from 
the upsilon spectrum [110] as discussed in Section 19, the heavy-heavy and 
heavy- light spectrum [111], and the heavy- light decay constants /s and 
/bs [109]. The drawback of NRQCD is that it is an effective theory and 
one cannot take a 0. In fact, to simulate charm and bottom quarks, one 
has to work in a narrow window 1 < 1/a < 2.5 GeV to keep niHCi > 1. The 
good news is that current data [110,111,109] show that the discretization 
and relativistic corrections are small and under control with the terms 
included in Eq. 13.7. 



13.4- Perfect action 



The development of a perfect Dirac action is not complete. The method is 
same as used for the gauge action and discussed in [90]. The bottleneck 
is that the saddle point integration gives an action that even in the free 
field limit involves a number of sites, and the gauge connections between 
any pair of sites is the average of a large number of paths [112]. First tests 
of such an action truncated to a 3^ hypercube have been carried out by 
Degrand and collaborators [100]. From these it is not clear whether the 
improvement, compared to other discretization schemes, justifies the extra 
cost of simulating such an action. Hopefully better ways to tune and test 
such actions will be formulated soon. 



INTRODUCTION TO LATTICE QCD 



79 



13.5. Fermilab action 

The Fermilab group has proposed a "Heavy Wilson" fermion action [113]. 
This action is meant to interpolate between the light 0{a) improved SW- 
clover action and the NRQCD action for heavy quarks. The goal is to 
remove not only 0{a) discretization corrections but also take into account 
all powers of ma errors so that charmonium and bottomonium systems can 
be simulated with the same reliability as light-light physics. This improve- 
ment is especially targeted towards improving heavy-light matrix elements. 
So far no tests have been reported with the proposed action, therefore I 
will not reproduce the action here but refer interested readers to [113] for 
details. 



13.6. Summary 

In analogy to discussions of improved gauge actions, having seen tangi- 
ble evidence for improvement in quenched lattice data with each of the 
above formulations, the question is, which action is the most efhcient for a 
comprehensive study? The deciding factor, in my opinion, will be dynami- 
cal simulations since the prcfactor and the scaling exponent (L^^") makes 
the cost of update of background configurations prohibitive without im- 
provement on today's computers. Dynamical simulations using improved 
actions have just begun, so it will be a few years before enough data can 
be generated to determine the advantages and disadvantages of different 
approaches. 

Finally, I should mention a recent novel approach domain wall 

fermions for improving the chiral properties of Wilson fermions [114]. 

In this approach the system (domain) is 5-dimensional and the fermions live 
on 4-dimcnsional surfaces that make up the opposite boundaries (walls). 
The 4-dimensional gauge fields on the two walls are identical, whereas the 
fermion zero-modes have opposite chirality. The overlap between these 
modes is exponentially suppressed by the separation in the fifth dimen- 
sion. Consequently, the chiral symmetry breaking contributions are also 
supressed, and results with good chiral properties for vector like theories 
like QCD can be obtained by adding the contribution from the two domain 
walls. First simulations of weak matrix elements show results with good 
chiral properties and in general look extremely promising [115]. Since I do 
not do justice to this exciting approach in these lectures, interested readers 
should consult [115] and references therein. 
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14. Confinement and Asymptotic Freedom 

For any theory to provide a successful description of strong interactions 
it should simultaneously exhibit the phenomena of confinement at large 
distances and asymptotic freedom at short distances. Lattice calculations 
support the hypothesis that for non-abelian gauge theories the two domains 
are analytically connected, and confinement and asymptotic freedom co- 
exist. Similary, one way to show that QCD is the correct theory of strong 
interactions is that the coupling extracted at various scales (using experi- 
mental data or lattice simulations) is unique in the sense that its variation 
with scale is given by the renormalization group. The data for is re- 
viewed in Section 19. In this section I will discuss what these statements 
mean and imply. 

14.1. Asymptotic Freedom 

The term asymptotic freedom implies that the running coupling g ^ as 
the momentum scale of the probe — > 00. This is characterized by the 
/3-function 

m|^=-/3(5) = -(/3o5'+/3i5' + ---) (14.1) 

where P{g) is positive for non-abelian gauge groups. In this case 5 = is 
a UV stable fixed point of the theory, and about this point the theory can 
be analyzed using standard perturbative expansions. This novel behavior 
was confirmed by 't Hooft, Politzer, and by Gross and Wilczek who showed 
that the perturbative /3-function satisfies Eq. 14.1 [116], i.e., for N colors 
and nf active flavors 

ft = (^-i^-^^i^)/(I6,V. (14.2) 

These two leading terms in the expansion of l3{g) are gauge and regular- 
ization scheme invariant. From these it is easy to see that P{g) is positive 
for n/ < 8. This result was essential in establishing QCD as the the- 
ory of strong interactions for it explained existing experimental data which 
showed that the strength of strong interactions decreases as the momentum 
exchanged in a process increases. 

Asymptotic freedom, Eq. 14.1, implies that QCD dynamically generates 
a mass scale. Integrating Eq. 14.1 from momentum scale ni to /X2 with 
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^2 > Ail and keeping only the /3o term to simplify the standard calculation, 
gives 

2/305^(^2) ~ 2/3o52(/ii) " ' ^^^"^^ 

i.e. the coupling constant of non-abelian gauge theories depends logarith- 
mically on the momentum scale of the process. Rewriting Eq. 14.3 

— log/x = constant 



r 1 1 M 



introduces Kqcd, the invariant scale of the theory with dimensions of mass. 
Thus QCD, a theory with a dimensionless coupling constant and no intrin- 
sic mass scale in the absence of quark masses, dynamically generates a mass 
scale. This happens because in order to specify g one has to also specify 
a momentum scale at which it is defined. Extending the above analysis to 
include /3i in Eq. 14.1 gives 

^1 

Kqcd = ^im ^ (^^)'^ exp[- ^^^] ^ ^ fMl^)) •(14-5) 

This 2-loop definition of Kqcd is not unique; the value of Kqcd depends 
on the the precise relation between g and fx. However, once the value 
of A is determined in one scheme it can be related to that in any other 
pcrturbative scheme. For example, in the lattice regularized theory Kiatuce 
is also defined by Eq. 14.5 but with /i replaced by 1/a. Then to 1-loop 



Kqcd ( 1 



fj,a exp 



Klattice I 2/3o 



1 1 



g^in) g^{a) 



(14.6) 



In perturbation theory the two coupling constants are related as 

g\f,) = g\a)ll-Po9\a)(\og{i^af-\ogcA+0{g^)\ (14.7) 



By substituting Eq. 14.7 into Eq. 14.6 one finds 

Kqcd = C Klattice (14.8) 

i.e. the two constants, Kqcd and Kiattice, arc related by a multiplicative 
constant. To calculate C requires knowing the finite part of the coupling 
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Tabic 6 



The relation between Amom and ^.nd ^latt ^ ^ function of 

the number of active flavors. The results for ^mom / ^latt with g 
defined by the triple gluon vertex are in Feynman gauge. 








1 


2 


3 


4 




28.8 


34.0 


41.1 


51.0 


65.5 


^MOM /Matt 


83.4 


89.4 


96.7 


105.8 


117.4 



constant renormalization to 1-loop in both the lattice and continuum reg- 
ularization schemes [117,36]. The results are listed in Table 6 for Amom 
and A ^g, . It is important to note that even though the above definition 
and estimate of Aqc d is made using perturbation theory, it is intrinsically 
a non-perturbative quantity. 

All dimcnsionful quantities in lattice simulations arc measured in units 
of the lattice spacing, for example one measures Ma and not M. Such 
dimensionless quantities vary with g as 

Mitt = Cifi{g) = Cih^„„^_j,^^^a . (14.9) 

where the Cj are constants representing the continuum value and /j are 
scaling functions. Lattice correlation lengths = l/MiU diverge as the 
fixed point at (7 = is approached. This divergence is an artifact - it is the 
unit of measurement a that is going to zero and not Mj. This is precisely 
what one wants to have happen to get rid of the scaffold (lattice) erected to 
regularize the theory. Non-perturbative renormalization consists of taking 
the continuum limit holding some physical quantity Mj fixed while allowing 
a ^ according to ,fi{g). The particular quantity Mj (string tension, or 
nucleon mass, etc.) held fixed defines the renormalization condition, and 
the constraint between a and g is the scaling relation. 

— Scaling: The hypothesis of scaling is that all fi{g) converge to a uni- 
versal scaling function f{g) = Knon-pertd- The renormalization condition 
is then independent of the state used to define it. For g < g scaling each 
lattice quantity is characterized by the non-perturbative number Cj and all 
dimensionless ratios of physical quantities are independent of g. 

— Asymptotic Scaling: Close to the fixed point at 5 = 0, perturbation 
theory is apposite and all /j(.g) fpertig) Aqcoa- 

Lattice simulations, therefore, provide continuum results to within some 
error tolerance from simulations done in a neighborhood of = 0. In this 
region scaling holds. For a fixed error criteria, the extent of this scaling 
window can be enlarged by improving the action and operators as discussed 
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in Section 8. For simulations done outside the scaling window it is imper- 
ative that a reliable extrapolation to a = be done. Obviously, one way 
to get an upper bound on the errors is to compare the measured scaling 
function with /pert (.9)- The good news is that simulations in the past few 
years show that corrections to scaling can be reduced to a few percent level 
for lattice spacings with 5 < ^/a < 10 ( where ^/a = 1/^/aa). This corre- 
sponds to2< l/a<4 GcV. Consequently, realistic results with physical 
light quarks can be obtained on lattices of side ~ 75 — 150, and for many 
quantities, like those involving the heavy quarks, much smaller lattices are 
sufficient. The goal of the lattice approach is twofold. First, improve the 
discretization of the action and operators to extend the scaling region to 
stronger coupling, and second do simulations at a number of points within 
this region and remove the remaining errors by extrapolation to a = 0. 

An implication of the above discussion is that if all states in QCD have 
zero mass in the continuum limit then g = would be a trivial fixed point 
as the renormalizcd coupling gn would also be zero. On the other hand one 
can ask whether some of the states of QCD can be massless while the others 
are massive? It is interesting to note that we can actually say something 
about this aspect of the spectrum of QCD once the scaling behavior is fixed 
by Eq. 14.9. In fact, under the following two assumptions, the pure gauge 
sector of QCD will have a mass gap: 

— There arc no zero mass states at any finite non-zero value of g, i.e. the 
dimensionless quantity Ma does not vanish at any non-zero value of the 
lattice spacing o in a region that lies in the same thermodynamic phase as 
the continuum limit at a = 0. 

— There exists only one relevant coupling g in QCD, and corresponding to 
it a single universal scaling function defined by say Eq. 14.5. In this case 
the scaling behavior of all observables is fixed and all mass-ratios have to 
stay finite as the continuum limit is taken. 

These two assumptions arc borne out by present numerical data, i.e. there 
is no evidence for a massless state in the pure gauge sector. Also there is 
no indication of a phase boundary separating the region where simulations 
have been done and g = 0. Therefore, combining lattice data with renor- 
malization group scaling, QCD predicts that the lightest glueball state is 
massive. 

This prediction does not change with the introduction of n/ flavors of 
massless quarks into the theory. Spontaneously broken chiral symmetry 
will give rise to — 1 massless pions in the spectrum. To take the chiral 
limit one has to tune the quark masses to zero to define the physical world. 
This tuning has to be done at each value of a. Thus one cannot form mass- 
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ratios with respect to these Goldstone states as the corresponding Ci{g) are 
tuned to zero. 

To conclude, once we have verified that scaling (perturbative or non- 

pcrturbativc) exists then the only unknowns needed to predict the spec- 
trum of QCD are the constants Cj. Asymptotic freedom simplifies QCD 
by providing an analytical prediction of the universal scaling behavior as 
g ^ 0. Current data show that the coefficients q, which are intrinsically 
non-perturbative, can be extracted using LQCD from simulations done at 
a ~ 0.1 fermi. 

14-2. Confinement 

There are two ways one can test whether a theory confines. One, demon- 
strate that the free energy of an isolated charge is infinite; two, show that 
the potential energy between two charges grows with distance. The first 
test shows that it takes infinite energy to add an isolated test charge to the 
system, while the second shows that it requires infinite energy to separate 
two charges by infinite distance. I shall show that the first is probed by the 
expectation value of the Wilson line {£) in the pure gauge theory, while the 
second by the expectation value of Wilson loops (W) or by the correlation 
function (£(r)£t(o)). 

14-2.1. Wilson Loops 

The extra action associated with an external charge placed in a gauge field 
is »Sj = / d'^xJ^A'Jj^. For a point charge J{x) = 6'^{x), Sj is given by the 
path ordered integral of the gauge field along the world line of the charge. 
Thus 

(W) = e-^^ ~ e-^(^)^ (14.10) 

for a i? X T rectangular Wilson loop. One can regard the expectation value 
of a Wilson loop as the creation of a qq pair at time T = at point R/2, 
separated instantaneously to R and 0, allowed to evolve for time T, and 
then allowed to annihilate. Amongst the many terms (instantcous creation, 
annihilation, . . .) that contribute to the action for this procedure is that 
due to the potential between two static charges separated by distance R 
and integrated for time T, i.e. V{R)T. Thus 

V{R) = hm -llog(W(i?,T)) . (14.11) 

T — *oo 1 

The simplest Ansdtz for a confining potential is the Cornell potential 
V{R) = Vo + aR-^, (14.12) 
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where (T is the string tension. At large distances aR dominates, while 
at short distances it is the Coulomb term a/R that is relevant. Such a 
potential, therefore, simultaneously exhibits confinement and asymptotic 
freedom. For such linear potentials large Wilson loops show an area law, 
log(W(i?, T)) oc RT, and thus probe confinement. A state-of-the-art calcu- 
lation by the Wuppertal collaboration [118] of V{R) for SU(3) gauge theory 
is shown in Fig. 10. It shows the expected behavior - a linear rise at large 
R and a Coulomb behavior at small R. 

One can also extract the potential from the correlation function 
(£(i?)£^(0)) of Wilson lines, defined in Section 14.2.3, pointing in the T 
direction and separated by distance R. It is easy to see that this correlation 
function is equivalent to a Wilson loop with the sides at T = 0, Nt removed. 
The Wilson lines are closed due to the periodic boundary conditions on the 
gauge links and thus individually gauge invariant. This correlation function 
measures the action due to two static test charges separated by distance 



Measurements of the potential provide a way to set the scale of lattice 
calculations. One approach is to use a. The string tension can be extracted 
from the slope of V{R) at large R or from the asymptotic behavior of Creutz 
ratios [120] 



It can be measured very accurately in lattice simulations, however, its 

physical value is either inferred from the Regge slope or from potential 
models. Regge phenomenology gives the relation Mf{l) = l/a' + c where 
Ml is the mass of a state on the Regge trajectory with orbital angular 
momentum I. To get a relation between a and a' two approaches are used: 
(1) string models which give a = l/27ra', and (ii) potential models which 
give a = l/8a' [123]. The value of a' depends on the quark mass. For 
example a' from the light hadrons, cc, and hb is 0.85 0.85/2 -^k, 0.85/4.5 
GeV~^. It is customary to choose the value a' = 0.85 GeV~^ from the light 
sector as it is expected to be most sensitive to the large distance part of 
the potential. Combining this with the string Ansdtz a = Xjl'Ka' yields 
y/a ~ 440 MeV [124,125]. A survey of results from potential models can 
be found in [126] and one finds a varies in the range 0.18 — 0.22 GeV^. In 
most lattice calculations the value ^/g = 440 MeV is used, however, one 
should be aware of the significant uncertainty in its extraction. 



R. 




(14.13) 



86 Rajan Gupta 




12 3 4 

R/Rq 

Fig. 10. The static qq potential in the quenched approximation obtained by the Wup- 
pertal collaboration [118]. The data at /3 = 6.0, 6.2, 6.4 and 6.8 has been scaled by Rg 
defined in Eq. 14.14, and normalized such that V{Ro) = 0. The collapse of the different 
sets of data on to a single curve after the rescaling by Rg is evidence for scaling. The 
linear rise at large R implies confinement. 



An even better quantity, proposed by Sommer [127], is the distance scale 
i?o defined through the force between two static charges 



,2dV{R) 



= 1.65 

R=Ro 



(14.14) 



INTRODUCTION TO LATTICE QCD 



87 



o 



6 



4 - 



2 - 







T — I — I — I — I — I — I — I — I — I — I — I — I — r 

X /? = 5.6, /c=0.156 
> ig = 5.6, /c=0.1575 



T — I — r 




11^=2 potential 



J I I L 



J I I L 



J I I L 



J I I L 







2 

R/Rq 



Fig. 11. The static qq potential with two flavors of dynamical fermions obtained by the 
Wuppertal collaboration [119]. The dynamical quark mass in the f} = 6.0, k = 0.156 
simulation is ^ 2ms while that at k = 0.1575 is ~ rus- The rest is same as in Fig. 10. 
These data lie lower than the quenched points, however, for these dynamical quark 
masses there is no evidence yet of string breaking at large R. 



The Cornell and Richardson potentials, which have been very successful in 
predicting levels of charmonium and Upsilon system, yield 



i?o ~ 0.49 fernii 



(14.15) 



Setting the scale using Rq has the following advantages: (i) its value is 
not very sensitive to the choice of the phenomenological potential, (ii) it is 
defined at an intermediate distance where phenomenological potentials are 
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highly constrained by the cc and bb spectra, (iii) lattice measurements of 
V{R) at these intermediate distances can be obtained with high precision, 
and (iv) this construction is vahd for both full and quenched QCD. The 
status of the current quenched data from three collaborations is shown in 
Fig. 12. The data agree, and it turns out that over this same range of 
/3 the relation Rq = l.lS/i/a holds to a very good approximation. (The 
a is extracted from the large R behavior of the qq potential as discussed 
above.) Thus, in subsequent analyses of lattice data I will assume this 
relation when using either a or i?o to set the scale. 

For full QCD, the linear rising potential is screened. The string between 
the qq breaks when the energy is large enough to pop an additional qq out 
of the vacuum and create two mesons. This phenomena of screening can be 
tested by plotting V{R) for different values of the dynamical quark mass 
and observing the flattening of the linear rise. The value of R/Rq at which 
the flattening should occur should decrease with ruq. A recent analyses, 
again by the Wuppertal collaboration, is shown in Fig. 11 [119]. The data 
show no flattening - presumably because the sea quark masses used in the 
update are not light enough to cause string breaking at these R/Rq. 



14-2.2. Strong Coupling Expansion 

Conflnement, for /3 ^ 0, can be demonstrated using strong coupling ex- 
pansions. I will illustrate this technique by two examples - string tension 
and the 0++ glueball mass. 

The strong coupling expansion for SU(N) gauge theory depends on the 
following identities for integration over link matrices 



dg 



1 

















I dg Ui^j^ . . . Ui^jj^ — — eii...ijy ^ji...jN ■ (14.16) 

The essential point is that the group integration gives a non-zero result 
only if each link occurs in a combination from which a color singlet can 
be formed. Eq. 14.16 shows this for the identity, "meson" and "baryon" 
configuration of links. 

The expectation value of the Wilson loop, at small [3 (large g) for the 
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plaquette action, can be expanded as follows. 

(Wrt) = I y rff/W,iTe'^/^^(<^+^") 

= ^J dUWRT{l + P/2N{W+ + Wfi) + . . . ) (14.17) 

where W^^, W^i are the two orientations of the plaquette and the trace 
over color incidcs in each loop is implicit. An inspection of Fig. 13A shows 
that the first non-zero contribution to the integral occurs when the loop 
Wrt is tiled by elementary plaquettes with the right orientation. Each 
such plaquette brings a factor of (3/2N from the expansion and another 
factor of 1/N from the integration. (Note that SU(2) is special and the 
combined factor is p/N^ as the two orientations of the loop have the same 
value, i.e. the trace is real.) Thus 

(W«^) = f^(l + ---) (14-18) 

where the leading corrections come from replacing any given tile with a 
pillbox. Now using Eqs. 14.10 and 14.12 gives 

cT = -log^ + 0{P). (14.19) 

Thus all gauge theories, including U(l), confine in the strong coupling limit. 
What distinguishes U(l) from non-abelian gauge theories is that the U(l) 
gauge theory has a phase transition at ~ 1 (see Section 15), and the 
physical theory (electrodynaniic;s) lies in the weak coupling region which is 
analytically disconnected from the confining strong coupling limit. 

As a second application of the strong coupling expansion consider the 
calculation of 0++ glueball mass. The minimum tiling of the connected 
plaquette-plaquette correlation function is the rectangular tube consisting 
of four sides of size 1 x T as shown in Fig. 13B. Thus 



27V2 



(Wn(T)Wn(0))~e-^o++^=(:^) (l + ...), (14.20) 
and 

Mo++ = -41og^ + ... . (14.21) 

The fact that the lowest order result Mq++ / ^/a = 4 is very close to present 
day estimates (see Section 18.3) is fortuitous. 
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(A) Minimum tiling of a 6x6 Wilson loop. 
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(B) Tiling of one face of a plaquette-plaquette correlation function. 

Fig. 13. Examples of minimum tiling of (A) a 6 X 6 Wilson loop, and (B) the plaquette- 
plaquette correlation function. 



Exercise: What is the minimal tiling for the 0'^'^ correlation function built 
out of a Rx R Wilson loop. Show that Eq. 14-20 holds independent of the 
size of the loop. 

14.2.3. Wilson/Polyakov Line 

The Wilson line, pointing in say the time direction, is defined by the path 
ordered product 



Nt 



{C){x,y,z) = Tr P J| [74(0;, y, z, t) 



(14.22) 
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It is gauge invariant due to periodic boundary conditions on gauge fields, 
U4{x,y,z, Nt + 1) = U4{x,y,z,l). Repeating the above arguments made 
for the Wilson loop one finds that the free energy of an isolated quark is 
given by the Wilson/Polyakov line 

(£) ~ e-^"^^ . (14.23) 

Thus, the possibility of having isolated charges in the theory requires Fq 
be finite. 

In addition to gauge invariancc, SU(N) gauge theories have a global 
Z(N) invariance. It consists of multiplying all links on a given time slice 
by the same element z of Z(N) without changing the action. Under this 
transformation {€) z{C). Such a global symmetry can be spontaneously 
broken. Thus we have two possibilities 

Broken : {£.) ^ Deconfined {Fg finite) 
Unbroken : = Confined {Fg infinite) (14.24) 

Thus, (C) is an order parameter, and tests for confinement in pure gauge 
theories. It is also used to probe the location and order of the finite tem- 
perature transition - a discontinuity in {£.) signals a first order transition, 
whereas a continuous change implies second order. 

Svetitsky and Yaffe [128] argued that for the purposes of determining 
the order of the finite temperature transition the relevant effective theory is 
a 3-dimensional spin model of short-ranged interactions (between Wilson 
lines reduced to spin variables) with a Z(N) symmetry. Then using the 
universality classes known from renormalization group analyses, they pre- 
dicted that the transition for SU(2) should be second order (Ising class), 
whereas for SU(N> 3) it should be first order (Potts class). Numerical 
results are consistent with these predictions, in particular for SU(2) the 
critical exponents have been determined using standard techniques of Sta- 
tistical Mechanics and found to agree with those of the 3-dimensional Ising 
model [129]. 

In the presence of dynamical quarks, it is easy to see that the discretized 
Dirac action is not invariant under the Z(N) symmetry. Consequently, 
(£) ^ for all values of /?, so (jC) is no longer an order parameter. Never- 
theless, changes (discontinuous or continuous) in it are used to search for 
the location of the finite temperature phase transition and its order. 
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15. Phase Transitions in the Lattice Theory 

The lattice theory with a finite cut-ofF ii = n/a can be regarded as an 
effective theory. The integration over momenta in the range {7r/a, 00} 
renormalizes the couplings and generates additional effective interactions 
(sec Section 8 and the lectures by M. Liischer for a detailed discussion of 
this approach). Thus, one can regard the lattice theory as a point in an 
infinite dimensional space of couplings, and taking the continuum limit as 
a flow in this space to the critical point at a = that defines the physical 
theory. In this way of thinking about LQCD as a Statistical Mechanics 
system one is automatically lead to the questions 

— What is the phase diagram in this extended coupling constant space? 

— Are there other fixed points, and if so what is the nature of the theory 
at those points? 

— What arc the order parameters that characterize different phases? 

The starting point of my discussion of these questions will be the pure 
gauge theory. As before, the generalized form of the gauge action is 



where the sum over i in WJj^l is over all possible Wilson loops and j is over 
all representations of the loop. (3 is the overall coupling and the Cij give 
the relative strengths of the different loops. 

The three simplest gauge invariant probes of the phase diagram are 
(gauge non-invariant quantities have zero expectation values as shown by 
Elitzur [37]) 

— Wilson loops: These, as discussed in Section 14.2.1, contain information 
on the potential between static quarks. 

— Wilson/Polyakov lines. These measure the free energy of an isolated 
quark. 

— The chiral condensate {tpip) ■ This provides information on vacuum align- 
ment under spontaneous chiral symmetry breaking. 

LQCD, in the limit g 00, can be analyzed using strong coupling ex- 
pansions as discussed in Section 14.2.2. Using these it was shown that, in 
the limit g( — > 00, the expectation value for Wilson loops has an area law 
for all gauge groups (U(l) and SU(N)). Since electrodynamics (U(l)) does 
not confine, this region of coupling space cannot correspond to the physi- 
cal theory. In other words, there should exist a phase transition separating 
the strong coupling and weak coupling phases of non-confining theories like 
U(l). This has been established for U(l), both analytically [130] and by 




(15.1) 
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Monte Carlo simulations [131]. For a theory like QCD, with both eon- 
finement and asymptotic freedom, we need to know if the two regions of 
parameter space are analytically connected. Analytical methods like strong 
and weak coupling expansions can be used provided the range of validity of 
the expansions gives a sufficiently large overlap, otherwise non-perturbative 
methods are necessary. 

To address these questions for SU(3) I show, in Fig. 14, the behavior of 
square loops, {Re Tr W^^"^}, as a function of the coupling /? for the simplest 
action - Wilson's plaquette action. Also shown for comparison are the 
strong coupling expansion to 0(/3^) [133], the weak coupling expansion to 
0{g^) [134] and its tadpole improved (TIPT) version. These numerical 
data show that there is a smooth though sharp transition between the 
weak coupling and strong coupling phases, i.e. no phase transitions. The 
crossover takes place in the window 5 < /3 ~ 6, with the larger loops 
becoming perturbative at slightly higher 0. Both strong coupling and weak 
coupling expansions fail in this region (for TIPT this is evident only for 
large loops), so numerical methods are essential. Such an analysis was the 
basis of the pioneering work of Creutz, and Creutz, Jacobs, and Rebbi [132], 
who showed that the string tension changes from its strong to weak coupling 
behavior without a phase transition. Their calculations yielded the first 
clue that it may be possible to get continuum physics with 0(20%) errors 
from simulations at /? > 6.0. 

The TIPT result shown in Fig. 14 is obtained as follows. The mean-field 
improved Wilson loop expectation value is defined to be 

(ly^x') = (15.2) 

where for uq I use the non-perturbative value obtained from the plaquette 
as it is more readily available. The improved expansion for (PF*^'), after 
absorbing the tadpoles using the plaquette, is given by 

/W'X') = ^ (15.3) 

where both terms on the right are first expanded in terms of = 
g'^/uoipert) = *° order g^, and then the ratio is truncated 

at order g"*. It is obvious that using from the plaquette guarantees a 
perfect fit for (ly^^^), however the improvement in even the (ly^^**), is 
remarkable. 

The first study of the phase diagram in the generalized coupling constant 
space was done using an action constructed from the plaquette in both the 
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5 10 15 

^ = 6/g' 

Fig. 14. The behavior of the expectation value of square Wilson loops as a function 
of coupling /3. The data, show a sharp crossover between 5^/3^ 6.0. Result of the 
strong coupling expansion to 0(/3^) and weak coupling expansion to 0(g*) are shown 
by the dashed line. Results of tadpole improved weak coupling expansion (TIPT) are 
shown by the solid lines. Since the improved coupling is defined in terms of the 1x1 
loop, 3I = / < ly^^l >, the agreement between the data and TIPT is a trivial result 
of the construction as explained in the text. The improvement in larger loops, however, 
demonstrates the efficacy of TIPT. 



fundamental and adjoint representations of the SU(N) gauge theory [135] 
5 = /3F^^ReWixi + /3a^-^|WixiP . (15.4) 
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For this action, the effective bare couphng, obtained using a Taylor expan- 
sion, is 

2N 

/3e// = — =/?f + 2/3a. (15.5) 

Seff 

Thus simulations can be done in the negative /3a region as long as 
Pf > 2 1 /3a I- This is a classical condition that avoids the weak coupling 
singularity. 

The resulting lines of first order bulk transitions, along with the location 

of the end point for three gauge groups, are show in Fig. 15. The point 
at /3a = 00 is the first order transition in the Z{N) gauge theory, while 
that at /3ir = corresponds to that in the SO{N) theory. The location of 
the third end-point with respect to the f3p axis depends on N. Current 
estimates of the location are [136-138] 

(3f 1.22 /3a = 1.25 SU(2) 

Pf = 4.00(7) /3a = 2.06(8) SU(3) 

/3f-12-15 /3a = (-1) - (-5) SU(4) (15.6) 

The fact that for < 3 this point lies above the (3f axis is consistent with 
the observation that there is no discontinuity in the behavior of Wilson 
loop data as show in Fig. 14. If one were to take the continuum limit along 
a line crossing the transition line above this point, for example along the 
dotted line A in Fig. 15, then there would be discontinuities in Wilson loops 
(and consequently in the string tension), specific heat, and glueball masses 
at the point of intersection with the transition line. At the end-point X 
the specific heat diverges and consequently the 0++ glueball mass goes to 
zero. (Note that the specific heat is the volume integral of the 0++ glueball 
correlation function with the plaquette as the interpolating operator. Thus 
the divergence in the specific heat implies a zero in the 0++ glueball mass.) 
The string tension and masses of glueball states with other spin-parity 
quantum numbers are non-zero and continuous at this point. Such phase 
transitions are lattice artifacts. For example, taking the continuum limit 
at X, i.e. setting a = 0, would give a free theory as all states other than 
the 0"'""'" glueball become infinitely heavy. 

One can take the continuum limit for the SU{3) gauge theory along any 
trajectory like A or or B. In case A it is necessary to keep Pf > Pp lo 
avoid the lattice artifacts. Along W, while there are no singularities, the 
artifact X could, around Pf = 5.7, cause significant deviations from the 
physics of the relevant fixed point at g = 0. Lattice data verify this scenario 
- the non-perturbative /3-function measured from observables other than 
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Fig. 15. The phase diagram of SU(N) gauge theories in the fundamental-adjoint couphng 
constant space. 



Mq++ shows a large dip around jSp — 5.7 [89]. Thus, to avoid these artifacts 
requires simulations for j3F ^ 6.0. On a trajectory like B the influence of 
X is expected to be less than along W , and it is possible that contact 
with continuum physics can be made on coarser lattices. If so then the 
trajectory B corresponds to an "improved action" in accordance with the 
third criteria enumerated in Section 8. 

This qualitative argument for taking the continuum limit along a trajec- 
tory like B is supported by Monte Carlo renormalization group (MCRG) 
estimates of the renormalized trajectory [87]. Such calculations of the RT 
yield a negative value for /3^. As explained above, the most sensitive probe 
of the expected improvement on including a (5a term in the action, is the 
behavior of Mq++. Until recently tests of this conjecture were limited by 
statistical accuracy [166]. The recent calculations on anisotropic lattices of 
glueball spectra by Morningstar and Peardon [139] validate this scenario. 
I shall discuss their results in Section 18.3. 

It is clear that there are additional phase transitions in the generalized 
lattice theory. For example one would get a similar diagram to that in 
Fig. 15 if the action was defined in terms of the 1x2 loop in the funda- 
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mental and adjoint representations. The question, therefore, is whether 
there exist other second order transition points at which all correlation 
lengths (measured in lattice units) diverge i.e. where a non-trivial contin- 
uum limit exists. The commonly held belief, supported by numerical data, 
is that all such points lie on the hypersurface defined by geff = 0. The long 
distance behavior of these theories is controlled by the fixed point defined 
with respect to a suitable renormalization group transformation. Taking 
the continuum limit of the lattice theory on this hypersurface provides a 
non-perturbative definition of QCD. This is the standard scenario. 

An alternative scenario is the presence of a non-trivial fixed point at 
geff ^ [140]. If this unlikely possibility turns out to be true then the 
perturbative relation between g^ff and a, based on asymptotic freedom, 
would change - the coupling would no longer run to zero as the momentum 
scale Q ^ oo. This could be exposed by precise measurements of as a 
function of Q^, or by a detailed search for non-trivial second order phase 
transitions in the lattice theory. 

In either case the non-perturbative lattice analyses of quantities such as 
the spectrum and weak matrix elements does not change. To get to the 
continuum limit we would still extrapolate the lattice data in a to a = 0. 
Only the PQCD relation between a and g would be invalid, consequently 
extrapolations using the perturbative scaling relation could not be used. 
My bottom line on this issue, having mentioned the heretical view, is to 
proceed assuming that the standard scenario is correct. 

16. Errors in Lattice Results 

Numerical simulations of lattice QCD introduce both statistical and a num- 
ber of systematic errors. I have collected together a brief description of the 
various errors here, and in the remainder of the article I will only give 
quantitative estimates. I will label a particular systematic uncertainty 
as negligible/important if it is negligible/comparable in magnitude to the 
statistical error for that observable as determined in a state-of-the-art cal- 
culation. 

16.1. Statistical errors 

The monte carlo method for doing the functional integral employs statis- 
tical sampling. The results, therefore, have statistical errors. The current 

understanding, based on agreement of results from ensembles generated 
using different update algorithms, different initial starting configuration, 
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and different random number generators in the Markov proeess, is that the 
functional integral for QQCD is dominated by a single region. Second, we 
find that configurations with non-trivial topology are properly represented 
in an ensemble generated using a Markov chain based on small changes to 
link variables. Based on these observations we conclude that the energy 
landscape is simple, consequently, the statistical errors should fall asl/VN, 
where N is the number of independent measurements. Tests like binning 
of the data and the calculation of auto-correlation times for various ob- 
servables are used to determine the monte-carlo time between independent 
configurations. 

Estimates of decorrrelation times in dynamical simulations are just be- 
coming available and the understanding of decorrelation times for various 
update algorithms and different fcrmion formulations is not complete. A 
recent study of decorrelations for update of Wilson fermions [142] using the 
hybrid Monte Carlo algorithm [141] (the algorithm of choice) presented very 
encouraging results. They found that decorrelation time based on tunnel- 
ing rate between different topological sectors (expected to be amongst the 
slowest modes) was within a factor of two of standard probes like large 
Wilson loops and hadron correlators, and « 80 units of molecular dynam- 
ics time for m^/Mp > 0.56. Thus trajectory lengths of 5000-10000 units, 
which are feasible on todays computers, provide a reasonable statistical 
sample. Even though details like how this decorrelation time grows with 
decreasing quark mass still need to be filled in, it is clear that current 
Monte Carlo algorithms do provide a reliable way of doing the functional 
integral for both QQCD and full QCD. 

16.2. Finite Volume 

This is the simplest error to visualize as it is associated with using finite 
lattices to represent an infinite system. Liischer has shown that for suf- 
ficiently large lattices of side L, the finite volume corrections in the mass 
M of a given state fall off as e~^^ [27]. This anaylsis assumes that the 
interactions with the mirror states on a periodic lattice are small. On small 
lattices the wave-function of the states is squeezed, resulting in an much 
more rapid increase of the mass with decreasing lattice size. This effect is 
seen clearly for L < 1.5 fermi, where the finite size behavior fits a power- 
law. ^ [143]. The goal, therefore, is to work on lattices large enough 
such that the finite size effects are exponentially suppressed. 

To determine the lattice size needed for observables involving quark 
propagation I choose the pion as the test state as it the lightest and thus 
has the largest correlation length. Results of quenched simulations show 
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that for MttL > 5 the exponential suppression apphes and the corrections 
are neghgible. For physical value of this translates into L > 7 fermi. 
Another way of stating the same criteria is that the quantum mechanical 
properties of hadrons, typically of size < 1 fermi, arc unaltered if the box 
size is larger than 7 fermi. Current simulations explore quark masses in 
the range > ms/4 or equivalently M^/Mp > 0.4. For these "heavier 
masses" the criteria M^^L > 5 translates into L > 3 fermi. The lattices used 
in the most recent calculations, for example the CP-PACS runs, discussed 
later, satisfy this condition. 

Another consequence of finite L is momentum resolution. For example, 
the lowest non-zero momentum, and spacing, is 27r/La = 393 MeV for 
L = 32 and 1/a = 2 GeV. One would like a much finer resolution when 
investigating the dispersion relation for mesons and baryons or when cal- 
culating matrix elements required to study form-factors in semi-leptonic 
and rare decays. This can only be done by increasing either L or a. In- 
creasing L is limited by computer resources, while increasing a increases 
discretization errors which are discussed next. 

16.3. Finite lattice spacing errors 

The lattice is merely a technical scaffold and the physical results are ob- 
tained by taking the a ^ limit. At finite a, lattice results have dis- 
cretization errors whose size depend on the degree to which the lattice 
action and operators have been improved (see Sections 8, 12, and 13). As 
discussed in Section 8, we employ two strategics for removing these errors. 
First, improve the lattice action and operators so that the errors at fixed 
a are small, and secondly, repeat the simulations at a number of values 
of a and extrapolate to a = 0. Since the functional form used to char- 
acterize the errors is usually truncated to just the leading term in a, the 
extrapolation has an associated systematic error. To get a handle on this 
uncertainty, extrapolations have been done for three different formulations 
(Wilson, clover, staggered) for which the leading corrections are different 
(0(a), 0{aa) — 0{a^), 0(a^)). The difference between the three results in 
the a = limit is a measure of the residual uncertainty. I shall illustrate 
this procedure for controlling the discretization errors in Section 20 where 
the calculation of quark masses is discussed. 

16.4- Chiral Extrapolations in the light quark masses 

The physical u and d quark masses arc too light to simulate on current 
lattices. For 1/a = 2 GeV, realistic simulations require L/a> 70 to avoid 
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finite volume effects, i.e. keeping Mt^L > 5. Current biggest lattice sizes 
are L/a = 32 for quenched and L/a = 24 for unquenched. Thus, to get 
results for quantities involving light quarks, one typically extrapolates in 
m„ = nid from the range m,,/4 — 2ms using fimctional forms derived using 
chiral perturbation theory but truncated to the first couple of terms. Since 
the range of the chiral extrapolation is large, it is important to quantify the 
size of the higher order chiral corrections. This requires precise data at a 
large number of values of quark masses. Most of the current fits are based 
on keeping just the lowest order terms as the quality of the data is not good 
enough to quantify the various higher order corrections, and/or to identify 
the quenched artifacts discussed in Section 16.8. The status of the current 
data with respect to resolving these higher order corrections is discussed 
in Section 18 where an analysis of SU(3) flavor symmetry breaking, for 
example the mass differences between mesons and baryons within SU(3) 
multiplets like the baryon octet and the decuplet, is presented. 

Current simulations also neglect isospin breaking effects and electromag- 
netic contributions as these are a few MeV in nature, and smaller than the 
present numerical resolution. The iso-spin symmetric mass is defined by 
rn = (to„ + md)/2. For an exploratory analysis of iso-spin breaking effects 
in the presence of a strong U(l) field see Ref. [144]. 

I expect that our understanding of the chiral corrections will change 
significantly in the next two years due to the factor of 10-100 increase 
in the computational power made available to LQCD. To study iso-spin 
breaking effects properly requires dynamical lattices with roughly physical 
values for m„ and m^. Such simulations are still a few years away. 

16.5. Discretization of Heavy Quarks 

Simulations of heavy quarks (c and b) have large discretization errors of 
0{ma) in addition to the 0{Aqcd(i) and 0{pa). This is because quark 
masses measured in lattice units, rrica and mta, are of order unity for 
2 GeV < l/a < 5 GeV. Data show that these discretization errors are 
large even for rric in the case of Wilson and clover actions. (Staggered 
fermions do not have any advantage over Wilson-like discretizations for 
heavy quarks and are hence not used to study heavy quarks.) A num- 
ber of alternate approaches are being investigated to address this issue 
(sec Section 13). These include the non-relativistic QCD (NRQCD) for- 
mulation [108], lattice versions of the heavy quark effective theory [145], a 
variant of the Dirac discretization that interpolates between NRQCD and 
clover discretizations [113], and non-pcrturbatively 0{a) improved SW- 
clover fermions [104]. Many collaborations are testing these formulations, 
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however, it is still too early to judge which approach will provide the best 
solution to the simulations of heavy quarks. 

16.6. Matching between lattice and continuum scheme (renormalization 
constants) 

Experimental data are analyzed using some continuum renormalization 
scheme like MS so, to make contact with phenomenology, results in the lat- 
tice regularization scheme have to be converted to the continuum scheme. 
The pcrturbativc relation between renormalized quantities in MS and the 
lattice scheme are, in almost all cases, known only to 1-loop. Data show 
that the 0{as) corrections can be large, ~ 10 — 50% depending on the 
quantity at hand. However, as explained by Parisi [77] and by Lepage 
and Mackenzie [78], these large corrections are mostly due to lattice ar- 
tifacts which, once identified as coming from tadpole diagrams, can be 
removed. The Lepage-Mackenzie prescription for reorganizing the lattice 
perturbation theory to remove these artifacts has significantly reduced the 
associated uncertainty. Nevertheless, the improvement is not uniform and 
in many cases the corrections arc still large, and the residual errors in TI 1- 
loop perturbative estimates are hard to ascertain. The final step in improv- 
ing the reliability of the matching factors is to develop non-perturbative 
methods [146-148]. Significant progress has been made in setting up these 
calculations, and once they are complete the uncertainties in existing re- 
sults due to using 1-loop Z's will be removed. An illustration of the effect of 
using perturbative versus non-perturbative Z's is presented in Section 20.1. 

16.7. Operator mixing 

As mentioned in Section 3, a major application of LQCD is to calculate 
the matrix elements of operators that arise in the effective weak Hamilto- 
nian between hadronic states. In addition to the question of the overall 
normalization of these operators, the operators can, in general, mix with 
operators of the same, higher, and lower dimensions. On the lattice this 
mixing arises due to quantum corrections and discretization errors. The 
set of lattice operators one has to consider is usually larger than those in 
the continuum theory because at finite a the symmetries of the lattice the- 
ory are smaller than those of the continuum theory, for example the hard 
breaking of chiral symmetry in Wilson fermions. In cases where there is 
mixing with lower dimensional operators, the mixing coefficients have to 
be known very accurately otherwise the power divergences overwhelm the 
signal as a ^ 0. Similarly, in cases where there is mixing with operators of 
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the same dimension but with different tensor structures, as in Wilson-hke 
actions due to the exphcit breaking of chiral symmetry, the chiral behavior 
may again be completely overwhelmed by these lattice artifacts if the co- 
efficients arc not known precisely. Examples where such mixing has posed 
serious computational challenges are the matrix elements of 4-fermion op- 
erators needed in the calculation of the A/ = 1/2 rule, Bk, Bq, and Bg. In 
these cases also non-pcrturbative methods for calculating the mixing coeffi- 
cients are essential. For a discussion of these methods see Prof. Martinelli's 
lectures. 

16.8. Quenched Approximation (QQCD) 

This approximation is computationally the hardest to shed and quantify. 
Since it plays a key role in today's simulations, I will discuss it in some 
detail. The approximation consists of ignoring the fermion contribution to 
the path integral, i.e. setting det M = constant in Eq. 4.3. The sole reason 
for this approximation is limitations of computer power. With known algo- 
rithms the computational requirements of full QCD simulations go up by 
factors of lO'' — 10^ depending on the quark mass. Since QQCD is confin- 
ing, asymptotically free, and shows spontaneous chiral symmetry breaking 



of the; l)ackgTOund gauge configurations, the physics analyses arc identical 
in most cases. It is therefore considered reasonable to do detailed quenched 
simulations in order to understand and control all other sources of errors 
while waiting for better algorithms or the computer technology needed to 
provide the required additional factor of 10"^ — 10^, i.e. 10-1000 teraflop 
capability. 

Physically, the quenched approximation corresponds to turning off vac- 
uum polarization effects of quark loops. This is illustrated in Fig. 16 for 
the pion correlator. One important consequence of neglecting vacuum po- 
larization effects is the behavior of the potential between a qq pair as a 
function of the separation. In the full theory the string breaks by the cre- 
ation of a qq pair out of the vacuum, therefore at large distances there is a 
screened potential between two mesons. For such screened potential large 
Wilson loops should show a perimeter law. In QQCD the string does not 
break and large Wilson loops show an area law. Thus the long distance 
behavior of the two theories is very different and one might be led to believe 
that QQCD is a bad approximation. However, because of confinement, the 
long distance scale that is relevant to hadronic physics has a natural cutoff 
of a few fermi. Thus, if one can match the potential between some fine 
scale (~ 0.01 fermi, below which asymptotic freedom makes the calcula- 




^ 0), and differs from full QCD only in the relative weighting 
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(A) Quenched QCD: quark loops neglected 




(B) Full QCD 

Fig. 16. An illustration of the difference between quenched QCD and QCD. QQCD has 
no vacuum polarization loops on virtual gluon lines. 

tions perturbative) and the typical size of hadrons 1 fermi) then QQCD 
could give reasonable estimates for the hadronic spectrum. This argument 
is really only applicable to heavy onia, where potential models work well. 
For light quarks one simply has to do the simulations to estimate the size 
of the distortions. In any case the community has proceeded by assuming 
that the two theories can be roughly matched by adjusting the quenched 
coupling constant at some scale like 0.5 fermi, or equivalently by adjust- 
ing the overall scale Aqqcd , and hopes that the quenching corrections are 
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small, 10 — 20%, for many cases. If this bears out then even QQCD would 
yield many phenomenologically useful predictions. 

In spite of the optimism, the fact remains that QQCD is not a uni- 
tary theory, so it is important to understand how it fails. I would like to 
illustrate its shortcomings with the following examples. 

We can calculate the piriT coupling in both QCD and QQCD by measur- 
ing the three point function (T[p7r7r]) illustrated in Fig 17A. The difference 
between the two values is a measure of the validity of QQCD. On the other 
hand the diagram shown in Fig 17B is absent in QQCD, and the p spectral 
function does not have a cut beginning at 2M^. Thus, the ptttt coupling 
cannot be extracted from an analysis of the p 2-point function in QQCD. 

A very important difference between full QCD and QQCD is the behav- 
ior of the 77'. In full QCD the singlet axial current is anomalous and the 77' 
acquires a large mass due to the summation of vacuum polarization graphs 
shown in Fig 18, 



p^-M^ p^-M^ ^p-'-M^ p'^-M^-M^ ^ ' 

where, for simplicity I have approximated the gluonic coupling by a 
momentum independent constant Mg. This is related to the topologi- 
cal susceptibility via the Witten-Veneziano relation Mq = 2n/Xt//^ = 
Mfy + - 2M|. [149]. In QQCD the absence of the vacimm polarization 
diagrams means that the rj' propagator has a single and a double pole, i.e. 
(p2 - M2)-i and the "hairpin" term {p^ - M'^)-'^M'^{p^ - M'^)-'^ , where 
is the Goldstonc pion mass M"^. The fact that the 77', in the limit 
mq ^ of QQCD, is massless has important consequences. Two groups, 
Sharpe and collaborators [150], and Bernard and Golterman [151], have 
investigated these by formulating a chiral perturbation theory for QQCD 
which retains rj' as an additional Goldstone boson. I illustrate the differ- 
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Fig. 18. The iteration of the r]' propagator in full QCD. Only the first two diagrams 
survive in QQCD. 
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(A) (B) 
Fig. 19. The 1-loop correction to the pion propagator in (A) full QCD (B) QQCD. 



ences between xPT and QxPT by considering the behavior of the 2-point 
correlation function of the pion. 

At 1-loop full QCD has the diagram shown in Fig 19A which is absent in 
QQCD. On the other hand, in QQCD, the hairpin term gives a contribution 
via the diagram shown in Fig 19B. Such diagrams are neglected in full QCD 
as they iterate to produce a massive 77'. The chiral expansions for M,r and 
f„ in full QCD have been calculated by Gasser and Leutwyler [152], 

1 16 

= 2mB[l + X„- -Xr,, + ^{2Ml + Ml){2Le - L4) 

1 fi 

+ Ml{2Ls-L,) + ...] 
Jo 

fl = /o [1 - -2Xk + ^ {2Ml + Ml)U + ^M^L, + . . .] (16.2) 

Jo Jo 
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where Li are theO(p^) coefficients, Xi = Mf log(-rP-)/(47r/o)2 are the chi- 

XSB 

ral logs, and = 131 McV. In QQCD, Sharpe and Bernard and Golterman 
[150,151] find that the leading terms are 

M^Q = 2mBQ[l-61og-j^ + ...] 

1 fi 

/.^ = /q[i-7tM2l5 + ...] (16.3) 

JQ 

where S = Mo/247r^/^ and Li are the quenched 0{p'^) coefficients. The 
term proportional to 5 is a quenched chiral log and is singular in the limit 

TOq = 0. 

The differences, in general, between quenched and unquenched expres- 
sions, as illustrated by the above expressions, are 

• All the chiral constants are different in the two theories. For example, 
B^Bq and /o ^ fQ. 

• The normal chiral logs are missing in QQCD. 

• The quenched expressions can be singular in the chiral limit due to 
the goldstone r/', i.e. through terms proportional to 6. These artifacts 
become dominant as niq — > 0. 

The first question raised by the analyses of Sharpe and Bernard and 
Golterman is - are the results of quenched x^T, an expansion about niq = 
0, meaningful if they are singular in that limit? We do not have a formal 
answer to this question. The approach adopted is to do the simulations and 
look for such chiral logs by doing fits with and without these terms. The 
status of this search has been reviewed in [153,154]. The data, while not 
conclusive, docs suggest that S w 0.15, consistent with the value obtained 
by using phenomenological values for Mq and fj^. 

The second question, assuming that predictions of quenched make 
sense, is how docs one extract phenomenologically useful results for quan- 
tities that require a chiral extrapolation to m? One approach is to fit the 
quenched data using the full QCD chiral expansions in a region of quark 
masses where the quenched artifacts arc small. Current data suggest that 
this region is niq > ms/4:. In this case quenching errors are, by definition, 
a combination of those due to the absence of quark loops and those due 
to fits being made at heavier quark masses. The second possibility is to 
fit using the QQCD expression but extrapolate keeping only the terms in 
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the full QCD expressions. This approach has the disadvantage of increas- 
ing the number of free parameters. In either case the hope is that there 
exist quantities for which the QQCD artifacts are small and the differences 
between QQCD and QCD coefficients arc small; in such cases quenched 
estimates are useful phenomenological predictions. One such example is 
the K'^K^ mixing parameter Bk [155]. 

In view of these various sources of systematic error, any quantity mea- 
sured on the lattice depends not only on the input parameters, i.e. the 
quark masses and the coupling Ug, but also on the lattice size L, a (dis- 
cretization errors are oc g"^a"^ with the powers depending on the order of 
improvement in the action and operators), the method used to determine 
the renormalization constants, and the number of dynamical quark flavors 
Hf. While the theoretical understanding of these errors is, in general, not 
complete, nevertheless, a lot is known and one can make consistency checks 
on the data. For example, the infinite volume continuum limit results for 
fixed nf simulations should be independent of the discretization scheme, 
the numerical approach used, and the definition of the renormalization 
constants. As of this writing, results for a number of observables show 
stability once extrapolated to the L ^ oo and a ^ limits, albeit in the 
quenched approximation (n/ = 0). We regard this consistency check as 
the first important step towards providing precise results. Our optimism 
stems from knowing that for many phenomenologically important quanti- 
ties {Bk, fo, Jd^ , Jb, ,Bb,C(s, and the quark masses) quenched lattice 
QCD results are already competitive with the best non-lattice methods, 
and future enhancements in methodology and computer technology will 
systematically improve their precision. 

17. Lattice Correlators 

In Section 4 we discussed the reduction of the pion 2-point correlation func- 
tion to a product of quark propagators and possibly gauge links, and its 
relation to sum over states from which and can be extracted. Here 
I enlarge the discussion to general 2-point and 3-point functions. Con- 
sider the calculation of the matrix element (^K'^\Vfj,\D^') which arises in the 
extraction of semi-leptonic form factors. The 3-point correlation function 
needed to calculate this ME is 




(17.1) 
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= E e"'^'-^+'""M^(0)75^(0)c(2/)7ps(2/)^(2^)75c(x)) 
= E (75^1^(2;, 0)75^1^(0, vh^S^Av, ^)) 

x.y 

^ exp{-EK{p + q)ty - ED{p){t^ ~ty)) 
(4EK{p + q)ED{p)) 
(0|m75s|X(p + q)) {K{p + q)\-s^,,c\D{p)) {D{p)\c-f5u\Q) . 

I have used (i) Wick contractions to express the correlation function in 
terms of quark propagators Sp (see Fig. 20); (ii) local bihnears for the 
source and sink operators and the vector current (to improve the signal 
one usually uses smeared operators and an improved current); and (3) the 
interpolating operator for the D meson and the current insertion are at 
definite momentum. Momentum conservation then fixes pk- The final 
form in Eq. 17.2 shows the behavior for large Euclidean time separations 
{tx <^ ty <^ to). Of the many terms in the right hand side, we are only 
interested in the matrix element {K{p + q)\cjf^s\D{p)) . The remaining 
factors can all be extracted from 2-point correlation functions. Thus, one 
can make individual fits to all the required correlators, or design a ratio 
in which the maximum number of unwanted factors cancel. In practice, to 
improve the numerical signal, a single fit is made to a ratio of 3-point to 
a combination of 2-point functions, {T[KVf,D])/{T[KK]){T[DD]), that 
gets rid of the exponential in time behavior. 

To make fits to such correlation functions projected on to definite spatial 
momenta, we have to address three questions: 

— Is the correlator even or odd in its time variables? 

— Is the correlator even or odd in its momenta? 

— Is the correlator real or imaginary? As mentioned before only the real 
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Table 7 



The signs for bilinear currents under T, V and CVH transformations. 





1 


7i 


74 




7»74 


7i75 


7475 


75 


TA 


+ 


+ 




+ 






+ 






+ 




+ 


+ 




+ 






XA 


+ 




+ 


+ 




+ 







part of Euclidean correlation functions has a signal. In practice we usually 
work in terms of bilinear operators given in Table 8 where factors of i have 
been dropped for brevity. In this basis of operators one needs to know 
whether a given n-point function has a signal in the real or imaginary part. 
To answer these questions we use the discrete symmetries, T, V, and CVH 
discussed in Table 2 and in Eq. 10.11. For each element of the Dirac algebra, 
define three signs, ta, t^a, and xa, associated with T, V, and CVH: 

C7475rA7574C-' = X^r^. (17.2) 

These signs are given in Table 7 and the application to 2-point and 3-point 
correlators is illustrated in the next two sub-sections. 

17.1. Two-point meson correlators 

Consider the general two-point correlator 

ip, t) = e-^^"^ (?2 ix)^B ^1 (x) (0)r^^2 (0)) , 

= -Ea.e-'P"{TT {S2{0,x)TbSi{x,0)Ta)}, (17.3) 

where Tr is trace over spin and color indices and t = X4. T, V and CVH 
then imply the following: 

- If tatb = ±1, C^'^(p, t) is even (odd) in t. 

- If ttattb = ±1, C^^{p,t) is even (odd) in p. 

- HxaXb = ±1, {p,t) is real (imaginary). 

Note that in Eq. 17.3 the second quark propagator is from x — > 0. Using 
the hermiticity property S{x, 0) = 75S'(0, a;)^75 we can write this in terms 
of S{x, 0). This simple property, the relation between the quark propagator 
from a given point and the anti-quark propagator to that point, leads to 
a huge computational saving. For example, if the momentum projection is 
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Table 8 

The local interpolating field operators for mesons and baryons in Wilson-like theories. 
Projection to zero momentum states is obtained by summing over the points x on 
a time slice. The C-parity is only relevant for flavor degenerate meson states. The 
in baryon operators denote spin trace. C = 7274 and the symmetry properties of 
flavor indices for nucleons are discussed in the text. The decuplet baryon operator is 
completely symmetric in flavor index 



State 




Operator 


Scalar(<T) 


i-(o++) 


u{x)d{x) 




l-{0++) 


u{x)'y4,d{x) 


Pseudoscalar 


i-(o-+) 


u{x)'y5d(x) 




i-(o-+) 


u{x)'Y4j5d{x) 


Vector 


1+(1— ) 


u{x)^id{x) 




1+(1--) 


u{x)-yij4d{x) 


Axial (ai) 


1-(1++) 


u{x)'yiiBd{x) 


Tensor(6i) 


1+(1+-) 


u{x)'ri'yjd{x) 


Nucleon octet 


2V2 / 


{ulCdbhaSce"^'' 




-(-~) 

2^2 


(«rC754)5ce"'"= 


Delta decuplet 




(«JC7id6)Sc€<'''= 



done by summing over the sink points x with weight exp(ip • x), then, 
because of the hermiticity property we need only one propagator inversion 
instead of + 1. 

17.2. Interpolating operators 

For Wilson-Hke fermions the interpolating field operators for mesons and 
baryons at zero-momentum are given in Table 8. The three flavors are 
labeled u, d, s for brevity. The baryon operators need further clarification 
as discussed in [156] and reproduced below. 

The spin-1/2 baryons are created by the interpolating operators 

Oirj)k = «.C75^6,,)^c,/ce'^''= , (17.4) 

where a, b and c label color, while i, j and k label flavor. It is simple 
to show that 0(^ij^k = ~^{ji)ki so that there are only nine independent 
operators- eight 5t/(3) octets and the singlet 'Y^^jk^^^'^^iok- One way to 
project against the singlet is to form 



^ijk — ^(ij)k + ^(ik)j — ^ikj ■ 



(17.5) 
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<du>s 



(dus) 



Fig. 21. The two different types of contractions for the baryon states. 



There are eight independent Bijh's, the relation to the usual states being 
exemplified by 

— ^duu '^^uud '^^udu 1 



S — Buds ~t~ Bdsu — Bsud • 
\/3A° = Buds - Bdsu ■ 



(17.6) 



The overall factor in these equations is arbitrary, while the relative nor- 
malization is fixed by SU (3) symmetry. 

All spin- 1/2 baryon correlators are built out of the two contractions 
shown in Fig. 21. The notation {DU)S = {UD)S corresponds to quarks 
of flavors U and D contracted into a closed loop, while the propagator for 
S carries the spin quantum numbers of the baryon. The notation (DU S) 
corresponds to a single ordered contraction of the three quarks. We consider 
two types of correlator, "S-like" and "A-like". The former is exemplified 
by that of the S° 



S{UD} = S{DU} : 



(Bsudix) BsudiO)) 

{US)D + {DS)U+ (USD) + (DSU) . 



(17.7) 



This equation defines the sign conventions for the contractions. The proton. 
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neutron, S+, E , and ^ correlators are also of this type: they are, 
respectively, D{UU}, U{DD}, S{UU}, S{DD}, U{SS} and D{SS}. The 
second type of correlator is that of the A° 



S[UD] = S[DU] = i 



{US)D + {DS)U + 4{UD)S 
-{USD) - {DSU) + 2{SUD) 
+2{SDU) + 2{UDS) + 2{DUS) 



(17.8) 



When rn„ ^ rud, there is also a non- vanishing A° — S° cross correlator, for 
which no useful results have been found [156]. 

Correlators of the form and A{AB} are not independent — they 

are related by 

A[AB] = \ (2,B\AA\ + B{AA}) , 

A{AB} = i {B[AA] + 3B{AA}) . (17.9) 

One can think of the results for the A{BB} and ^[SB] masses as being 

those for the E and A, respectively, with m,, = mA and m„ = rud = tub- 
Unlike in the real world, there is nothing to stop niA = rns- Note, however, 
that in this case the E and A are also degenerate, i.e. M{A{AA}) = 
M{A[AA]). Indeed, the contractions in the two cases arc identical. 

The interpretation of the results for the completely non-degenerate cor- 
relators, j4[BC] and A{BC}, is more complicated. Because isospin is bro- 
ken, the E^- and A-like states mix, with both correlators containing con- 
tributions from both physical states. Let M+ and M_ be the masses of 
the heavier and lighter states, respectively, and SM the mass difference. 
At long times, the effective mass for both correlators will asymptote to 
M_. However, at times short compared to the inverse mass difference, i.e. 
^imax 1, there will be an approximate plateau at a value which is 
a weighted average of the two masses. The surprising result, derived at 
lowest order in the chiral expansion in [156], is that the masses extracted 
from these short time plateau are insensitive to the isospin breaking and 
can be interpreted as those for A[DD] and A{DD} where D = {B + C)/2 
is the mean mass for the two quarks. 

17.3. Three-point correlators 

A general three-point correlation function of meson operators and bilinear 

currents has the form 

C^''^{q,ty;p,t,) 
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= _^e-*(«-^+f-)(Tr {Ss{0,y)rcS2iy,x)rBS,{x,0)TA)), (17.10) 
where tx = X4 and ty = yi. T,V and CVH imply 

— If tatbtc = ±1, C^^'^{q, ty;p, tx) is even (odd) under the simultaneous 
inversion of tx and ty. 

— If ttattbttc = ±1, C^^^{q, ty;p, tx) is even (odd) under the simultaneous 
inversion of p and q. 

— If XaXbXc = ±1, C'^^^{q,ty;p,tx) is real (imaginary). 
174. CVS Symmetry 

The theory has an additional discrete symmetry, CPS, where S is the 

symmetry between switching s and d quarks [157]. This symmetry holds 
both on the lattice and in the continuum for = vris- For ma 7^ Ws the 
symmetry is softly broken, i.e., violations are proportional to nig — ma. 
This symmetry plays a very useful role in the calculation of weak matrix 
elements for it (i) restricts the set of possible Wick contractions, and (ii) 
restricts the set of operators that mix with each other. 

For a simple illustration consider (s7^iu)(u7^Ld), a local 4-quark op- 
erator which is even under CVS 
■p 

{sjij,Lu){u'yij,Ld) (s747^L74M)(u747p£74d) 

(s^C" ^747^L74Cm^ ) (uC " ^ 747^ L^^Cd^ ) 
= (s^747j-'^74u'^)(w'^747j^74rf^) 
= (a^ji?^tl^)(n7ji?^rf^) 
= {uR'Y^s)(dRj^u) 
= {u'yi^Ls)(d'yij,Lu) 

^ {u^^Ld){sj^Lu) (17.11) 

The brackets denote trace over spin and color and the final reordering of the 
two traced terms is allowed as time-ordering is implicit when calculating 
expectation values. 

This symmetry has a further generalization as discussed by Bernard and 

Soni in [2,6]. For 4-fermion operators of the form {iJjiTAijj2){ip3^Bip4), there 
are two general switchings: (i) CVS' corresponding to 1 ^ 2 and 3^4, 
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and (ii) CVS'' corresponding to 1 ^ 4 and 2^3. Bernard and Soni show 
in [2,6] how these switching symmetries are used to restrict (i) the possible 
contractions in the decay K ^ mr via the operator {s'y^Lu){u'y^Ld) , and 
(ii) the mixing of 0± = {s"f^Ld){ujfj_Lu) ± {s"/f^Lu){u"/f^Ld) — (u c) with 
other operators. I leave the details of these calculations as an exercise for 
the reader. 



18. Lattice Calculations of the Hadron Spectrum 

Successful calculations of the spectrum will test QCD. The basic steps in 
the lattice analysis are discussed in Sections 4 and 17. For a given set of 

input parameters, the mass of a given state is determined from the rate of 
exponential fall-off of the connected 2-point correlation function 

r(.) . iO,ir)Om)-{Of) m = E^^^^|lg^e---.(18.1) 

n 

We shall assume that the sum over n is restricted to the desired state by 
(i) optimizing Oi and Oj to get a large overlap with the wave function, 
i.e. make {^desired state\O\0) large and the rest small, (ii) increasing the 
statistics so that the signal extends to large enough r at which any re- 
maining contamination from higher states is negligible; and (3) making a 
zero-momentum projection on either O; or Oj. I shall further assume that 
the spectrum of possible states is discrete and has a mass-gap. Then, in 
order to extract M it is useful to define an effective mass 

Me//(r) = log(^^^^) (18.2) 

which, in the limit r ^ oo, converges to the desired value. The anatomy 

of the behavior of Mg/ / (t) for two different states (pion and nucleon) and 
for two different choices of Oj and Of are shown in Figs. 22 and 23. The 
points to note are: 

— The convergence to the asymptotic value M can be from above or below 
depending on the choice of O. Only for O/ = 0| is the correlation function 
positive definite and the convergence is monotonic and from above. 

- For Of ^ 0|, the convergence depends on subtle cancellations between 
contributions of various states. The relative signs are given by the product 
of the two amplitudes (n|O/|0)(n|Oi|0). Since this sign for state |n) can 
be plus or minus, the convergence is neither guaranteed to be monotonic 
nor from above. 
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Fig. 22. Effective mass plots for M,r versus r. The horizontal lines denote the estimate 
of the asymptotic value of M and the range of the plateau used in the fits, while the 
dotted lines show the error. The squares denote data for smeared-smeared interpolating 
operators {Of = 0\), while fancy crosses show data with wall-local operators. 



— The contribution of excited states is most obvious at small r, where 
Me// varies rapidly before flattening out. For r in this flat (called plateau) 
region one assumes that, within the numerical precision, only the lowest 
state contribution is significant. 

— The location in r of the onset and the length of the plateau region de- 
pends on the interpolating operators. 

— For a finite statistical sample the plateau region is not flat but the data 
usually exhibit correlated fluctuations [158]. Therefore, unless the plateau 
is long enough that the fit can average over a few cycles, the extracted 
Me// can be biased by these statistical fluctuations. 
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Fig. 23. Effective mass plots for Afjv versus t. The rest is same as in Fig. 22. 



— The statistical errors grow with r, except for the case of the pion. The 
reason for this was elucidated by Lepage in [159] and a synopsis of his 
argument is as follows. For a state with N valence quark lines, where 
N — 2(3) for mesons (baryons), the errors are controlled by the square of 
the correlator which has 2N lines. The state of lowest energy in the squared 
correlator is, therefore, two pions for mesons and three pions for baryons. 
Thus, while the signal decreases as exp (—Mr), the noise only decreases as 
exp(-M.;rT X NI2). Consequently, all states for which M > the 
errors will increase with r. The only exception to this condition is the pion 
for which the signal and noise have the same t dependence. 

— In cases where there is no clear plateau (the signal does not extend 
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Fig. 24. Effective mass plots for Mp versus r for the wall-local source. Since there is 
no clear plateau, a fit including two states is made to 3 < r < 26. From such fits the 
masses of p and p* are obtained. The mass of the p from the 2-state fit is consistent 
with that from a single state fit to data at r > 15 (blue dotted line). 



far enough in r), one can make a fit to a larger range of r by including 
the contribution of higher (radial) states in the sum in Eq. 18.1. This is 
illustrated in Fig. 24 for a two state fit which involves four free parameters, 
two amplitudes and two masses. The amplitude and mass of the excited 
state are less well constrained as they are trying to mimic the effect of a 
number of states, consequently such fits may show significant dependence 
on the range of t selected. One can refine the excited state parameters 
by iteratively improving single state fits to T{t) — Tq{t) at short t and 
to the ground state Tq{t) at large r until the whole range is covered. A 
better approach to gain confidence in ground state parameters is to make 
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simultaneous fits to a matrix of correlation fimctions created with different 
source and sink operators with the same quantum numbers. For a matrix 
with m (n) distinct sources (sinks) the simultaneous fit is more constraining 
as it involves m + n amplitudes but a single common mass. 

To summarize, we mostly extract M from fits to an intermediate range 
of T, the "plateau" region, where Afe// is essentially independent of r. 
At small T the influence of excited states is large, while at large r the 
statistical noise overwhelms the signal. Accuracy of the result depends on 
the presence of a plateau extending over a large range of r that averages a 
few cycles of correlated fluctuations. The data in Figs. 22, 23, and 24 show 
that, with 170 lattices of size 32^ x 64 at 1/a = 2.3 GeV, this is true for the 
pion but marginal for other states like the vector mesons and baryons. Our 
present estimate is that roughly 1000 independent lattices are necessary to 
get a reliable signal that leads to 0.1 — 0.4% accuracy in baryon masses 
with quark masses in the range — ms/4 respectively. 

Assuming that the above steps have been carried out and we have reliable 
estimates for hadron masses as a function of quark masses, a, L, and the 
number n/ of dynamical flavors, we are finally in a position to compare 
them to experimental data. 

18.1. Status of the Light Meson and Baryon Spectrum 

The first test of QCD is to reproduce the masses of pseudoscalar mesons 
{VS = {tTjKjT]}), vector mesons (V = {p,K*,u},(p}), octet baryons 
(^8 = {iV,S,S,A}), and decuplet baryons (Bw = {A, E*, S*, f]}) from 
simulations with four input parameters mu,md,ma, and as- I will assume 
that the lattices are large enough that finite size effects can be neglected. 
The remaining important issues are 

— Simulations have mainly been done in the quenched approximation. 
Rather than testing QCD, these calculations quantify the accuracy of 
QQCD. 

— Isospin breaking effects and electromagnetic contributions are neglected. 

— Simulations have been done with unphysical values of quark masses, typi- 
cally in the range [2rn,s — ms/4:]. To get physical results the data are extrap- 
olated to m = [mu + md) /2 using functional forms predicted by (quenched) 
chiral perturbation theory. 

The first step in the analyses is to extrapolate the data to physical values 
of m and m^. The simplest Ansdtze are the lowest order chiral expressions 

M^ga^ = a^(mi+m2)/2-t-..., 
Mya = Av + Bv a{mi + TO2)/2 + . . . , 
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Mioa 



An + 4F ma + 2{F - D) rusa + . . . , 

An + 4(F - ^) ma + 2{F + y ) m,a + . . . , 

Aio + -Bio a{mi + m2 + ms)/^ + . . . . 



(18.3) 



From these fits we determine the constants ^v,^Ar,^io and B-ps, By, 
F, D, BiQ. Then using a suitable set, for example Ay, B-pg, By or Aiq, B-pg, 
Bio, we fix the three inputs parameters 1/a, m, mg. Unless otherwise stated 
the following discussion will assume that m, a, and rus have been fixed 
using Mj^jMp, and (or equivalently Mk*) respectively. The question 
then is - how well do the masses of the remaining states compare with 
experimental numbers? 

This simple analysis to test QCD (or to quantify the quenched approx- 
imation) can fall short in two ways. First, the data could show significant 
deviations from the linear relations Eq. 18.3. In that case the analysis 
becomes much more elaborate due to the increase in the number of free pa- 
rameters; the possible leading non-linear terms include O(m^) corrections 
and chiral logarithms. To disentangle these requires very precise data at 
a number of values of niq in the range ms/4 — 3ms, however, such a pos- 
sibility of physics richer than that predicted by the "naive" quark model 
is exciting. The second possible failing is that different sets of states used 
to fix a.m. m,, give different results. There could be three reasons for this 
failure: (i) The coefficients extracted from simulations at the significantly 
heavier values of ruq change as the quark masses are decreased; (ii) the dif- 
ferences are artifacts of discretization errors and vanish in the limit a = 0; 
(3) quenching errors. The current data show evidence for non-linearities 
in fits in rUq [156,160]. Also, the a,m,,ms extracted from different sets are 
significantly different, and all three causes listed above contribute. As a 
result the focus of current analyses has been to first disentangle and quan- 
tify the contributions of the various artifacts. I will briefly illustrate these 
effects. 

Let me start with evidence for non-linear terms in the chiral fits. Fig 25 
is an example of the behavior of the vector meson mass as a function of 
the quark mass [156]. The plot makes it clear that in order to quantify 
the deviations one needs very precise data over a large range of rUg. With 
a statistical accuracy of 0.5%, a linear fit works for 0.02 < ma < 0.04 
corresponding to 50 < m < 100 MeV, but not for the full range. The 
leading chiral corrections are proportional to m^^'^ from chiral loops, and 
m^. Including either term allows fits over the full range, and the data 
are not precise enough to distinguish between them. The effect of these 
non-linearities seems to be small, for example the change in Mp, after 




extrapolation to m, is < 1% on including either of the higher order terms 
in the fit. 

The first test of the quenched approximation is whether it can reproduce 
the hyperfine splittings in the meson sector. Figure 26 shows some current 
data. The only input quantity is the scale 1/a which is used to convert 
lattice masses to physical units. If M^r and Mp are used to set m and 1/a 
respectively, then the curves, by construction, have to pass through the 
experimental point at M-pg = Mt^. Different choices for the scale setting 
quantity will change 1/a and shift the curves as a whole. On the other 
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Fig. 26. Hyperfine splittings in the meson sector. Quenched data from (i) LANL collabo- 
ration with Wilson fermions (blue octagons) [156], (ii) LANL collaboration with tadpole 
improved clover fermions (red crosses) [162], and (iii) CP-PACS collaboration Wilson 
runs at /3 = 5.9, 6.1, 6.25, 6.47 (green fancy crosses, pluses, diamonds, and squares) [161]. 
The experimental points are shown by symbol burst. 



hand, the slope at any given point is independent of the lattice scale, and 
is therefore the more robust lattice prediction. For this reason the quantity 
J = MydMv /dMpg is commonly used as a test of the hyperfine splittings. 
Lattice data (for a review see [161]) gives values for J that are ~ 20% too 
small. Another manifestation of this discrepancy will appear in Section 20 
where I show that the estimates of to^ fixed using Mk versus Mjf. (or M^) 
differ by ~ 15%. 
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The data shown in Fig. 26 allow two comparisons. First, the octagons 
(blue points) and crosses (red) are LANL data from the same 170 quenched 
32^ X 64 lattices at /3 = 6.0 with the Wilson and tadpole improved clover 
Dirac actions [156,162]. One docs not find any significant improvement with 
clover fermions. Second, the fancy crosses, pluses, diamonds, and squares 
(green points) are Wilson fermion data from CP-PACS collaboration at /? = 
5.9, 6.1, 6.25, 6.47 [161]. The data at the four /3 values arc consistent. Both 
comparisons imply that the smaller splittings are not due to discretization 
errors. 

The data in Fig. 26 also show a difference between the CP-PACS and 
LANL Wilson fermion results at quark masses heavier than m^. This is 
surprising since in this regime both collaborations have reliable measure- 
ments of the vector and pscudoscalar masses. The culprit is the lattice 
scale. It turns out that the estimate of 1/a from LANL data is 3% larger 
compared to the value obtained by an interpolation of the CP-PACS data. 
Adjusting the LANL scale shifts their data as a whole and makes all the 
results overlap. 

The conclusion, from this and other quenched data, is that the splittings 
come out too small. Even though the size of discretization errors is, in my 
opinion, not fully resolved, the deviations are thought to be mostly due to 
the use of the quenched approximation. 

The most extensive discussion of the mass-splittings in the baryon octet 
and decuplet using the Wilson action is given in [156]. Note that since this 
calculation was done at only one 1/a « 2.3 GeV, therefore discretization 
errors arc not resolved. There arc two important observations, (i) The mass 
splittings in the baryon octet show significant non-linearities, for example 
the dependence of Ms — M\ on the quark mass is shown in Fig. 27. The 
non-linear terms significantly increase the splittings, and extrapolations 
to m give results roughly consistent with experimental values, (ii) The 
decuplet baryon masses show a behavior linear in mi -|- m2 -|- ms, and the 
splittings come out too small for either choice ms{MK) or m,s{M^). 

The state-of-the-art Wilson fermion data have been obtained by the CP- 
PACS collaboration. Data have been obtained at four values of the lattice 
spacing so a reliable extrapolation to a = can be made. Unfortunately, 
the data, in particular the baryon mass-splittings have not yet been fully 
analyzed so some of the observations presented here should be considered 
prehminary. Fig. 28, taken from Yoshie's review at LATTICE 97 [161], is 
CP-PACS's version of the chiral fits. They find that in the range mg/A—ms 
only the nucleon and the A particles show any significant curvature. Even 
though data at more values of the quark mass are required to resolve the 
form of the leading higher order correction, the conclusion of their analyses 
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Fig. 27. Quadratic fit to Ms — MA/{rn — ms), including baryons composed of completely 
non-degenerate quarks [156]. The value extrapolated to the physical point is shown by 
the burst symbol at the extreme left. The linear relations shown in Eq. 18.3 would 
predict a constant value. 



is that the negative curvature found in their fits to nuclcon data is suffi- 
ciently large to change the quenched estimate of /Mp to almost perfect 
agreement with the experimental value ~ 1.22. (This negative curvature 
was not seen/resolved in calculations reported before 1997, consequently a 
higher value in the range 1.3 — 1.4 was usually reported.) Their extrapola- 
tion to a = of results with different chiral fits is shown in Fig. 29. 

Second, they confirm that the decuplet masses are linear in the quark 
masses, and the splittings, even after extrapolation to a = 0, come out too 
small with either tUs^Mk) or TOs(M0). Their full analyses on splittings in 
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Fig. 28. (A) Chiral fit to M„^cleon data by the CP-PACS collaboration using the Wilson 
gauge and fcrmion actions. The data show negative curvature, though more points are 
needed to resolve between the cubic fit M^v = co + ciM-^ + C2M^ + c^M^, and xPT 
(ci = 0), QxPT (c3 = 0), xPT (fix) (ci = -0.5.S) fits. (B) Continuum extrapolation of 
the nucleon mass obtained using the different chiral expressions are shown in (A). The 
figures are reproduced from [161]. 



the baryon octet will appear soon. 

A final summary of the quenched results taken from [161] using the 
Wilson action is given in Fig. 30. This figure compares the recent CP- 
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Fig. 29. Continuum extrapolation of the nucleon, A, A, and Q masses obtained by the 
CP-PACS collaboration. The scale a is set by Mp. The previous best results by the 
GFll collaboration (unfilled symbols) are also shown. Experimental values are shown 
by the symbol diamond. Figure reproduced from [161] . 



PACS results extrapolated to a = with the previous best calculation by 
the GFll collaboration. The main lesson to be learned from this plot is 
that there does not seem to be a value of the strange quark mass that 
allows a matching of all the meson and baryon states. I shall come back 
to this point when discussing the extraction of quark masses. The obvious 
question is - is this deviation an artifact of the quenched approximation? 
My view is that we still do not have full control over the quenched data, and 
therefore cannot point to quenching as the sole culprit. For example, the 
definition of ruq itself has corrections of 0{mqa), incorporating which would 
change the nature of the chiral fits. To fully understand the discretization 
errors we need to wait for similarly good quality data with the staggered 
and 0{a) improved clover formulations before claiming definitive quenched 
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Fig. 30. State-of-the-art results from the CP-PACS collaboration for the meson and 
baryon masses after o = extrapolation and using Wilson fermions. The scale a is set 
using Mp and m using M„ . Results are shown for two ways of setting lUs - Mk and 
and it is clear that neither give correct splittings in the baryon octet and decuplet. 
The data are reproduced from [161]. 
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values. 

As of this writing dynamical simulations have begun in earnest, but 
precise quantitative results at a number of values of a and with different 
actions are lacking. Thus, I feel it is too early to draw any conclusions 
from them. For a status report on n/ = 2 results see the recent reviews by 
Yoshie and Gottlieb [161,160]. 

18.2. Decay Constants 

In addition to the extraction of masses from the rate of exponential de- 
cay of the 2-point correlation function, the amplitudes {|n|O|0) give decay 
constants. For example, the axial current matrix element (tt] ^^^^(a;) |0) = 
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if-nPfi cxp{ip ■ x) gives the pion decay constant. Some of the recent resuhs 
for decay constants have been reviewed by Prof. Martinelh at this school, 
and other recent reviews can be found in [163,164]. The technical points I 
would like to explain here are the reasons why the errors in decay constants 
are much larger than those in the determination of the masses. First, from 
Eq. 18.1 one notes that for a 1% error in the determination of the mass, the 
fractional error in the decay amplitude is w exp(O.OlmT) — 1 w 0.005mr. 
Since fits typically have mr « 5 for 1/a ~ 2 GeV lattices, the error in 
the decay constant is > 3%. Second, for any given configuration, the 
correlation function is highly correlated in t. Plotting log(r(T)) for each 
individual configuration one notes that the dominant fluctuation is in the 
overall normalization, i.e., the amplitude. A very rapid growth in these 
fluctuations with decreasing is, in fact, the signal for exceptional con- 
figuration as discussed in Section 10. These two sources of error combine 
to produce a statistical error in decay constants that is much larger than 
in the extraction of Me// . To this one has to add a systematic uncertainty 
coming from the renormalization constants for the local currents. 

18.3. Lattice Calculations of the Glueball Spectrum 

The existence of glueballs, hadronic states composed mainly of gluons, is 
a major untested prediction of QCD. Assuming QCD is the correct the- 
ory, we can determine the quantum numbers of a vast number of possible 
glueball states, but cannot yet calculate their masses, or their mixing with 
qq, qqqq, . . . states, or understand in detail the production and decay mech- 
anisms. Various models like bag models, flux tube model, sum rules etc. 
give estimates of masses that differ significantly from each other and the 
reliability of these calculations is poor. A somewhat dated summary of the 
theoretical expectations, compiled by Burnett and Sharpe, can be found in 
[165]. The lack of knowledge of properties of glueballs does not help the 
experimentalists who have to isolate these states from the myriad of meson 
states in the 1 — 2.5 GeV region. Clearly, LQCD calculations can play a 
vital role by calculating the spectrum and the mixing with qq mesons. 

The first goal of LQCD has been to calculate the quenched spectrum 
in a world in which the mixing with quark states is turned off. I shall 
concentrate on just the O"*""*" and 2++ states as these have been studied 
most extensively. The glueball operators O are constructed out of Wilson 
loops. Consider a n x m Wilson loop W'^y"' lying in the xy plane. Then 
the rotationally invariant sum W"jj"* + Wy^^ + W"^;"* projects onto 0++ 
state, while W^^f^ — Wy^"^ and its permutations project onto 2++. The 
signal in the 2-point correlator can be improved by taking a suitable linear 
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combination of different sized loops that maximize the overlap with physical 
glueball states. One such strategy is discussed in [166]. Unfortunately, the 
statistical signal in glueball correlation functions is much weaker than that 
for meson and baryon correlation functions. Even with a sample of O(IO^) 
independent configurations, the signal does not extend beyond r ~ 10 or 
Mr ~ 6 and one barely has a "plateau" in the effective mass plot. The 
reason for this poor signal in glueball correlation functions is that they 
measure the fluctuations in the vacuum corresponding to the creation and 
propagation of a virtual glueball, whereas in meson and baryon correlators 
the quarks in the propagators arc added as external sources. 

The more recent, high statistics results for Mq++ and M2++ states have 
been taken from Refs. [166-169]. The data and the extrapolation to a = 
are shown in Fig. 31. In this plot I have chosen to set the lattice scale 
using i?o = l-lS/^/a, where a is the string tension [170]. The figure 
makes it clear that the lattice data from the various collaborations are 
completely consistent, though the statistical quality of the GFll data 
[169] is much higher. A fit in a^, the leading discretization error, gives 
M0++R0 = 4.36(9) and M2++R0 = 6.13(16) for the continuum limit of the 
quenched theory. To obtain masses in physical units one needs to choose 
a value for Rq or equivalently a. This brings us back to the overall scale 
uncertainty inherent in all quenched calculations. There are two questions 
relevant to this analysis (i) what value should one choose for cr, and (ii) 
by how much would the results change if one had used a different quan- 
tity like Mp or Aqcd to set the scale? The two most quoted results are 
the following. The UKQCD-WUPPERTAL collaboration use l/Ra = 373 
MeV based on choosing ^ = 440 MeV and get Mo++ = 1626(34) and 
Mo++ = 2286(60) MeV [170]. The updated GFll collaboration estimates 
are Mo++ = 1648(58) and M2++ = 2273(116) GeV using A^^^ = 235(6) 
MeV which is based on Mp to set the scale [171]. Thus, the estimates 
of the two groups are consistent and agree with the fit shown in Fig. 31. 
Therefore, in my view, based on these fits, the current combined lattice 
estimate is Mo++ = 1640(30)(160) MeV and M2++ = 2280(60)(230) MeV, 
where I have added a second systematic error of 10% to take into account 
the scale uncertainty inherent to quenched simulations. 

18.4- Improved action Results on Anisotropic lattices 

To alleviate the problem that the signal in glueball 2-point correlation func- 
tions falls rapidly, Peardon and Morningstar have carried out simulations 

on anisotropic hypcrcubic lattices, at << a,, [172]. The advantage of this 
trick is the following. Assume that there exists an approximate plateau in 
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Fig. 31. Extrapolation of 0+"*" and 2++ glueball mass data to a = 0. The data are from 
UKQCD-WUPPERTAL collaboration (crosses), GFll collaboration (diamond), and Los 
Alamos collaboration (squares). The lattice scale is set using Rg which is related to the 
string tension by I/-R0 = v'<'^/l-18- 



Mef f (r) between r = 2 — 4 in an isotropic lattice simulation before the sig- 
nal degrades to noise. Such a plateau is hard to infer from three time-slices. 
On the other hand for ag/at — 5 lattices, one has roughly 11 times-slices to 
play with, thereby greatly increasing the confidence in the estimate. Using 
anisotropic lattices does introduce an additional parameter into the analy- 
ses, as /(It or the velocity of light. This they fix non-perturbatively using 
the measured potential in spatial and temporal directions. 
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Fig. 32. Glueball mass estimates in terms of Ro against the lattice spacing (as/Ro)^. 
Results from the ^ = 5 simulations for the lattice irreps A++ (0++), E++ , T++ (2++), 
and (!"' ) are labeled o,n,0, and A (blue), respectively. The corresponding 

symbols in red indicate the results from the 5 = 3 simulations. Data from Wilson action 
simulations taken from Refs. [167—169] are shown using green crosses. The dashed, 
dotted, and dash-dotted curves indicate extrapolations to the continuum limit obtained 
by fitting to the ^ = 3 data, the § = 5 data, and all data, respectively. The solid line 
indicates the extrapolation of the Wilson action data to the continuum limit. 



The second improvement they make is to tune the action. They tran- 
scribe the tadpole improved Liischer-Weisz glue action, Eq. 12.14, on to 
anisotropic lattices and include the plaquette in both the fundamental and 
adjoint coupling space. The motivation for the latter is to avoid the phase 
transition in the fundamental-adjoint space as discussed in Section 15. 
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Their new data [139] confirms that including a negative adjoint coupling 
does improve the scaling in Mq++. 

Their results, prior to those including a negative adjoint coupling, are 
shown in Fig. 32. The advantage of improved action and anisotropic lattices 
is clearly visible. First, the calculations can be performed on much coarser 
lattices. Second, the discretization errors are much smaller. Lastly, the 
values after extrapolation to a = arc consistent with conventional meth- 
ods. Their estimates M0++R0 = 3.98(15) and M2++R0 = 5.85(2), when 
combined with their estimate i?o = 410(20) MeV, are consistent with the 
values from conventional simulations given in the previous section. Thus, 
their results provide another confirmation of the quenched estimates. 

Let me end this discussion of gluey states with pointers to other interest- 
ing calculations. Results for other spin-parity glueball states can be found 
in [168,172,173]; and for hybrid states in [173-176]. 

18.5. Mixing between scalar quarkonia and glueballs 

The phenomenologically important question is whether the lattice es- 
timates can help pin down which of the two experimental candidates, 
/o(1500) or /o(1710), is the scalar glueball. The answer, using just the 
above quenched estimates, is no. The reason, in addition to the fact that 
the above estimates overlap with both candidates, is that effects of mixing 
with scalar quarkonia have not been included. Lee and Wcingarten have 
initiated a study of mixing by calculating the two-pint functions (»S(t)<S(0)) 
and (G'(r)5(0)) where S is the scalar density {uu + dd) / \/2 or ss and G is 
the 0+"'" glueball operator [171]. The statistical quality of these very hard 
calculations is not yet good. Lee and Weingarten, therefore, diagonalize 
an approximate 3x3 matrix that is obtained mainly from a phcnomcno- 
logical analyses. They then find that mixing raises Mq++ by ~ 64 MeV, 
i.e. M(,++ changes from 1640 1704 MeV. Also, the associated wave- 
functions of the mixed states explain some of the unexpected features of 
decays, like why the KK decay mode of /o(1500) is suppressed and why 
r(J/V ^ 7/o(1710)) > r(J/V ^ 7/o(1390)) > T{J/ij ^ 7/o(1500)). 
Consequently, Lee and Weingarten predict that /o(1710) is dominantly 
(~ 75%) a giuebah. 

In obtaining these predictions, they neglected mixing with other (multi- 
particle) states and decay widths in the mixing matrix and quenching errors 
in the lattice calculations. These, Lee and Weingarten claim on the basis 
of phenomenological arguments, are small and do not change their central 
conclusion. Hopefully, future LQCD calculations will provide a much better 
quantitative validation of their prediction. 
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18.6. Mass Inequalities 

In addition to getting hard numbers, LQCD (actually Euclidean field the- 
ory) has also prompted the derivation of rigorous mass- inequalities. The 
first such was derived by Weingarten showing that the pion is the lightest 
meson state of non-zero isospin in QCD [177]. Consider the 2-point corre- 
lation function for mesons, i.e. O = uTd where T can be any one of the 16 
elements of the Clifford algebra. Then 

rr(r) = (0|(iru(r,y) urrf(0,a;)|0) 

= {TvTj5S''J{T,y;0,x)^5rSUr,y;0,x)). (18.4) 

For the special case of the pion, = 1, consequently Tt^{t) has the max- 
imum value on each configuration. Furthermore, for to„ — nid it is the 
absolute square of the propagator. The other spin-parity states involve 7 
matrices. Since these have both plus and minus signs, the various spinor 
terms in the trace have cancellations. The second property needed to trans- 
late this property of correlation functions into an inequality on masses is 
that the integration measure for QCD is positive, i.e. each configuration 
contributes with the same sign. Thus, r^(T) > rr(T) for all r. Since, 
r(r) ~ e^"'^, the inequality implies that the pion is the lightest meson. 

Concurrently, Witten and Nussinov derived further inequalities. Exam- 
ples of these are 

— The electromagnetic mass shift in the pion is positive [179]. Specifi- 
cally what is shown is that {V^{k)V^{-k) - (A;)^^ (-/c)) > for any 
Euclidean k. 

— For mesons with valence flavors A and B the inequality 2M^g > 
M^-j + Mg-g holds in the limit that annihilation into gluonic states can be 
neglected in the flavor singlet correlators [179,178]. 

Recently West has applied/extended these methods to glueballs [181]. 
He shows that 0++ is the lightest glueball state. 

19. The strong coupling constant 

Spectroscopy studies, in addition to testing QCD, allow the determination 
of some of the fundamental parameters of the SM, i.e, the coupling and 
the quark masses. To determine ag, one needs to know not only its strength 
but also the energy scale Q at which the coupling is measured. The calcu- 
lation of Us, in general, proceeds as follows. One chooses a quantity (O) 
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that can be measured very accurately, either in high energy experiments or 
via non-perturbative methods on the lattice, and for which a perturbative 
expansion is known, i.e., 



where the An are the perturbative coefficients. Inverting this relation gives 



One promising lattice approach for calculating has been pioneered by 
the NRQCD collaboration [110]. They examine several short-distance (O), 
like the expectation value of small Wilson loops, that can be calculated very 

accurately on the lattice and for which perturbative expansions are known 
to at least O(a^). The scale Q characterizing these lattice observable is 
typically of the form Q = C/a where a is the lattice spacing and C is a 
constant that needs to be fixed. They determine C relevant to (O) by 
evaluating the mean momentum flow in the perturbative expressions for 
these quantities. The scale 1/a is fixed using level splittings in quarkonia. 
In particidar they use the spin-averaged S-P and the 1S-2S splittings in the 
Upsilon system. The advantage of using these splittings is that they show 
negligible sensitivity to the precise tuning of the bottom quark mass, and 
the corrections due to light quark loops are expected to be small. 

Finally, to convert from a,, in the lattice scheme to a ^g, they use the 

two loop relation, and include reasonable estimates of the term in the 
error analysis. Their final estimate 



includes quenching uncertainties. This remarkably precise result is consis- 
tent with the world average of experimental determinations, = 0.118(3), 
as shown in Fig. 33. 

An alternate method for measuring ag has been carried out by the AL- 
PHA collaboration [182]. In their method, called the Schrodinger func- 
tional approach, is defined by measuring the response of the lattice 
system to an external field. They fix the lattice scale a using the force 
at some fixed distance (i?o = 0.5 fermi). They evolve their coupling to a 
very high scale using a non-perturbative renormalization group stepping 
procedure. At this high scale they define A.^^ using the two loop relation 
and include the three loop uncertainty, which is small, in error estimates. 
They have completed the calculation for the quenched theory with the re- 
sult Ag^ = 251(21) MeV, and are now doing the n/ = 2 simulations [182]. 
The corresponding value Q;-g-^(90 GeV) = 0.0818(10) is lower than in na- 
ture. This difference is attributed to having neglected sea quark effects. 




(19.1) 



n 



asiQ). 



c^MsiMz) = 0.1174(24) 



(19.2) 
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A^^^ (in MS scheme, in GeV) 
0.04 0.1 0.2 0.5 0.8 



Deep inelastic 
GLS sum rule 
Deep inelastic HERA 
Jets at HERA 
QQ lattice 
X decay 
T decay 
Fragmentation 
Rie+e-) 



e+e~ event shape 
Z width* 
AVERAGE 



J I I I I I I I I I I I L. 



0.10 0.12 0.14 0.16 

Fig. 33. Summary of the values of as{Mz) and A^®' (the QCD scale for 5 active 
quark flavors) extracted from various processes. These are ordered from top to bottom 
by increasing energy scale of measurements. This figure is reproduced from the 1996 
Review of Particle Physics, updated to show the latest lattice result, aj^{Mz) = 
0.1174 it 0.0024. It is a remarkable confirmation of QCD that the Os extracted from 
processes at different scales has the same value. 



A similar low value is obtained in the NRQCD analyses if no corrections 

arc made for the quenched approximation. For details of the Schrodinger 
functional method see the lectures by Prof. Liischer at this school. 
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20. Quark Masses 

Quarks are not asymptotic states of nature, so their masses cannot be 
measured in experiments. They are generally defined in some pcrturbative 
scheme like MS and at some scale /x. Their values remain amongst the 
least well known of the Standard Model parameters. In this section I will 
discuss their extraction from lattice data. 

The extraction of the up, down, and strange quark masses from the 
hadron spectrum is done by inverting Eqs. 18.3 once the coefficients 
etc. are determined from fits to the lattice data. As I have mentioned 
before, current lattice simulations neglect iso-spin breaking effects, so we 
can only determine m = {mu + md)/2 and rus- I will discuss these first and 
then briefly mention the extraction of the heavy quark masses rric and m&. 

In the quenched approximation, the most straightforward way to es- 
timate m is to extrapolate M^/X"^ to the physical value. Here X is a 
dimcnsionful quantity like /jr or Mp or Mn or Ma that does not depend 
on nis, and de facto sets the lattice scale 1/a. Of the above quantities, 
using Mp is the most controversial as it decays via strong interactions and 
has a large width, 170 MeV. On the other hand the lattice data for Mp is 
most readily available and has good signal. Also, in the quenched approx- 
imation the p does not decay. Using Mp to set a in QQCD versus some 
other quantity can therefore give an overall scale uncertainty. The spread 
in 1/a extracted from different quantities is found to be « 10%, and this 
uncertainty, inherent to the quenched approximation, will be treated as an 
overall systematic error. Similarly, one can set by extrapolating the 
ratio Ms/X, where Ms is any state that depends on nis like Mk, or Mk', 
or M^ or any of the baryons containing a strange quark. The world data 
show that different choices of A4s give rise to a w 20% spread in mg. It is 
important to note that if one assumes the relation M^^ = B-psmq (which 
LQCD has not validated as true for the whole range [m, m^]), and uses Mk 
to fix ms, then the ratio ms/m has to be 25.9 if m is set using M^. 

In the full theory X depends on m^, and similarly Ms depends on m, 
through sea quark effects. Thus, m, a, and nis have to be fixed by solving 
the coupled set of equations. 

Simulations of the hadron spectrum have been done by a number of col- 
laborations around the world, and there exists a large amount of data for 
the masses of hadrons as a function of the quark masses and lattice spacing. 
(Almost all lattice collaborations, for example the APE, Columbia Univer- 
sity, CP-PACS, Fermilab, GFll, HEMCGC, JLQCD, Los Alamos, SCRI, 
staggered, UKQCD, and the Wuppertal University collaborations, have 
contributed to these calculations [183].) In early analyses the estimates of 
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quark masses varied by almost a factor of two between different formula- 
tions of fermions, and the pattern of the discretization and quenching errors 
was not obvious from the work of any one collaboration. Clarification of 
this issue has been obtained by putting together all the world data for 
Wilson, clover, and staggered fermions [183,184]. The final picture for m, 
based on the most recent, high statistics, and large lattices data is shown 
in Fig. 34 for the three fermion formulations. This figure shows clearly 
that the factor of ~ 2 difference at 1/a ~ 2 GeV shrinks to ~ 15% after 
extrapolation to o = 0, i.e. most of the difference can be explained as due 
to discretization errors. 

The behavior of versus a is very similar. A summary of the quenched 
results in MeV at scale /U = 2 GeV, based on an analysis of the 1997 world 
data [185], is 





Wilson 


TI Clover 


Staggered 




4.1(1) 


3.8(1) 


3.5(1) 


ms{MK) 


107(2) 


99(3) 


91(2) 




139(11) 


117(8) 


109(5) . 



These estimates were obtained using the following Ansatze. 

— Chiral fits were made keeping only the lowest order terms shown in 
Eqs. 18.3. 

— M^/Mp is used to fix m. 

— Results for rus are given for two choices, M^/Mp and M^/Mp. (Using 

Mk* /Mp gives values almost identical to those with M^/Mp. This implies 
that for the choice ms{M^) one could equivalently have set the scale using 
Mk' instead of Mp.) 

— Since linear chiral fits have been used, ms(Af^) = 25.9m. 

— The matching factor, Z^, between lattice and MS schemes is the tadpole 
improved 1-loop result. 

— The extrapolation to a = is done using the lowest order discretization 
error, i.e. 0{a) and ©(a^a) for Wilson and clover, and a quadratic in 
for staggered. 

The final difference in estimates between Wilson, tadpole improved 
clover (TI Clover), and staggered results is ~ 15% as shown in Fig. 34. 
This could be due to the neglected higher order discretization errors and/or 
due to the difference between non-perturbative and 1-loop estimates of Z^s. 
Similarly, the « 20% variation in rUs with the state used to extract it, Mk 
versus M^ (or equivalently Mk*) could be due to the quenched approxi- 
mation or partly an artifact of keeping only the lowest order discretization 
correction in the extrapolations. To disentangle these discretization and 
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Fig. 34. The behavior of mu + md as a function of the lattice spacing a for three 
discretization schemes - Wilson, clover, and staggered. Linear extrapolations to a = 
are shown for Wilson and clover formulations, while for staggered both 0(a?) and ©(a*) 
terms are incorporated in the fit. The theoretically expected form for discretization 
errors in the clover formulation is O(asa). The extrapolated value using this form is 
shown by the symbol burst. 

quenching errors requires precise unquenched data with roughly physical 
values for light quarks. 

For best estimates of quenched results I average the data and use the 
spread as the error. To these, I add a second systematic error of 10% due to 
the uncertainty in the determination of the scale 1/a, which partly reflects 
the uncertainty due to the quenched approximation. The results, in the 
MS scheme and evaluated at 2 GeV, are [185] 



TO = 3.8(4)(4) MeV 
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ms = 110(20)(ll) MeV. (20.1) 
20.1. Light Quark Masses from Ward Identities 

A second method for extracting quark masses is based on axial and vector 
current Ward identities 

ZAd^Af^i = (mi + m2)RZp'P , 

Zvd^V^ = (mi - m2)RZsS , (20.2) 

where 1,2 are the flavor indices and P and S are the lattice pseudoscalar 
and scalar densities. From these one can construct ratios of 2-point corre- 
lation functions 

(dMr) V) = (mi + m2)R^{V{T) V{0)) 

(dMr) S{0)) = (mi - m2)R^{S{T) 5(0)) , (20.3) 

which give the renormalized quark masses once the various renormalization 

constants Zi arc known. Since the Ward identities are operator identities, 
Eq. 20.3 holds for all t except at very short times where there is influence 
of contact terms. Thus it is not necessary for the signal to persist to 
asymptotic r as needed in spectroscopy method to isolate the lowest state. 

The two estimates of quark masses, from spectroscopy and from the 
Ward identities, should agree in the continuum limit. A comparison of the 
two data sets for the Wilson action and using 1-loop pcrturbativc Z's is 
shown in Fig. 35. The two estimates are significantly different at finite a, 
however, a linear extrapolation works for both sets of data, and the ex- 
trapolated values are in surprisingly good agreement. One could therefore 
attribute the differences at finite a to discretization errors. 

The recent non-perturbative calculations of the Z's [186] challenge this 
nice picture. Fig. 35 also shows the same data analyzed using these non- 
perturbative Z's. At fixed a the non-perturbative Z's increase both the 
hadron spectroscopy and WI results to roughly the same value. This sug- 
gests that most of the difference is due to the O(a^) corrections to the 1- 
loop Z's and not due to the 0(a) discretization errors as mentioned before. 
However, the non-perturbative calculation of Z's themselves have 0{a) er- 
rors that have not been resolved in the results of [186]. As a result we do 
not have a clear separation between the effects of the 0{a) discretization 
errors and the 0(0^) errors in these Z's. 

The bottom line is that, even though we have not resolved the two 
kinds of errors, and while both the value at fixed a and the slope in a 




Fig. 35. The behavior of m = (rriu +mj)/2 as a function of the lattice spacing a for Wil- 
son fermions. Data from the hadron spectroscopy and Ward Identity method are shown 
for two ways of estimating the renormalization constants, i.e. one loop perturbative and 
non-pert urbative. 



increases significantly on changing from 1-loop to non-perturbative Z's, 
the extrapolated value is roughly unchanged. 

To conclude, I believe that we have finally obtained values for light quark 
masses with reliable estimates of statistical and all systematic errors other 
than quenching. Calculations are already underway to significantly increase 
the number of data points with Uf = 2 and 4 flavors of dynamical quarks to 
reduce the quenching uncertainty. Also, new ways to calculate the renor- 
malization constants non-perturbatively are being developed. Therefore, I 
expect that in the next few years there will be a significant increase in the 
precision of both QQCD and full QCD estimates of light quark masses. 
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20.2. Schrodinger Functional Approach 

In the above discussion the term non-perturbative Z's referred to only the 
lattice scheme. In the required matching constants one still has to use per- 
turbative expressions for the renormalization constants in the continuum 
scheme and make a good estimate of the matching point. The ALPHA 
Collaboration has proposed a method for removing this last vestige of de- 
pendence on perturbation theory [187]. The details of their method have 
been discussed by Prof. Liischer at this school, so I will just summarize 
the key points. 

(i) Calculate the quark mass using the axial WI method at some infrared 
lattice scale. At the same time calculate Zp and Za non-perturbatively 
in the Schrodinger functional scheme. 

(ii) Evolve this lattice result to some very high scale using the step scaling 
function calculated non-perturbatively. 

(iii) At this high scale, where as is small, define the renormalization group 
independent mass 

m = limmfl(2/3off')-^°/''''' (20.4) 

This fh is scheme independent, thus subsequent conversion to a running 
mass in some scheme and at some desired scale can be done in the con- 
tinuum using the most accurate versions (4-loop) of the scale evolution 
equations [188]. 

To reiterate, the advantage of this method is that the lattice calcu- 
lations directly predict the scale and scheme invariant quantity fh. The 
calculations by the ALPHA Collaboration are still in progress and at this 
time estimates for fh have not been released. It will be interesting to see 
how they compare with conventional hadron spectroscopy or WI methods 
discussed above. My prejudice is that they will be very similar. 

20.3. Role of quark masses in CP violation 

A very opportune application of the precise determination of the light quark 

masses is in the Standard Model prediction of direct CP violation in kaon 
decays characterized by the ratio e'/e (see lectures by Prof. Buras and 
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[189]). The reason is that e'/e is, to a good approximation, proportional to 
l/(md + m,)2 [190], 

.A = A{c..M-+.Br][^i^]'}. (20.5, 

where all quantities are to be evaluated at the scale ruo = 1.3 GeV. Eq. 20.5 

highlights the dependence on the light quark masses and the bag param- 

1/2 3/2 

eters Bg and B^' which LQCD is also expected to provide. The coef- 
ficients Co, cq, and cg have been calculated to next-to- next- to-leading or- 
der, and with Mtop now known the main uncertainty in their extraction 
comes from the value of Aqcd- The quantity A = | V^;, 1 1 Vc6 1 sin 5 in the 
parameterization of Eq. 2.4, requires input of other SM parameters like 
fB,\Vub\,\Vcb\,BK, \/BBafBa- Using the central values quoted by Buras at 
scale rric = 1.3 GeV [190] one gets A = 1.29 x 10"^, cq = -1.4, cq = 7.9, 
cg = —4.0. Unfortunately, the errors in these quantities make the resulting 
theoretical prediction, e'/e = 3.6(3.4) x 10"'' [190], very loose. 

Current experimental estimates are 7.4(5.9) x 10"'' from Fermilab E731 
[191] and 23(7) x 10"* from CERN NA31 [192]. The new generation of 
experiments, Fermilab E832, CERN NA48, and DA$NE KLOE, will reduce 
the uncertainty to « 1 x 10"''. First results from these experiments (see 
lectures by L. Fayard) should be available in the next couple of years. Thus, 
it is very important to tighten the theoretical prediction. Clearly, lower 
estimates of quark masses suggested by LQCD would significantly increase 
the size of direct CP violation, and make its experimental determination 
and validation much easier. 

20.4- Heavy quark masses rric and rrib 

The masses of the heavy quarks, rric and mf,, have been calculated using the 
charmonium and upsilon spectroscopy. For heavy quarks there exist two 

different definitions of the quark mass the pole and MS masses. The 

pole mass is defined to be 771^^'°''=) = {Monia - Etinding)/'2, where Monia is 
taken from experimental data and Ebinding is calculated on the lattice. This 
calculation has been done using two approaches, NRQCD [193] and heavy 
quark effective theory (HQET) [194], and the results are consistent. The 
MS mass, which has the advantage of being free of renormalon ambiguities 
that are inherent in m^P°'-'^\ is obtained from rn^^"'*^ using perturbation 
theory. The present results are 

mf^(mf^) =4.15(5)(20) GeV (APE) 

= 4.16(15) GeV (NRQCD). (20.6) 
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These estimates are consistent with determinations from other non-lattice 

methods. 

The agreement is not as good for various lattice estimates of m^^ . 
Current estimates lie between 1.2 — 1.5 GeV. For example, three different 
methods investigated by the APE collaboration [195,186], Fermilab collab- 
oration [196,197], and by Bochkarev and Forcrand [198] who evaluate the 
correlation functions arising in QCD sum-rules on the lattice, give 



"ic (^c ) = 1.525(40)(125) GeV 
= 1.33(8) GeV 
= 1.22(5) GeV 



Fermilab 



B&F . 



APE 



(20.7) 



At this stage I do not consider the difference significant as the various 
systematic errors have not been resolved equally well in all the calculations. 
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Conclusions and Acknowledgements 

I hope that these lectures succeed in providing a pedagogical introduc- 
tion to Lattice QCD, and to the analyses of numerical data. Undoubtedly, 
some topics have been left out and many are sketchy. My purpose was 
to give a broad overview. I have concentrated on providing the basics and 
have illustrated the analyses of Monte Carlo data by the calculations of the 
hadron spectrum, strong coupling constant a^, and quark masses. Prof. 
Martinelli will, among other things, discuss the extraction of matrix ele- 
ments of weak Hamiltonian, while Prof. Liischer will provide a panoramic 
view of what future calculations will look like and the improvement possible 
with better discretization schemes. Together, I hope these set of lectures 
provide sufficient details and excitement to motivate some of you to probe 
deeper into the field. 

I express my gratitude to all the lecturers, students, and staff who made 
this Les Houches school fun and exciting. It, most certainly, was a mem- 
orable six weeks for me. In preparing these lectures I would like to thank 
G. Bali, U. Heller, C. Michael, M. Peardon, the CP-PACS collaboration, 
in particular S. Hashimoto and T. Yoshic. for providing unpublished data. 
Special thanks go to my collaborators T. Bhattacharya, D. Daniel, G. Kil- 
cup, A. Patel, and S. Sharpe from whom I have, over the years, learned so 
much. Finally, I would like to thank T. Bhattacharya, D. Daniel, W. Lee, 
and S. Sharpe for proof reading these notes and for their suggestions. 
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